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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Geometry Effects on Multipole Components and Beam Optics in High-Velocity Multi-Spoke Cavities

Velocity-of-light, multi-spoke cavities are being proposed to accelerate electrons in a compact light-source. There are strict requirements on the beam quality which require that the linac have only small non-uniformities in the accelerating field. Beam dynamics simulations have uncovered varying levels of focusing and defocusing in the proposed cavities, which is dependent on the geometry of the spoke in the vicinity of the beam path. Here we present results for the influence different spoke geometries have on the multipole components of the accelerating field and how these components, in turn, impact the simulated beam properties.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗
DOE OSTI2018

Structural Basis for Binding of Allosteric Drug Leads in the Adenosine A 1 Receptor

Despite intense interest in designing positive allosteric modulators (PAMs) as selective drugs of the adenosine A 1 receptor (A 1 AR), structural binding modes of the receptor PAMs remain unknown. Using the first X-ray structure of the A 1 AR, we have performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) technique to determine binding modes of the A 1 AR allosteric drug leads. Two prototypical PAMs, PD81723 and VCP171, were selected. Each PAM was initially placed at least 20 Å away from the receptor. Extensive GaMD simulations using the AMBER and NAMD simulation packages at different acceleration levels captured spontaneous binding of PAMs to the A 1 AR. The simulations allowed us to identify low-energy binding modes of the PAMs at an allosteric site formed by the receptor extracellular loop 2 (ECL2), which are highly consistent with mutagenesis experimental data. Furthermore, the PAMs stabilized agonist binding in the receptor. In the absence of PAMs at the ECL2 allosteric site, the agonist sampled a significantly larger conformational space and even dissociated from the A 1 AR alone. In summary, the GaMD simulations elucidated structural binding modes of the PAMs and provided important insights into allostery in the A 1 AR, which will greatly facilitate the receptor structure-based drug design.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

A Linac Simulation Code for Macro-Particles Tracking and Steering Algorithm Implementation

In this paper, a linac simulation code written in Fortran90 is presented and several simulation examples are given. This code is optimized to implement linac alignment and steering algorithms, and evaluate the accelerator errors such as RF phase and acceleration gradient, quadrupole and BPM misalignment. It can track a single particle or a bunch of particles through normal linear accelerator elements such as quadrupole, RF cavity, dipole corrector and drift space. One-to-one steering algorithm and a global alignment (steering) algorithm are implemented in this code.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Summary of Jefferson Lab LDRD on FFA@CEBAF Beam Dynamics Simulations

As Thomas Jefferson National Accelerator Facility (Jefferson Lab) looks toward the future, we are considering expanding our energy reach by using Fixed-Field Alternating Gradient (FFA) technology. Significant efforts have been made to design a hybrid accelerator which combines conventional recirculating electron LINAC design with permanent magnet-based FFA technology to increase the number of beam recirculations, and thus the energy. In an effort to further this progress, Jefferson Lab awarded a Laboratory Directed Research and Development (LDRD) grant to focus not on the design, but on detailed simulations of the designs created by the larger collaboration. This document will summarize the work performed during this LDRD, and direct the reader to other proceedings which describe elements of the work in greater detail.

Deitrick, K.↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
DOE OSTI2015

Particle acceleration and plasma dynamics during magnetic reconnection in the magnetically dominated regime

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional (3D) kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron–positron plasmas starting with a magnetically dominated, force-free current sheet (σ≡B{sup 2}/(4πn{sub e}m{sub e}c{sup 2})≫1). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra f∝(γ−1){sup −p} and approaches p = 1 for sufficiently large σ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to nonrelativistic regimes. In the 3D simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications for a variety of high-energy astrophysical problems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-T-635: Quantitative Study On Beam Flatness Variation with Beam Energy Change

Purpose: Beam flatness check has been proposed for beam energy check for photon beams with flattering filters. In this work, beam flatness change with beam energy was investigated quantitatively using the Monte Carlo method and its significance was compared with depth dose curve change. Methods: Monte Carlo simulations for a linear accelerator with flattering filter were performed with different initial electron energies for photon beams of 6MV and 10MV. Dose calculations in a water phantom were then perform with the phase space files obtained from the simulations. The beam flatness was calculated based on the dose profile at 10 cm depth for all the beams with different initial electron energies. The percentage depth dose (PDD) curves were also analyzed. The dose at 10cm depth (D10) and the ratio of the dose at 10cm and 20cm depth (D10/D20) and their change with the beam energy were calculated and compared with the beam flatness variation. Results: It was found that the beam flatness variation with beam energy change was more significant than the change of D10 and the ratio between D10 and D20 for both 6MV and 10MV beams. Half MeV difference on the initial electron beam energy brought in at least 20% variation on the beam flatness but only half percent change on the ratio of D10 and D20. The change of D10 or D20 alone is even less significant. Conclusion: The beam energy impact on PDD is less significant than that on the beam flatness. If the PDD is used for checking the beam energy, uncertainties of the measurement could possibly disguise its change. Beam flatness changes more significantly with beam energy and therefore it can be used for monitoring the energy change for photon beams with flattering filters. However, other factors which may affect the beam flatness should be watched as well.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1985

Simulation study of Type 2 counterstreaming electrons along auroral field lines

The production of counterstreaming electrons associated with parallel fields along auroral field lines is examined through the use of computer simulation. A 2 1/2-dimensional (two spatial and three velocity dimensions) electrostatic particle algorithm and auroral boundary conditions are used to set up a self-consistent V potential structure. The simulation produces signatures of counterstreaming electrons resembling those observed by the Dynamics Explorer 1 satellite. The main signatures are as follows: (1) the phase space contours of the electron distribution function are elongated along the V-parallel axis, and (2) the energy of electrons streaming in the upward direction is comparable to the energy of the accelerated electron beam. The simulation indicates that a portion of the accelerated electron beam is trapped by large amplitude electrostatic waves produced through the two-stream instability. Strong wave-particle interactions then thermalize the trapped electrons to produce suprathermal electrons streaming in the direction opposite to that of the accelerated electron beam. These results suggest a possible mechanism of producing counterstreaming electron fluxes through nonlinear processes of the two-stream instability.

Wagner, J. S.↗
NASA NTRS2013

Numerical Simulation of the Anomalous Transport of High-Energy Cosmic Rays in Galactic Superbubble

A continuous-time random-walk (CTRW) model to simulate the transport and acceleration of high-energy cosmic rays in galactic superbubbles has recently been put forward (Barghouty & Schnee 2102). The new model has been developed to simulate and highlight signatures of anomalous transport on particles' evolution and their spectra in a multi-shock, collective acceleration context. The superbubble is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. This work concentrates on the effects of the bubble's assumed astrophysical characteristics (cf. geometry and roughness) on the particles' spectra.

Barghouty, A. F.↗
DOE OSTI2023

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗
DOE OSTI2015

SU-E-T-461: Impact of the Radiation Sensitivity of Tumor and Geometric Design of the GRID Block On the Clinical Response of Grid Therapy

Purpose: The spatially fractionated radiation therapy (Grid therapy) is a novel technique that has been recently introduced for treatments of the patients with advanced bulky tumors. The purpose of this study is to investigate the effect of the radiation sensitivity of tumor and Grid block design on the clinical response of the Grid therapy. Methods: The Geant4 toolkit is used to simulate a Varian2100C linear accelerator with the block-based Grid. The dose profile has been utilized to calculate the Therapeutic Ratio (TR) and Equivalent Uniform Dose (EUD) for different types of tumors. These calculations were performed with respect to their radiation sensitivities using Linear Quadratic model (LQ). Clinical responses of few hundred patients from different publications that were treated with Grid therapy, have been used to validate correlation between tumor radiation sensitivity and clinical response of Grid therapy. Moreover, the influence of Grid design on clinical response of this radiation technique was investigated with using Monte Carlo simulation of a Grid block with different hole-diameters of 0.5 cm, 0.75 cm, 1.25 cm, and 1.5 cm. The Geometrical Sparing Factor (GSF) was considered to compare the results from different Grid designs. Results: The higher clinical response of Grid therapy for more radio-resistant tumors was confirmed in both theoretical assumptions and clinical outcomes. In addition, the results indicate that the Grid hole diameter have more effect on increasing the TR in comparison with the hole center-to-center distance effect. From this Monte Carlo study, a Grid design with the hole diameter between 1.25 cm and 1.5 cm with GSF value of 0.90 is recommended. Conclusion: Grid Therapy is more effective for more radio-resistant tumors.The optimum results of the Grid therapy depend on selection of the appropriate Grid design.This study is a first step toward clarifying some important aspects in a clinical response of spatially fractionated radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. II. Connection with γ-Rays

Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Investigation of longitudinal proton acceleration in exploded targets irradiated by intense short-pulse laser

It was recently shown that a promising way to accelerate protons in the forward direction to high energies is to use under-dense or near-critical density targets instead of solids. Simulations have revealed that the acceleration process depends on the density gradients of the plasma target. Indeed, under certain conditions, the most energetic protons are predicted to be accelerated by a collisionless shock mechanism that significantly increases their energy. We report here the results of a recent experiment dedicated to the study of longitudinal ion acceleration in partially exploded foils using a high intensity (∼5 × 10{sup 18} W/cm{sup 2}) picosecond laser pulse. We show that protons accelerated using targets having moderate front and rear plasma gradients (up to ∼8 μm gradient length) exhibit similar maximum proton energy and number compared to proton beams that are produced, in similar laser conditions, from solid targets, in the well-known target normal sheath acceleration regime. Particle-In-Cell simulations, performed in the same conditions as the experiment and consistent with the measurements, allow laying a path for further improvement of this acceleration scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Design and Simulation of IOTA - a Novel Concept of Integrable Optics Test Accelerator

The use of nonlinear lattices with large betatron tune spreads can increase instability and space charge thresholds due to improved Landau damping. Unfortunately, the majority of nonlinear accelerator lattices turn out to be nonintegrable, producing chaotic motion and a complex network of stable and unstable resonances. Recent advances in finding the integrable nonlinear accelerator lattices have led to a proposal to construct at Fermilab a test accelerator with strong nonlinear focusing which avoids resonances and chaotic particle motion. This presentation will outline the main challenges, theoretical design solutions and construction status of the Integrable Optics Test Accelerator (IOTA) underway at Fermilab.

43 PARTICLE ACCELERATORS↗