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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗
DOE OSTI2018

Acceleration of tropical cyclogenesis by self-aggregation feedbacks

Idealized simulations of tropical moist convection have revealed that clouds can spontaneously clump together in a process called self-aggregation. This results in a state where a moist cloudy region with intense deep convection is surrounded by extremely dry subsiding air devoid of deep convection. Because of the idealized settings of the simulations where it was discovered, the relevance of self-aggregation to the real world is still debated. Here in this study, we show that self-aggregation feedbacks play a leadingorder role in the spontaneous genesis of tropical cyclones in cloudresolving simulations. Those feedbacks accelerate the cyclogenesis process by a factor of 2, and the feedbacks contributing to the cyclone formation show qualitative and quantitative agreement with the self-aggregation process. Once the cyclone is formed, wind-induced surface heat exchange (WISHE) effects dominate, although we find that self-aggregation feedbacks have a small but nonnegligible contribution to the maintenance of the mature cyclone. In conclusion, our results suggest that self-aggregation, and the framework developed for its study, can help shed more light into the physical processes leading to cyclogenesis and cyclone intensification. In particular, our results point out the importance of the longwave radiative cooling outside the cyclone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Early Experiences Porting the NAMD and VMD Molecular Simulation and Analysis Software to GPU-Accelerated OpenPOWER Platforms, In: ISC High Performance 2016: High Performance Computing

All-atom molecular dynamics simulations of biomolecules provide a powerful tool for exploring the structure and dynamics of large protein complexes within realistic cellular environments. Unfortunately, such simulations are extremely demanding in terms of their computational requirements, and they present many challenges in terms of preparation, simulation methodology, and analysis and visualization of results. We describe our early experiences porting the popular molecular dynamics simulation program NAMD and the simulation preparation, analysis, and visualization tool VMD to GPU-accelerated OpenPOWER hardware platforms. We report our experiences with compiler-provided autovectorization and compare with hand-coded vector intrinsics for the POWER8 CPU. We explore the performance benefits obtained from unique POWER8 architectural features such as 8-way SMT and its value for particular molecular modeling tasks. Finally, we evaluate the performance of several GPU-accelerated molecular modeling kernels and relate them to other hardware platforms.

Computer Science↗
DOE OSTI2015

A KINETIC TRANSPORT THEORY FOR PARTICLE ACCELERATION AND TRANSPORT IN REGIONS OF MULTIPLE CONTRACTING AND RECONNECTING INERTIAL-SCALE FLUX ROPES

Simulations of particle acceleration in turbulent plasma regions with multiple contracting and merging (reconnecting) magnetic islands emphasize the key role of temporary particle trapping in island structures for the efficient acceleration of particles to form hard power-law spectra. Statistical kinetic transport theories have been developed that capture the essential physics of particle acceleration in multi-island regions. The transport theory of Zank et al. is further developed by considering the acceleration effects of both the mean and the variance of the electric fields induced by the dynamics of multiple inertial-scale flux ropes. A focused transport equation is derived that includes new Fokker-Planck terms for particle scattering and stochastic acceleration due to the variance in multiple flux-rope magnetic fields, plasma flows, and reconnection electric fields. A Parker transport equation is also derived in which a new expression for momentum diffusion appears, combining stochastic acceleration by particle scattering in the mean multi-flux-rope electric fields with acceleration by the variance in these electric fields. Test particle acceleration is modeled analytically considering drift acceleration by the variance in the induced electric fields of flux ropes in the slow supersonic, radially expanding solar wind. Hard power-law spectra occur for sufficiently strong inertial-scale flux ropes with an index modified by adiabatic cooling, solar wind advection, and diffusive escape from flux ropes. Flux ropes might be sufficiently strong behind interplanetary shocks where the index of suprathermal ion power-law spectra observed in the supersonic solar wind can be reproduced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling

Kinetic Monte Carlo (KMC) simulation provides insights into catalytic reactions unobtainable with either experiments or mean-field microkinetic models. Sensitivity analysis of KMC models assesses the robustness of the predictions to parametric perturbations and identifies rate determining steps in a chemical reaction network. Stiffness in the chemical reaction network, a ubiquitous feature, demands lengthy run times for KMC models and renders efficient sensitivity analysis based on the likelihood ratio method unusable. We address the challenge of efficiently conducting KMC simulations and performing accurate sensitivity analysis in systems with unknown time scales by employing two acceleration techniques: rate constant rescaling and parallel processing. We develop statistical criteria that ensure sufficient sampling of non-equilibrium steady state conditions. Our approach provides the twofold benefit of accelerating the simulation itself and enabling likelihood ratio sensitivity analysis, which provides further speedup relative to finite difference sensitivity analysis. As a result, the likelihood ratio method can be applied to real chemistry. We apply our methodology to the water-gas shift reaction on Pt(111).

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

Phase Space Reconstruction from Accelerator Beam Measurements Using Neural Networks and Differentiable Simulations

Characterizing the phase space distribution of particle beams in accelerators is a central part of accelerator understanding and performance optimization. However, conventional reconstruction-based techniques either use simplifying assumptions or require specialized diagnostics to infer high-dimensional (> $2D$) beam properties. In this Letter, we introduce a general-purpose algorithm that combines neural networks with differentiable particle tracking to efficiently reconstruct high-dimensional phase space distributions without using specialized beam diagnostics or beam manipulations. Furthermore, we demonstrate that our algorithm accurately reconstructs detailed 4D phase space distributions with corresponding confidence intervals in both simulation and experiment using a single focusing quadrupole and diagnostic screen. This technique allows for the measurement of multiple correlated phase spaces simultaneously, which will enable simplified 6D phase space distribution reconstructions in the future.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Electron acceleration by magnetic islands in a dynamically evolved coronal current sheet

This work simulated the electron acceleration by magnetic islands in a drastically evolved solar coronal current sheet via the combined 2.5-dimensional (2.5D) resistive Magnetohydrodynamics (MHD) and guiding-center approximation test-particle methods. With high magnetic Reynolds number of 105, the long–thin current sheet is evolved into a chain of magnetic islands, growing in size and coalescing with each other, due to tearing instability. The acceleration of electrons is studied in one typical phase when several large magnetic islands are formed. The results show that the electrons with an initial Maxwell distribution evolve into a heavy-tailed distribution and more than 20% of the electrons can be accelerated higher than 200 keV within 0.1 second and some of them can even be energized up to MeV ranges. The most energetic electrons have a tendency to be around the outer regions of the magnetic islands or to be located in the small secondary magnetic islands. We find that the acceleration and spatial distributions of the energetic electrons is caused by the trapping effect of the magnetic islands and the distributions of the parallel electric field E{sub p}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Electron Acceleration during Macroscale Magnetic Reconnection

Here, the first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations, the spectra of energetic electrons take the form of power laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A strong guide field suppresses the production of nonthermal electrons by weakening the Fermi drive mechanism. For a weak guide field the total energy content of nonthermal electrons dominates that of the hot thermal electrons even though their number density remains small. Our results are benchmarked with the hard x-ray, radio, and extreme ultraviolet observations of the X8.2-class solar flare on September 10, 2017.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

CompF2: Theoretical Calculations and Simulation Topical Group Report

This report summarizes the work of the Computational Frontier topical group on theoretical calculations and simulation for Snowmass 2021. We discuss the challenges, potential solutions, and needs facing six diverse but related topical areas that span the subject of theoretical calculations and simulation in high energy physics (HEP): cosmic calculations, particle accelerator modeling, detector simulation, event generators, perturbative calculations, and lattice QCD (quantum chromodynamics). The challenges arise from the next generations of HEP experiments, which will include more complex instruments, provide larger data volumes, and perform more precise measurements. Calculations and simulations will need to keep up with these increased requirements. The other aspect of the challenge is the evolution of computing landscape away from general-purpose computing on CPUs and toward special-purpose accelerators and coprocessors such as GPUs and FPGAs. These newer devices can provide substantial improvements for certain categories of algorithms, at the expense of more specialized programming and memory and data access patterns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

SIMULATION OF ENERGETIC PARTICLE TRANSPORT AND ACCELERATION AT SHOCK WAVES IN A FOCUSED TRANSPORT MODEL: IMPLICATIONS FOR MIXED SOLAR PARTICLE EVENTS

We use numerical solutions of the focused transport equation obtained by an implicit stochastic differential equation scheme to study the evolution of the pitch-angle dependent distribution function of protons in the vicinity of shock waves. For a planar stationary parallel shock, the effects of anisotropic distribution functions, pitch-angle dependent spatial diffusion, and first-order Fermi acceleration at the shock are examined, including the timescales on which the energy spectrum approaches the predictions of diffusive shock acceleration theory. We then consider the case that a flare-accelerated population of ions is released close to the Sun simultaneously with a traveling interplanetary shock for which we assume a simplified geometry. We investigate the consequences of adiabatic focusing in the diverging magnetic field on the particle transport at the shock, and of the competing effects of acceleration at the shock and adiabatic energy losses in the expanding solar wind. We analyze the resulting intensities, anisotropies, and energy spectra as a function of time and find that our simulations can naturally reproduce the morphologies of so-called mixed particle events in which sometimes the prompt and sometimes the shock component is more prominent, by assuming parameter values which are typically observed for scattering mean free paths of ions in the inner heliosphere and energy spectra of the flare particles which are injected simultaneously with the release of the shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Direct observations of shape fluctuation in long-time atomistic simulations of metallic nanoclusters

Metallic nanoclusters are functional materials whose unique physical and chemical properties are sensitively controlled by their shapes and structures. An in-depth understanding of their morphological stability is therefore of crucial importance. It has been well documented by transmission electron microscopy (TEM) studies that metallic nanoclusters can interconvert between different isomers. However, the relevant mechanisms remain elusive because the timescales of such shape fluctuations are too short to be resolved experimentally and yet too long for conventional atomistic simulations. By employing massively parallel accelerated molecular dynamics (AMD) simulations reaching timescales of milliseconds, here we provide a clear description of the dynamical processes leading to shape fluctuation in metallic nanoclusters of platinum, silver, copper, and gold at different sizes and temperatures. In all these materials, direct transformations between face-centered-cubic and fivefold symmetric structures (decahedron or icosahedron) are observed away from the melting point. These transitions occur following either a slip-mediated twinning mechanism or a surface-reconstruction driven process. The identified pathways are shown to be reversible, allowing for genuine shape fluctuation processes that do not involve melting or other external factors.

36 MATERIALS SCIENCE↗
DOE OSTI2015

NONTHERMALLY DOMINATED ELECTRON ACCELERATION DURING MAGNETIC RECONNECTION IN A LOW-β PLASMA

By means of fully kinetic simulations, we investigate electron acceleration during magnetic reconnection in a nonrelativistic proton–electron plasma with conditions similar to solar corona and flares. We demonstrate that reconnection leads to a nonthermally dominated electron acceleration with a power-law energy distribution in the nonrelativistic low-β regime but not in the high-β regime, where β is the ratio of the plasma thermal pressure and the magnetic pressure. The accelerated electrons contain most of the dissipated magnetic energy in the low-β regime. A guiding-center current description is used to reveal the role of electron drift motions during the bulk nonthermal energization. We find that the main acceleration mechanism is a Fermi-type acceleration accomplished by the particle curvature drift motion along the electric field induced by the reconnection outflows. Although the acceleration mechanism is similar for different plasma β, low-β reconnection drives fast acceleration on Alfvénic timescales and develops power laws out of thermal distribution. The nonthermally dominated acceleration resulting from magnetic reconnection in low-β plasma may have strong implications for the highly efficient electron acceleration in solar flares and other astrophysical systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Laser–solid interaction studies enabled by the new capabilities of the iP2 BELLA PW beamline

The new capabilities of the short focal length, high intensity beamline, named iP2, at the BELLA Center will extend the reach of research in high energy density science, including accessing new regimes of high gradient ion acceleration and their applications. This 1 Hz system will provide an on-target peak intensity beyond 1021 W/cm2 with a temporal contrast ratio of <10−14 that will be enabled by the addition of an on-demand double plasma mirror setup. An overview of the beamline design and the main available diagnostics are presented in this paper as well as a selection of accessible research areas. As a demonstration of the iP2 beamline's capabilities, we present 3D particle-in-cell simulations of ion acceleration in the magnetic vortex acceleration regime. The simulations were performed with pure hydrogen targets and multi-species targets. Proton beams with energy up to 125 MeV and an approximately 12° full angle emission are observed as preplasma scale length and target tilt are varied. The number of accelerated protons is on the order of 109/MeV/sr for energies above 60 MeV.

Hakimi, Sahel (ORCID:000000022493564X)↗