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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-T-523: On the Radiobiological Impact of Lateral Scatter in Proton Beams

Introduction: In proton therapy, justified concern has been voiced with respect to an increased efficiency in cell kill at the distal end of the Bragg peak. This coupled with range uncertainty is a counter indication to use the Bragg peak to define the border of a treated volume with a critical organ. An alternative is to use the lateral edge of the proton beam, obtaining more robust plans. We investigate the spectral and biological effects of the lateral scatter . Methods: A general purpose Monte Carlo simulation engine (MCNPX 2.7c) installed on a Scientific Linux cluster, calculated the dose deposition spectrum of protons, knock on electrons and generated neutrons for a proton beam with maximal kinetic energy of 200MeV. Around the beam at different positions in the beam direction the spectrum is calculated in concentric rings of thickness 1cm. The deposited dose is converted to a double strand break map using an analytical expression.based on micro dosimetric calculations using a phenomenological Monte Carlo code (MCDS). A strict version of RBE is defined as the ratio of generation of double strand breaks in the different modalities. To generate the reference a Varian linac was modelled in MCNPX and the generated electron dose deposition spectrum was used . Results: On a pristine point source 200MeV beam the RBE before the Bragg peak was of the order of 1.1, increasing to 1.7 right behind the Bragg peak. When using a physically more realistic beam of 10cm diameter the effect was smaller. Both the lateral dose and RBE increased with increasing beam depth, generating a dose deposition with mixed biological effect. Conclusions: The dose deposition in proton beams need to be carefully examined because the biological effect will be different depending on the treatment geometry. Deeply penetrating proton beams generate more biologically effective lateral scatter.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

WE-G-17A-03: MRIgRT: Quantification of Organ Motion

Purpose: To develop an MRI-based methodology and tools required for the quantification of organ motion on a dedicated MRI-guided radiotherapy system. A three-room facility, consisting of a TrueBeam 6X linac vault, a 1.5T MR suite and a brachytherapy interventional room, is currently under commissioning at our institution. The MR scanner can move and image in either room for diagnostic and treatment guidance purposes. Methods: A multi-imaging modality (MR, kV) phantom, featuring programmable 3D simple and complex motion trajectories, was used for the validation of several image sorting algorithms. The testing was performed on MRI (e.g. TrueFISP, TurboFLASH), 4D CT and 4D CBCT. The image sorting techniques were based on a) direct image pixel manipulation into columns or rows, b) single and aggregated pixel data tracking and c) using computer vision techniques for global pixel analysis. Subsequently, the motion phantom and sorting algorithms were utilized for commissioning of MR fast imaging techniques for 2D-cine and 4D data rendering. MR imaging protocols were optimized (e.g. readout gradient strength vs. SNR) to minimize the presence of susceptibility-induced distortions, which were reported through phantom experiments and numerical simulations. The system-related distortions were also quantified (dedicated field phantom) and treated as systematic shifts where relevant. Results: Image sorting algorithms were validated for specific MR-based applications such as quantification of organ motion, local data sampling, and 4D MRI for pre-RT delivery with accuracy better than the raw image pixel size (e.g. 1 mm). MR fast imaging sequences were commissioning and imaging strategies were developed to mitigate spatial artifacts with minimal penalty on the image spatial and temporal sampling. Workflows (e.g. liver) were optimized to include the new motion quantification tools for RT planning and daily patient setup verification. Conclusion: Comprehensive methods were developed and validated for the quantification of organ motion with applications in MRI-guided RT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-C-304-03: Experimental Investigation On the Accuracy of Plastic Scintillation Dosimeters in Small Fields

Purpose: To investigate the accuracy of the Exradin W1 (SI) and of an “in-house” plastic scintillation dosimeter (CHUQ PSD) in small radiation fields. Methods: Output factor (OF) measurements with the W1 and CHUQ PSD were performed for field sizes of 0.5 x 0.5, 1 x 1 and 2 x 2 cm{sup 2}. Both detectors were placed parallel to the central axis (CAX) in water. The spectrum discrimination calibration method was performed in each set-up to account for the Cerenkov (CRV) signal created in the fiber. The OFs were compared to the expected field factors in water derived using i) Monte Carlo (MC) simulations of an accurate accelerator model and ii) microLion (PTW) and D1V diode (SI) OFs. MC-derived correction factors were applied to both the microLion and D1V OFs. For the CHUQ PSD the calibration was repeated in water (// CAX), solid water (perpendicular to CAX) and under a shielded configuration. The signal was collected using a spectrometer (wavelength range = 185–1100 nm). Spectral analysis was performed to evaluate potential changes of the spectral distributions under the various calibration set-up configurations. Results: The W1 OFs presented an over-response for the 0.5 x 0.5 cm{sup 2} in the range of 3 – 4.1% relative to the expected field factor. The CHUQ PSD presented an under-response in the range of 1.5 – 2.7%, without accounting for volume averaging. The CRV spectra under the various calibration procedures appeared similar to each other and only minor changes were observed to the respective OFs. Conclusion: The W1 and CHUQ PSD can be used in small fields down to a 1 x 1 cm{sup 2} field size. Discrepancies were encountered between the two detectors for the smallest field size of 0.5 x 0.5 cm{sup 2} with the CHUQ PSD exhibiting a closer agreement to the expected field factor. Funding sources: 1) Alexander S. Onassis Public Benefit Foundation in Greece and 2) CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council, Grant number: 432290, RGPIN-2014-06475.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Machine learning for data-driven discovery in solid Earth geoscience

Understanding the behavior of Earth through the diverse fields of the solid Earth geosciences is an increasingly important task. It is made challenging by the complex, interacting, and multiscale processes needed to understand Earth’s behavior and by the inaccessibility of nearly all of Earth’s subsurface to direct observation. Substantial increases in data availability and in the increasingly realistic character of computer simulations hold promise for accelerating progress, but developing a deeper understanding based on these capabilities is itself challenging. Machine learning will play a key role in this effort. We review the state of the field and make recommendations for how progress might be broadened and accelerated.

Bergen, Karianne J.↗
DOE OSTI2022

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗
DOE OSTI2018

Numerical Simulations of Space Charge Compensation with an Electron Lens

The future high energy physics program at Fermilab requires that the proton complex operate with beam bunch intensities four times larger than is currently handled. At these intensities space charge nonlinear defocussing effects cause unacceptable particle losses especially in the low energy rapid-cycling-synchrotron (RCS) Booster. Focusing electron lens elements may offer a solution by providing partial space charge compensation but there is a need for detailed simulations as this technique has not been demonstrated. We report on high fidelity numerical 6D space charge simulations in a model accelerator lattice with a record high space charge tune shift approaching unity.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Simulation of Transparent Spin Experiment in RHIC

The transparent spin mode has been proposed as a new technique for preservation and control of the spin polarization of an ion beam in a synchrotron. The ion rings of the proposed Jefferson Lab Electron-Ion Collider adopted this technique in their figure-8 design. The transparent spin mode can also be setup in a racetrack with two identical Siberian snakes. There is a proposal to test the predicted features of the spin transparent mode in RHIC, which already has all of the necessary hardware capabilities. We have earlier analytically estimated the setup parameters and developed a preliminary experimental plan. In this paper we describe simulations of the proposed experiment using a Zgoubi model of RHIC. We present simulation of the polarization preservation during acceleration in the transparent spin mode. We also demonstrate polarization control by simulating an adiabatic polarization adjustment. A procedure to measure the strength of the spin transparent resonance is illustrated.

Morozov, Vasiliy↗
DOE OSTI2016

ER@CEBAF: Modeling code developments

A proposal for a multiple-pass, high-energy, energy-recovery experiment using CEBAF is under preparation in the frame of a JLab-BNL collaboration. In view of beam dynamics investigations regarding this project, in addition to the existing model in use in Elegant a version of CEBAF is developed in the stepwise ray-tracing code Zgoubi, Beyond the ER experiment, it is also planned to use the latter for the study of polarization transport in the presence of synchrotron radiation, down to Hall D line where a 12 GeV polarized beam can be delivered. This Note briefly reports on the preliminary steps, and preliminary outcomes, based on an Elegant to Zgoubi translation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Simulation of an FFA return channel ERL in a RHIC-based electron-ion collider

The linac-ring electron-ion collider designed at BNL is based on a 20 GeV ERL comprised of a 1.3 GeV linac and two 3.8km FFA lattice return loops placed along the RHIC heavy ion collider. This document introduces the simulation of the FFA ERL in Zgoubi, including a start-to-end 6-D simulation of the acceleration-deceleration cycle of a polarized electron bunch.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

[Office of Basic Energy Sciences (BES)] Roundtable on Producing and Managing Large Scientific Data with Artificial Intelligence and Machine Learning

To identify specific Priority Research Opportunities (PROs) for artificial intelligence and machine learning (AI/ML) at the US Department of Energy’s scientific user facilities, the Office of Basic Energy Sciences (BES) convened a roundtable of facility experts encompassing the fields of physics, chemistry, materials synthesis science, computational science, detector and accelerator technology, theory, modeling, simulation, and atomic-scale characterization techniques. The roundtable met on October 22–23, 2019, to identify coordinated, long-term AI/ML research efforts that will drive major advances in neutron, photon, and nanoscale sciences. This report describes the four PROs identified at the roundtable: PRO 1 on how AI/ML can extract high-value information from the large datasets; PRO 2 on how AI/ML can use such information in real time to maximize the facilities’ scientific output; PRO 3 on using AI/ML virtual laboratories (i.e., computational models of experimental facilities) to aid the facilities and user community in design and control of machine parameters and the design and execution of experiments, including training AI/ML models for PROs 1 and 2; and PRO 4 on how a shared scientific data infrastructure can provide tools to assemble and analyze the totality of data coming from user facilities. A section of the end of the report provides a summary on computer science and mathematics that highlights where enhanced AI/ML capabilities could be particularly impactful for BES user facilities.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2014

Modeling Long-Term Changes in Climate, Ice Sheets and Sea Level: Using the Paleo Record to Understand Possibilities for the Future

During the project, we completed the development and testing of the infrastructure for coupled CESM-CISM paleo-simulations. Our development and testing efforts have focused on (1) CISM2, a parallel, scalable code with a suite of higher-order ice-flow solvers, (2) conservative, two-way coupling of CISM and CLM, (3) an efficient spinup method for coupled CESM-CISM simulations, and (4) inclusion of processes important for long, millennial-scale simulations (dynamic topography, orbital acceleration, asynchronous coupling to dynamic vegetation) and their evaluation (isotopes).

58 GEOSCIENCES↗
DOE OSTI2021

Predictive Battery Lifetime Modeling at NREL [Slides]

Battery lifetime models are used to extrapolate data from accelerated aging tests to simulate degradation in real-world applications such as electric vehicles and battery energy storage systems. Methods developed at NREL utilize both expert domain-knowledge and machine-learning to identify models, using statistical methods such as cross-validation and bootstrap resampling to interrogate model performance and quantify uncertainty. These models can be utilized in systems level simulations to predict battery performance or technoeconomic models to estimate the lifetime cost of battery systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗
DOE OSTI2016

( Fermi Fermi)/LAT study of gamma-ray emission in the direction of the Monoceros Loop supernova remnant

We present an analysis of the gamma-ray measurements by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope in the region of the supernova remnant (SNR) Monoceros Loop (G205.5+0.5). The brightest gamma-ray peak is spatially correlated with the Rosette Nebula, which is a molecular cloud complex adjacent to the southeast edge of the SNR. After subtraction of this emission by spatial modeling, the gamma-ray emission from the SNR emerges, which is extended and fit by a Gaussian spatial template. The gamma-ray spectra are significantly better reproduced by a curved shape than a simple power law. The luminosities between 0.2 and 300 GeV are ∼4×10{sup 34} erg s{sup −1} for the SNR and ∼3×10{sup 34} erg s{sup −1} for the Rosette Nebula, respectively. We argue that the gamma-rays likely originate from the interactions of particles accelerated in the SNR. The decay of neutral pions produced in nucleon–nucleon interactions of accelerated hadrons with interstellar gas provides a reasonable explanation for the gamma-ray emission of both the Rosette Nebula and the Monoceros SNR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Structure of Energetic Particle Mediated Shocks Revisited

The structure of collisionless shock waves is often modified by the presence of energetic particles that are not equilibrated with the thermal plasma (such as pickup ions [PUIs] and solar energetic particles [SEPs]). This is relevant to the inner and outer heliosphere and the Very Local Interstellar Medium (VLISM), where observations of shock waves (e.g., in the inner heliosphere) show that both the magnetic field and thermal gas pressure are less than the energetic particle component pressures. Voyager 2 observations revealed that the heliospheric termination shock (HTS) is very broad and mediated by energetic particles. PUIs and SEPs contribute both a collisionless heat flux and a higher-order viscosity. We show that the incorporation of both effects can completely determine the structure of collisionless shocks mediated by energetic ions. Since the reduced form of the PUI-mediated plasma model is structurally identical to the classical cosmic ray two-fluid model, we note that the presence of viscosity, at least formally, eliminates the need for a gas sub-shock in the classical two-fluid model, including in that regime where three are possible. By considering parameters upstream of the HTS, we show that the thermal gas remains relatively cold and the shock is mediated by PUIs. We determine the structure of the weak interstellar shock observed by Voyager 1 . We consider the inclusion of the thermal heat flux and viscosity to address the most general form of an energetic particle-thermal plasma two-fluid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

PARTICLE ACCELERATION AND HEATING BY TURBULENT RECONNECTION

Turbulent flows in the solar wind, large-scale current sheets, multiple current sheets, and shock waves lead to the formation of environments in which a dense network of current sheets is established and sustains “turbulent reconnection.” We constructed a 2D grid on which a number of randomly chosen grid points are acting as scatterers (i.e., magnetic clouds or current sheets). Our goal is to examine how test particles respond inside this large-scale collection of scatterers. We study the energy gain of individual particles, the evolution of their energy distribution, and their escape time distribution. We have developed a new method to estimate the transport coefficients from the dynamics of the interaction of the particles with the scatterers. Replacing the “magnetic clouds” with current sheets, we have proven that the energization processes can be more efficient depending on the strength of the effective electric fields inside the current sheets and their statistical properties. Using the estimated transport coefficients and solving the Fokker–Planck (FP) equation, we can recover the energy distribution of the particles only for the stochastic Fermi process. We have shown that the evolution of the particles inside a turbulent reconnecting volume is not a solution of the FP equation, since the interaction of the particles with the current sheets is “anomalous,” in contrast to the case of the second-order Fermi process.

79 ASTRONOMY AND ASTROPHYSICS↗