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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Acceleration of relativistic electrons by magnetohydrodynamic turbulence: Implications for non-thermal emission from black hole accretion disks

We use analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low-luminosity systems such as Sgr A* and M87.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Analytical solutions for radiation-driven winds in massive stars. I. The fast regime

Accurate mass-loss rate estimates are crucial keys in the study of wind properties of massive stars and for testing different evolutionary scenarios. From a theoretical point of view, this implies solving a complex set of differential equations in which the radiation field and the hydrodynamics are strongly coupled. The use of an analytical expression to represent the radiation force and the solution of the equation of motion has many advantages over numerical integrations. Therefore, in this work, we present an analytical expression as a solution of the equation of motion for radiation-driven winds in terms of the force multiplier parameters. This analytical expression is obtained by employing the line acceleration expression given by Villata and the methodology proposed by Müller and Vink. On the other hand, we find useful relationships to determine the parameters for the line acceleration given by Müller and Vink in terms of the force multiplier parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE MAGNETIZATION DEGREE OF THE OUTFLOW POWERING THE HIGHLY POLARIZED REVERSE-SHOCK EMISSION OF GRB 120308A

GRB 120308A, a long duration γ-ray burst (GRB) detected by Swift, was distinguished by a highly polarized early optical afterglow emission that strongly suggests an ordered magnetic field component in the emitting region. In this work, we model the optical and X-ray emission in the reverse and forward shock scenario and show that the strength of the magnetic field in the reverse-shock region is ∼10 times stronger than that in the forward shock region. Consequently, the outflow powering the highly polarized reverse-shock optical emission was mildly magnetized at a degree of σ ∼ a few percent. Considering the plausible magnetic energy dissipation in both the acceleration and prompt emission phases of the GRB outflow, the afterglow data of GRB 120308A provides us with compelling evidence that, at least for some GRBs, a nonignorable fraction of the energy was released in the form of Poynting flux, confirming the finding first made in the reverse-forward shock emission modeling of the optical afterglow of GRB 990123 by Fan et al. in 2002 and Zhang et al. in 2003.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE EFFECT OF TURBULENCE INTERMITTENCE ON THE EMISSION OF SOLAR ENERGETIC PARTICLES BY CORONAL AND INTERPLANETARY SHOCKS

Major solar energetic particle events are associated with shock waves in solar corona and solar wind. Fast scattering of charged particles by plasma turbulence near the shock wave increases the efficiency of the particle acceleration in the shock, but prevents particles from escaping ahead of the shock. However, the turbulence energy levels in neighboring magnetic tubes of solar wind may differ from each other by more than one order of magnitude. We present the first theoretical study of accelerated particle emission from an oblique shock wave propagating through an intermittent turbulence background that consists of both highly turbulent magnetic tubes, where particles are accelerated, and quiet tubes, via which the accelerated particles can escape to the non-shocked solar wind. The modeling results imply that the presence of the fast transport channels penetrating the shock and cross-field transport of accelerated particles to those channels may play a key role in high-energy particle emission from distant shocks and can explain the prompt onset of major solar energetic particle events observed near the Earth's orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Domain-specific compilers for dynamic simulations of quantum materials on quantum computers

Abstract Simulation of the dynamics of quantum materials is emerging as a promising scientific application for noisy intermediate-scale quantum (NISQ) computers. Due to their high gate-error rates and short decoherence times, however, NISQ computers can only produce high-fidelity results for those quantum circuits smaller than some given circuit size. Dynamic simulations, therefore, pose a challenge as current algorithms produce circuits that grow in size with each subsequent time-step of the simulation. This underscores the crucial role of quantum circuit compilers to produce executable quantum circuits of minimal size, thereby maximizing the range of physical phenomena that can be studied within the NISQ fidelity budget. Here, we present two domain-specific (DS) quantum circuit compilers for the Rigetti and IBM quantum computers, specifically designed to compile circuits simulating dynamics under a special class of time-dependent Hamiltonians. The compilers outperform state-of-the-art general-purpose compilers in terms of circuit size reduction by around 25%–30% as well as wall-clock compilation time by around 40% (dependent on system size and simulation time-step). Drawing on heuristic techniques commonly used in artificial intelligence, both compilers scale well with simulation time-step and system size. Code for both compilers is open-source and packaged into a full-stack quantum simulation software with tutorials included for ease of use for future researchers wishing to perform dynamic simulations of quantum materials on quantum computers. As our DS compilers provide significant improvements in both compilation time and simulation fidelity, they provide a building block for accelerating progress toward physical quantum supremacy.

Physics↗
DOE OSTI2024

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A general Bayesian algorithm for the autonomous alignment of beamlines

Autonomous methods to align beamlines can decrease the amount of time spent on diagnostics, and also uncover better global optima leading to better beam quality. The alignment of these beamlines is a high-dimensional expensive-to-sample optimization problem involving the simultaneous treatment of many optical elements with correlated and nonlinear dynamics. Bayesian optimization is a strategy of efficient global optimization that has proved successful in similar regimes in a wide variety of beamline alignment applications, though it has typically been implemented for particular beamlines and optimization tasks. In this paper, we present a basic formulation of Bayesian inference and Gaussian process models as they relate to multi-objective Bayesian optimization, as well as the practical challenges presented by beamline alignment. We show that the same general implementation of Bayesian optimization with special consideration for beamline alignment can quickly learn the dynamics of particular beamlines in an online fashion through hyperparameter fitting with no prior information. We present the implementation of a concise software framework for beamline alignment and test it on four different optimization problems for experiments on X-ray beamlines at the National Synchrotron Light Source II and the Advanced Light Source, and an electron beam at the Accelerator Test Facility, along with benchmarking on a simulated digital twin. We discuss new applications of the framework, and the potential for a unified approach to beamline alignment at synchrotron facilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

Addendum to the AAPM's TG-51 protocol for clinical reference dosimetry of high-energy photon beams

An addendum to the AAPM's TG-51 protocol for the determination of absorbed dose to water in megavoltage photon beams is presented. This addendum continues the procedure laid out in TG-51 but new k{sub Q} data for photon beams, based on Monte Carlo simulations, are presented and recommendations are given to improve the accuracy and consistency of the protocol's implementation. The components of the uncertainty budget in determining absorbed dose to water at the reference point are introduced and the magnitude of each component discussed. Finally, the consistency of experimental determination of N{sub D,w} coefficients is discussed. It is expected that the implementation of this addendum will be straightforward, assuming that the user is already familiar with TG-51. The changes introduced by this report are generally minor, although new recommendations could result in procedural changes for individual users. It is expected that the effort on the medical physicist's part to implement this addendum will not be significant and could be done as part of the annual linac calibration.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-T-128: Dosimetric Evaluation of MLC Modeling in Pinnacle V9.2 for Varian TrueBeam STx

Purpose: To evaluate the effects of MLC modeling after commissioning the Varian TrueBeam LINAC in Pinnacle version 9.2. Methods: Stepand-shoot IMRT QAs were investigated when we observed our measured absolute dose results using ion chamber (Capintec PR-05P) were uncharacteristically low; about 4–5% compared to doses calculated by Pinnacle{sup 3} (Phillips, Madison, WI). This problem was predominant for large and highly modulated head and neck (HN) treatments. Intuitively we knew this had to be related to shortcomings in the MLC modeling in Pinnacle. Using film QA we were able to iteratively adjust the MLC parameters. We confirmed results by re-testing five failed IMRT QA patients; and ion chamber measurements were verified in Quasar anthropomorphic phantom. Results: After commissioning the LINAC in Pinnacle version 9.2, the MLC transmission for 6X, 10X and 15X were 2.0%, 1.7% and 2.0%, respectively, and additional Interleaf leakage for all three energies was 0.5%. These parameters were obtained from profiles scanned with an Edge detector (Sun Nuclear, Melbourne, FL) during machine commissioning. A Verification testing with radiographic EDR2 film (Kodak, Rochester, NY) measurement was performed by creating a closed MLC leaf pattern and analyzing using RIT software (RIT, Colorado Springs, CO). This reduced MLC transmission for 6X, 10X and 15X to 0.7%, 0.9% and 0.9%, respectively; while increasing additional Interleaf leakage for all three energies to 1.0%. Conclusion: Radiographic film measurements were used to correct MLC transmission values for step and shoot IMRT fields used in Pinnacle version 9.2. After adjusting the MLC parameters to correlate with the film QA, there was still very good agreement between the Pinnacle model and commissioning data. Using the same QA methodology, we were also able to improve the beam models for the Varian C-series linacs, Novalis-Tx, and TrueBeam M-120 linacs.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-T-250: Determining VMAT Machine Limitations of An Elekta Linear Accelerator with Agility MLC for Accurate Modeling in RayStation and Robust Delivery

Purpose: To implement VMAT in RayStation with the Elekta Synergy linac with the new Agility MLC, and to utilize the same vendor softwares to determine the optimum Elekta VMAT machine parameters in RayStation for accurate modeling and robust delivery. Methods: iCOMCat is utilized to create various beam patterns with user defined dose rate, gantry, MLC and jaw speed for each control point. The accuracy and stability of the output and beam profile are qualified for each isolated functional component of VMAT delivery using ion chamber and Profiler2 with isocentric mounting fixture. Service graphing on linac console is used to verify the mechanical motion accuracy. The determined optimum Elekta VMAT machine parameters were configured in RayStation v4.5.1. To evaluate the system overall performance, TG-119 test cases and nine retrospective VMAT patients were planned on RayStation, and validated using both ArcCHECK (with plug and ion chamber) and MapCHECK2. Results: Machine output and profile varies <0.3% when only variable is dose rate (35MU/min-600MU/min). <0.9% output and <0.3% profile variation are observed with additional gantry motion (0.53deg/s–5.8deg/s both directions). The output and profile variation are still <1% with additional slow leaf motion (<1.5cm/s both direction). However, the profile becomes less symmetric, and >1.5% output and 7% profile deviation is seen with >2.5cm/s leaf motion. All clinical cases achieved comparable plan quality as treated IMRT plans. The gamma passing rate is 99.5±0.5% on ArcCheck (<3% iso center dose deviation) and 99.1±0.8% on MapCheck2 using 3%/3mm gamma (10% lower threshold). Mechanical motion accuracy in all VMAT deliveries is <1°/1mm. Conclusion: Accurate RayStation modeling and robust VMAT delivery is achievable on Elekta Agility for <2.5cm/s leaf motion and full range of dose rate and gantry speed determined by the same vendor softwares. Our TG-119 and patient results have provided us with the confidence to use VMAT clinically.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-340: Dosimetry of a Small Field Electron Beam for Innovative Radiotherapy of Small Surface Or Internal Tumors

Purpose: An innovative small high intensity electron beams with energies from 6 to 12 MeV is being developed at Argonne National Laboratory to deliver an absorbed dose via a catheter to small malignant and nonmalignant lesions. This study reports on the initial dosimetric characteristics of this electron beam. These include output calibration, percent depth dose, beam profiles and leakage through the catheter. Methods: To simulate the narrow electron beam, the Argonne Wakefield Accelerator is used to produce high energy electron beams. The electron beam from the accelerator is monitored by measuring the current through a transmission coil while the beam shape is observed with a fluorescent screen. The dosimetry properties of the electron beam transmitting through bone and tissue-like materials are measured with nanodot optically stimulated luminescent dosimeters and EDR radiographic film. The 6 MV photon beam from a Varian True beam linac is used to calibrate both the OSLDs and the film. Results: The beam characteristics of the 12 MeV beam were measured. The properties of the small diameter, 5 mm, beam differs from that of broad clinical electron beams from radiotherapy linacs. Due to the lack of scatter from the narrow beam, the maximum dose is at the surface and the depth of the 50% depth dose is 35 mm compared to 51 mm for a clinical 12 MeV. The widths of the 90% isodose measured at the surface and depths of 2, 6, 12, and 16 mm varied from 6.6 to 8.8 mm while the widths of the FWHM isodose varied from 7.8 to 25.5 mm. Conclusion: Initial beam measurements show favorable dosimetric properties for its use in treating either small surface or internal lesions, particularly to deliver radiation at the time of surgery to maximize the dose to the lesion and spare normal tissue.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-360: End-To-End Dosimetric Testing of a Versa HD Linear Accelerator with the Agility Head Modeled in Pinnacle3

Purpose: The Versa HD incorporates a variety of upgrades, primarily including the Agility head. The distinct dosimetric properties of the head from its predecessors combined with flattening-filter-free (FFF) beams require a new investigation of modeling in planning systems and verification of modeling accuracy. Methods: A model was created in Pinnacle{sup 3} v9.8 with commissioned beam data. Leaf transmission was modeled as <0.5% with maximum leaf speed of 3 cm/s. Photon spectra were tuned for FFF beams, for which profiles were modeled with arbitrary profiles rather than with cones. For verification, a variety of plans with varied parameters were devised, and point dose measurements were compared to calculated values. A phantom of several plastic water and Styrofoam slabs was scanned and imported into Pinnacle{sup 3}. Beams of different field sizes, SSD, wedges, and gantry angles were created. All available photon energies (6 MV, 10 MV, 18 MV, 6 FFF, 10 FFF) as well four clinical electron energies (6, 9, 12, and 15 MeV) were investigated. The plans were verified at a calculation point (8 cm deep for photons, variable for electrons) by measurement with a PTW Semiflex ionization chamber. In addition, IMRT testing was performed with three standard plans (step and shoot IMRT, small and large field VMAT plans). The plans were delivered on the Delta4 IMRT QA phantom (ScandiDos, Uppsala, Sweden). Results: Homogeneous point dose measurement agreed within 2% for all photon and electron beams. Open field photon measurements along the central axis at 100 cm SSD passed within 1%. Gamma passing rates were >99.5% for all plans with a 3%/3mm tolerance criteria. The IMRT QA results for the first 23 patients yielded gamma passing rates of 97.4±2.3%. Conclusion: The end-to-end testing ensured confidence in the ability of Pinnacle{sup 3} to model photon and electron beams with the Agility head.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-298: The Impact of Modeling the Treatment Couch On Patient Specific VMAT QA

Purpose: The aim of this work is to quantify the impact of modeling the treatment couch on the passing rate of ion chamber measurements during VMAT quality assurance. Methods: For the initial characterization, attenuation and surface dose measurements were performed following the guidelines of TG-176 for the Civco Universal couch top using an Elekta VersaHD accelerator at an energy of 6 MV. A simulation CT was performed to aid in the creation of contours for representing the shape and size of the couch top in the treatment planning system (TPS). A uniform value of density for the couch wall was determined by comparing the ratios of ion chamber measurements made in a 30×30×11 cm3 water phantom with the TPS dose values of a plan with the same geometry. At our institution, patient specific quality assurance is performed using a Sun Nuclear ArcCheck with a multi-plug for chamber measurements, a 0.125cc PTW TN31010 chamber, and a Sun Nuclear 1010 electrometer. Ten VMAT plans were transferred into the phantom geometry created in the TPS with two settings: with and without the couch. The chamber measurements were compared to both treatment plans. Results: A maximum attenuation of 3.6% was observed when the gantry angle was set to 120 and 240 degrees, passing obliquely through the couch. A uniform density of 0.6 g/cm3 for the couch wall was determined in the TPS by comparison with measured data. The VMAT ion chamber measurement/plan ratios systematically improved by 1.79% ±0.53% for all patients when the couch was included in the calculation. Conclusion: The attenuation and surface dose changes produced by the Civco couch can generate observable dose difference in VMAT plans. Including a couch model in the phantom plan used for patient specific VMAT QA can improve the ionization chamber agreement by up to ∼2%.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-656: Monte Carlo Study On Air Activation Around a Medical Electron Linac

Purpose: In high energy photon therapy, several radiation protection issues result from photonuclear reactions. The activation of air - directly by photonuclear reactions as well as indirectly by capture of photoneutrons generated inside the linac head - is a major point of concern for the medical staff. The purpose of this study was to estimate the annual effective dose to medical workers due to activated air around a medical high energy electron linac by means of Monte Carlo simulations. Methods: The treatment head of a Varian Clinac in 18 MV-X mode as well as the surrounding concrete bunker were modeled and the radiation transport was simulated using the Monte Carlo code FLUKA, starting from the primary electron striking the bremsstrahlung target. The activation yields in air from photo-disintegration of O-16 and N-14 nuclei as well as from neutron capture on Ar-40 nuclei were obtained from the simulations. The activation build-up, radioactive decay and air ventilation were studied using a mathematical model. The annual effective dose to workers was estimated by using published isotope specific conversion factors. Results: The oxygen and nitrogen activation yields were in contrast to the argon activation yield found to be field size dependent. The impact of the treatment room ventilation on the different air activation products was investigated and quantified. An estimate with very conservative assumptions gave an annual effective dose to workers of < 1 mSv/a. Conclusion: From the results of this study it can be concluded that the contribution of air activation to the radiation exposure to medical workers should be negligible in modern photon therapy, especially when it is compared to the dose due to prompt neutrons and the activation of heavy solid materials such as the jaws and the collimators inside the linac head.

60 APPLIED LIFE SCIENCES↗