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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-T-673: Effects of Cardiac Induced Brain Pulsations On Proton Minibeams

Purpose: To quantify the dosimetric impact of internal motion within the brain on spatially modulated proton minibeam radiation therapy (pMRT) for small animal research. Methods: The peak-to-valley dose ratio (PVDR) is an essential dosimetric factor for pMRT. Motion of an animal brain caused by cardiac-induced pulsations (CIP) can impact dose deposition. For synchrotron generated high dose rate X-ray microbeams this effect is evaded due to the quasi-instantaneous delivery. By comparison, pMRT potentially suffers increased spread due to lower dose rates. However, for a given dose rate it is less susceptible to beam spread than microbeams, due to the spatial modulation being an order of magnitude larger. Monte Carlo simulations in TOPAS were used to model the beam spread for a 50.5MeV pMRT beam. Motion effects were simulated for a 50mm thick brass collimator with 0.3mm slit width and 1.0mm center-to-center spacing in a water phantom. The maximum motion in a rat brain due to CIP has been reported to be 0.06mm. Motion was simulated with a peak amplitude in the range 0–0.2mm. Results: The impact of 0.06mm peak motion was minimal and reduced the PVDR by about 1% at a depth of 10mm. For 0.2mm peak motion the PVDR was reduced by 16% at a depth of 10mm. Conclusion: For the pMRT beam the magnitude of cardiac-induced brain motion has minimal impact on the PVDR for the investigated collimator geometry. For more narrow beams the effect is likely to be larger. This indicates that delivery of pMRT to small animal brains should not be affected considerably by beamlines with linac compatible dose rates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Using a particle-in-cell model for accelerator control room applications

Many accelerator control rooms rely on envelope models to simulate beam dynamics because they are fast and accurate at tracking the beam core. Particle-in-Cell (PIC) models, however, can track particles inside and outside the core and, with the improvements of computers, are now fast enough to be used in control rooms. The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is currently developing a tool to use a Particle-in-Cell model for control room applications. This report covers the progress so far and the future goals of using PyORBIT, a Particle-in-Cell simulation model, in the SNS control room.

Cathey, Brandon↗
DOE OSTI2013

Start-2-End Simulation in FACET Linac

Start-2-end simulations are essential to evaluate various beam dynamic aspects. In this note, start-2-end simulation studies are presented for FACET electron linac accelerator, from the damping ring exit point to the final interaction point.

Accelerators↗
DOE OSTI2024

Simulating the CMS High Granularity Calorimeter with ML

Detector simulation is a key component of physics analysis and related activities in CMS. In the upcoming High Luminosity LHC era, simulation will be required to use a smaller fraction of computing in order to satisfy resource constraints. At the same time, CMS will be upgraded with the new High Granularity Calorimeter (HGCal), which requires significantly more resources to simulate than the existing CMS calorimeters. This computing challenge motivates the use of generative machine learning models as surrogates to replace full physics-based simulation. We study the application of state-of-the-art diffusion models to simulate particle showers in the CMS HGCal. We will discuss methods to overcome the challenges posed by the high-dimensional, irregular geometry of the HGCal. The quality of the showers produced by the diffusion model will be assessed by comparison to the full GEANT4-based simulation. The increase in simulation throughput will be quantified and methods to accelerate the diffusion model inference will also be discussed.

Amram, Oz↗
DOE OSTI2016

Solving the {sup 56}Ni puzzle of magnetar-powered broad-lined type Ic supernovae

Broad-lined Type Ic supernovae (SNe Ic-BL) are of great importance because their association with long-duration gamma-ray bursts (LGRBs) holds the key to deciphering the central engine of LGRBs, which refrains from being unveiled despite decades of investigation. Among the two popularly hypothesized types of central engine, i.e., black holes and strongly magnetized neutron stars (magnetars), there is mounting evidence that the central engine of GRB-associated SNe (GRB-SNe) is rapidly rotating magnetars. Theoretical analysis also suggests that magnetars could be the central engine of SNe Ic-BL. What puzzled the researchers is the fact that light-curve modeling indicates that as much as 0.2--0.5 M{sub ⊙} of {sup 56}Ni was synthesized during the explosion of the SNe Ic-BL, which is unfortunately in direct conflict with current state-of-the-art understanding of magnetar-powered {sup 56}Ni synthesis. Here we propose a dynamic model of magnetar-powered SNe to take into account the acceleration of the ejecta by the magnetar, as well as the thermalization of the injected energy. Assuming that the SN kinetic energy comes exclusively from the magnetar acceleration, we find that although a major fraction of the rotational energy of the magnetar is to accelerate the SN ejecta, a tiny fraction of this energy deposited as thermal energy of the ejecta is enough to reduce the needed {sup 56}Ni to 0.06 M {sub ⊙} for both SN 1997ef and SN 2007ru. We therefore suggest that magnetars could power SNe Ic-BL in aspects both of energetics and of {sup 56}Ni synthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Turbulent Transport in a Three-dimensional Solar Wind

Turbulence in the solar wind can play essential roles in the heating of coronal and solar wind plasma and the acceleration of the solar wind and energetic particles. Turbulence sources are not well understood and thought to be partly enhanced by interaction with the large-scale inhomogeneity of the solar wind and the interplanetary magnetic field and/or transported from the solar corona. To investigate the interaction with background inhomogeneity and the turbulence sources, we have developed a new 3D MHD model that includes the transport and dissipation of turbulence using the theoretical model of Zank et al. We solve for the temporal and spatial evolution of three moments or variables, the energy in the forward and backward fluctuating modes and the residual energy and their three corresponding correlation lengths. The transport model is coupled to our 3D model of the inhomogeneous solar wind. We present results of the coupled solar wind-turbulence model assuming a simple tilted dipole magnetic configuration that mimics solar minimum conditions, together with several comparative intermediate cases. By considering eight possible solar wind and turbulence source configurations, we show that the large-scale solar wind and IMF inhomogeneity and the strength of the turbulence sources significantly affect the distribution of turbulence in the heliosphere within 6 au. We compare the predicted turbulence distribution results from a complete solar minimum model with in situ measurements made by the Helios and Ulysses spacecraft, finding that the synthetic profiles of the turbulence intensities show reasonable agreement with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Balmer Filaments in Tycho’s Supernova Remnant: An Interplay between Cosmic-ray and Broad-neutral Precursors

We present Hα spectroscopic observations and detailed modeling of the Balmer filaments in the supernova remnant (SNR) Tycho (SN 1572). We used GH α FaS (Galaxy Hα Fabry–Pérot Spectrometer) on the William Herschel Telescope with a 3.′4 × 3.′4 field of view, 0.″2 pixel scale, and σ{sub instr}=8.1 km s{sup −1} resolution at 1″ seeing for ∼10 hr, resulting in 82 spatial−spectral bins that resolve the narrow Hα line in the entire SN 1572 northeastern rim. For the first time, we can therefore mitigate artificial line broadening from unresolved differential motion and probe Hα emission parameters in varying shock and ambient medium conditions. Broad Hα line remains unresolved within spectral coverage of 392 km s{sup −1}. We employed Bayesian inference to obtain reliable parameter confidence intervals and to quantify the evidence for models with multiple line components. The median Hα narrow-line (NL) FWHM of all bins and models is W{sub NL}=(54.8±1.8) km s{sup −1} at the 95% confidence level, varying within [35, 72] km s{sup −1} between bins and clearly broadened compared to the intrinsic (thermal) ≈20 km s{sup −1}. Possible line splits are accounted for, significant in ≈18% of the filament, and presumably due to remaining projection effects. We also find widespread evidence for intermediate-line emission of a broad-neutral precursor, with a median W{sub IL}=(180±14) km s{sup −1} (95% confidence). Finally, we present a measurement of the remnant’s systemic velocity, V{sub LSR}=−34 km s{sup −1}, and map differential line-of-sight motions. Our results confirm the existence and interplay of shock precursors in Tycho’s remnant. In particular, we show that suprathermal NL emission is near-universal in SN 1572, and that, in the absence of an alternative explanation, collisionless SNR shocks constitute a viable acceleration source for Galactic TeV cosmic-ray protons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The extent of power-law energy spectra in collisionless relativistic magnetic reconnection in pair plasma

Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron–positron plasmas, for a wide range of upstream magnetizations σ and system sizes L. The particle spectra are well-represented by a power law ${\gamma }^{-\alpha }$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to σ and L, respectively. As a result, for large L and σ, the power-law index α approaches about 1.2.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

ImpactX v0.1

ImpactX is the next generation of the IMPACT-Z code. It is a s-based simulation code for modeling intense beams in particle accelerators using symplectic tracking methods and includes collective effects. It is multi-node parallel and supports modern compute hardware such as GPUs, modern algorithms such as mesh-refinement and realistic geometries (embedded boundaries).

Huebl, Axel↗
DOE OSTI2016

Spin Tracking of Polarized Protons in the Main Injector at Fermilab

The Main Injector (MI) at Fermilab currently produces high-intensity beams of protons at energies of 120 GeV for a variety of physics experiments. Acceleration of polarized protons in the MI would provide opportunities for a rich spin physics program at Fermilab. To achieve polarized proton beams in the Fermilab accelerator complex, shown in Fig.1.1, detailed spin tracking simulations with realistic parameters based on the existing facility are required. This report presents studies at the MI using a single 4-twist Si-berian snake to determine the depolarizing spin resonances for the relevant synchrotrons. Results will be presented first for a perfect MI lattice, followed by a lattice that includes the real MI imperfections, such as the measured magnet field errors and quadrupole misalignments. The tolerances of each of these factors in maintaining polariza-tion in the Main Injector will be discussed.

Xiao, M.↗
DOE OSTI2017

Roofline Analysis in the Intel® Advisor to Deliver Optimized Performance for applications on Intel® Xeon Phi™ Processor

In this session we show, in two case studies, how the roofline feature of Intel Advisor has been utilized to optimize the performance of kernels of the XGC1 and PICSAR codes in preparation for Intel Knights Landing architecture. The impact of the implemented optimizations and the benefits of using the automatic roofline feature of Intel Advisor to study performance of large applications will be presented. This demonstrates an effective optimization strategy that has enabled these science applications to achieve up to 4.6 times speed-up and prepare for future exascale architectures. # Goal/Relevance of Session The roofline model [1,2] is a powerful tool for analyzing the performance of applications with respect to the theoretical peak achievable on a given computer architecture. It allows one to graphically represent the performance of an application in terms of operational intensity, i.e. the ratio of flops performed and bytes moved from memory in order to guide optimization efforts. Given the scale and complexity of modern science applications, it can often be a tedious task for the user to perform the analysis on the level of functions or loops to identify where performance gains can be made. With new Intel tools, it is now possible to automate this task, as well as base the estimates of peak performance on measurements rather than vendor specifications. The goal of this session is to demonstrate how the roofline feature of Intel Advisor can be used to balance memory vs. computation related optimization efforts and effectively identify performance bottlenecks. A series of typical optimization techniques: cache blocking, structure refactoring, data alignment, and vectorization illustrated by the kernel cases will be addressed. # Description of the codes ## XGC1 The XGC1 code [3] is a magnetic fusion Particle-In-Cell code that uses an unstructured mesh for its Poisson solver that allows it to accurately resolve the edge plasma of a magnetic fusion device. After recent optimizations to its collision kernel [4], most of the computing time is spent in the electron push (pushe) kernel, where these optimization efforts have been focused. The kernel code scaled well with MPI+OpenMP but had almost no automatic compiler vectorization, in part due to indirect memory addresses and in part due to low trip counts of low-level loops that would be candidates for vectorization. Particle blocking and sorting have been implemented to increase trip counts of low-level loops and improve memory locality, and OpenMP directives have been added to vectorize compute-intensive loops that were identified by Advisor. The optimizations have improved the performance of the pushe kernel 2x on Haswell processors and 1.7x on KNL. The KNL node-for-node performance has been brought to within 30% of a NERSC Cori phase I Haswell node and we expect to bridge this gap by reducing the memory footprint of compute intensive routines to improve cache reuse. ## PICSAR is a Fortran/Python high-performance Particle-In-Cell library targeting at MIC architectures first designed to be coupled with the PIC code WARP for the simulation of laser-matter interaction and particle accelerators. PICSAR also contains a FORTRAN stand-alone kernel for performance studies and benchmarks. A MPI domain decomposition is used between NUMA domains and a tile decomposition (cache-blocking) handled by OpenMP has been added for shared-memory parallelism and better cache management. The so-called current deposition and field gathering steps that compose the PIC time loop constitute major hotspots that have been rewritten to enable more efficient vectorization. Particle communications between tiles and MPI domain has been merged and parallelized. All considered, these improvements provide speedups of 3.1 for order 1 and 4.6 for order 3 interpolation shape factors on KNL configured in SNC4 quadrant flat mode. Performance is similar between a node of cori phase 1 and KNL at order 1 and better on KNL by a factor 1.6 at order 3 with the considered test case (homogeneous thermal plasma).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

DEVELOPMENT OF ENSEMBLE-BASED PROBABILISTIC FLOOD MAPS FOR PROBABLE MAXIMUM FLOOD

The magnitude and frequency of hydro-meteorological extremes are expected to increase in a changing environment that threaten the security of US energy-water assets. This includes probable maximum precipitation (PMP) and probable maximum flood (PMF) which are used as hydraulic design standards for critical infrastructures such as major hydropower reservoirs and nuclear power plants. To assess the potential flood hazards due to PMP/PMF, an integrated high resolution process-based hydro-meteorologic modeling framework was utilized to develop ensemble-based probabilistic flood maps based on best-available historic observations and future climate projections. A graphical process unit accelerated 2-dimensional hydrodynamic model was used to simulate the surface inundation areas corresponding to a total of PMF hydrographs. These ensemble-based PMF maps were compared with flood maps obtained from the conventional deterministic PMP/PMF approach, revealing added information of conditional probability of flooding. Further, a relative sensitivity test was conducted to explore the effects of various factors in the framework such as meteorological forcings, antecedent hydrologic conditions, reservoir storage, and flood model input resolution and parameters. The proposed framework better illustrates the uncertainties associated with model inputs, parameterization and hydro-meteorological factors allowing more informed decision-making for future emergency preparation.

13 HYDRO ENERGY↗
DOE OSTI2013

Accelerator shield design of KIPT neutron source facility

Argonne National Laboratory (ANL) of the United States and Kharkov Institute of Physics and Technology (KIPT) of Ukraine have been collaborating on the design development of a neutron source facility at KIPT utilizing an electron-accelerator-driven subcritical assembly. Electron beam power is 100 kW, using 100 MeV electrons. The facility is designed to perform basic and applied nuclear research, produce medical isotopes, and train young nuclear specialists. The biological shield of the accelerator building is designed to reduce the biological dose to less than 0.5-mrem/hr during operation. The main source of the biological dose is the photons and the neutrons generated by interactions of leaked electrons from the electron gun and accelerator sections with the surrounding concrete and accelerator materials. The Monte Carlo code MCNPX serves as the calculation tool for the shield design, due to its capability to transport electrons, photons, and neutrons coupled problems. The direct photon dose can be tallied by MCNPX calculation, starting with the leaked electrons. However, it is difficult to accurately tally the neutron dose directly from the leaked electrons. The neutron yield per electron from the interactions with the surrounding components is less than 0.01 neutron per electron. This causes difficulties for Monte Carlo analyses and consumes tremendous computation time for tallying with acceptable statistics the neutron dose outside the shield boundary. To avoid these difficulties, the SOURCE and TALLYX user subroutines of MCNPX were developed for the study. The generated neutrons are banked, together with all related parameters, for a subsequent MCNPX calculation to obtain the neutron and secondary photon doses. The weight windows variance reduction technique is utilized for both neutron and photon dose calculations. Two shielding materials, i.e., heavy concrete and ordinary concrete, were considered for the shield design. The main goal is to maintain the total dose outside the shield boundary at less than 0.5-mrem/hr. The shield configuration and parameters of the accelerator building have been determined and are presented in this paper. (authors)

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Rich-Text based PSA Model Transformation and Verification in RiskA

Nuclear Energy is one of the most possible alternative to fossil energy in large-scale. Nuclear safety is the lifeline of nuclear energy development, and has been an important part of national security, especially after the accident of Fukushima. On the basis of computer codes, reliability engineering and Probabilistic Risk Assessment (PSA) is effective methods to improve the safety of nuclear power plants. After many years of development, the methodologies, applications and relevant soft programs have been quite sufficient in the area of nuclear safety. However, it is still difficult for the current PSA technologies and software to meet the needs of dealing with the problems of quantitatively public risk assessment in the digital society. Especially, there are so many software that adopted different modeling and analyzing methods, which brings difficulties to the comparison, evaluation or verification for different software and models. By using its own academic advantages, the advanced nuclear energy research team, i.e., FDS Team has invested a few hundreds person-years of independent research and development of a reliability and probabilistic safety assessment software system RiskA. The system consists of three main functions: device reliability optimization management, probabilistic safety assessment and instantaneous risk monitoring, and the comprehensive public risk assessment in digital society. RiskA establishes an integrated risk evaluation system for the device reliability, nuclear power plant (NPP) risk monitoring and public risk assessment. Currently, RiskA has successfully been applied in more than 600 organizations and institutes from a variety of fields in 39 countries, involving nuclear, aerospace, defense and so on. For example, RiskA has been used in the safety analysis of International Thermonuclear Experimental Reactor (ITER), FDS series fusion reactor, and the Strategic Priority Research Program of Chinese Academy of Sciences 'the Future Advanced Nuclear Fission Energy - Accelerator Driven Sub-critical (ADS) transmutation system'. Also, RiskA has been used in public risk communication and risk assessment. Among them, Third Qinshan nuclear power plant Risk Monitor (TQRM), which is the China's first domestic nuclear safety analysis software used in the actual online operation, has been highly apprised by renowned experts in the field. Meanwhile, it has been identified as the first domestic nuclear power plant risk monitoring system with completely independent intellectual property rights, and has reached the international advanced level. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

The Myth of Fungible FTE: A Quantitative Assessment of Matrixed Resource Allocation

Matrix organizations allow scientific facilities to share specialized personnel across projects, operations, maintenance, and strategic initiatives. Nominal staffing allocations, however, may not capture the schedule consequences of fragmented individual commitments, limited access to specialist groups, and intermittent availability of key decision makers. We developed a stochas- tic, daily-time-step simulation of a hypothetical medium-sized accelerator-facility project com- prising sequential phases and parallel tasks. Each task requires role-specific work measured in FTE-days. Ordinary personnel may be unavailable because they contribute concurrently to other institutional activities, while designated key roles have independently specified daily un- availability probabilities. An organization-wide priority factor scales the number of people from each functional group who can effectively contribute to the project. It is interpreted as a composite proxy for project access and workforce fragmentation across competing commit- ments. We examined project completion time as a function of this factor and Project Lead unavailability using 100 Monte Carlo runs per condition. Increasing priority factor from 0.1 to 1.0 reduced median completion time from 1708.5 days (interquartile range 1681.5–1735.25) to 390 days (interquartile range 379–399). At priority factor = 0.1, increasing Project Lead unavailability from 0.5 to 0.9 increased median completion time from 1713.5 days (interquartile range 1691–1733.25) to 4,417 days (interquartile range 4271.75–4550.5). The model quantifies the commonly expected sensitivity of project schedules to fragmented resource commitments and limited coordination availability. Within this model, the results also indicate a possible threshold regime in which small increases in workforce availability yield only modest sched- ule improvements until sufficient capacity becomes accessible, after which project performance improves sharply. With further validation and calibration, this quantitative framework could support resource-allocation decisions during initial project planning and subsequent schedule rebaselining.

Bai, Mei [SLAC National Accelerator Laboratory (SL↗