Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,801 records · Page 100

Excitation of transverse dipole and quadrupole modes in a pure ion plasma in a linear Paul trap to study collective processes in intense beams

Transverse dipole and quadrupole modes have been excited in a one-component cesium ion plasma trapped in the Paul Trap Simulator Experiment (PTSX) in order to characterize their properties and understand the effect of their excitation on equivalent long-distance beam propagation. The PTSX device is a compact laboratory Paul trap that simulates the transverse dynamics of a long, intense charge bunch propagating through an alternating-gradient transport system by putting the physicist in the beam's frame of reference. A pair of arbitrary function generators was used to apply trapping voltage waveform perturbations with a range of frequencies and, by changing which electrodes were driven with the perturbation, with either a dipole or quadrupole spatial structure. The results presented in this paper explore the dependence of the perturbation voltage's effect on the perturbation duration and amplitude. Perturbations were also applied that simulate the effect of random lattice errors that exist in an accelerator with quadrupole magnets that are misaligned or have variance in their field strength. The experimental results quantify the growth in the equivalent transverse beam emittance that occurs due to the applied noise and demonstrate that the random lattice errors interact with the trapped plasma through the plasma's internal collective modes. Coherent periodic perturbations were applied to simulate the effects of magnet errors in circular machines such as storage rings. The trapped one component plasma is strongly affected when the perturbation frequency is commensurate with a plasma mode frequency. The experimental results, which help to understand the physics of quiescent intense beam propagation over large distances, are compared with analytic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TeV Neutrinos and Hard X-Rays from Relativistic Reconnection in the Corona of NGC 1068

Abstract The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is dn p / dE p ∝ E p − 1 below the break and dn p / dE p ∝ E p − s above the break, with IceCube neutrino observations suggesting s ≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.

Astronomy & Astrophysics↗

Concepts of association between cancer and ionising radiation: accounting for specific biological mechanisms

The probability that an observed cancer was caused by radiation exposure is usually estimated using cancer rates and risk models from radioepidemiological cohorts and is called assigned share (AS). This definition implicitly assumes that an ongoing carcinogenic process is unaffected by the studied radiation exposure. However, there is strong evidence that radiation can also accelerate an existing clonal development towards cancer. In this work, we define different association measures that an observed cancer was newly induced, accelerated, or retarded. The measures were quantified exemplarily by Monte Carlo simulations that track the development of individual cells. Three biologically based two-stage clonal expansion (TSCE) models were applied. In the first model, radiation initiates cancer development, while in the other two, radiation has a promoting effect, i.e. radiation accelerates the clonal expansion of pre-cancerous cells. The parameters of the TSCE models were derived from breast cancer data from the atomic bomb survivors of Hiroshima and Nagasaki. For exposure at age 30, all three models resulted in similar estimates of AS at age 60. For the initiation model, estimates of association were nearly identical to AS. However, for the promotion models, the cancerous clonal development was frequently accelerated towards younger ages, resulting in associations substantially higher than AS. This work shows that the association between a given cancer and exposure in an affected person depends on the underlying biological mechanism and can be substantially larger than the AS derived from classic radioepidemiology.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Hydrodynamic instability growth and mix experiments at the National Ignition Facility

Hydrodynamic instability growth and its effects on implosion performance were studied at the National Ignition Facility [G. H. Miller, E. I. Moses, and C. R. Wuest, Opt. Eng. 443, 2841 (2004)]. Implosion performance and mix have been measured at peak compression using plastic shells filled with tritium gas and containing embedded localized carbon-deuterium diagnostic layers in various locations in the ablator. Neutron yield and ion temperature of the deuterium-tritium fusion reactions were used as a measure of shell-gas mix, while neutron yield of the tritium-tritium fusion reaction was used as a measure of implosion performance. The results have indicated that the low-mode hydrodynamic instabilities due to surface roughness were the primary culprits for yield degradation, with atomic ablator-gas mix playing a secondary role. In addition, spherical shells with pre-imposed 2D modulations were used to measure instability growth in the acceleration phase of the implosions. The capsules were imploded using ignition-relevant laser pulses, and ablation-front modulation growth was measured using x-ray radiography for a shell convergence ratio of ∼2. The measured growth was in good agreement with that predicted, thus validating simulations for the fastest growing modulations with mode numbers up to 90 in the acceleration phase. Future experiments will be focused on measurements at higher convergence, higher-mode number modulations, and growth occurring during the deceleration phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Structural Lattices for a Davis Gun Earth Penetrator Impact Experiment

The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.

36 MATERIALS SCIENCE↗

Laser Wakefield Acceleration: Structural and Dynamic Studies. Final Technical Report ER40954

Particle accelerators enable scientists to study the fundamental structure of the universe, but have become the largest and most expensive of scientific instruments. In this project, we advanced the science and technology of laser-plasma accelerators, which are thousands of times smaller and less expensive than their conventional counterparts. In a laser-plasma accelerator, a powerful laser pulse exerts light pressure on an ionized gas, or plasma, thereby driving an electron density wave, which resembles the wake behind a boat. Electrostatic fields within this plasma wake reach tens of billions of volts per meter, fields far stronger than ordinary non-plasma matter (such as the matter that a conventional accelerator is made of) can withstand. Under the right conditions, stray electrons from the surrounding plasma become trapped within these “wake-fields”, surf them, and acquire energy much faster than is possible in a conventional accelerator. Laser-plasma accelerators thus might herald a new generation of compact, low-cost accelerators for future particle physics, x-ray and medical research. In this project, we made two major advances in the science of laser-plasma accelerators. The first of these was to accelerate electrons beyond 1 gigaelectronvolt (1 GeV) for the first time. In experimental results reported in Nature Communications in 2013, about 1 billion electrons were captured from a tenuous plasma (about 1/100 of atmosphere density) and accelerated to 2 GeV within about one inch, while maintaining less than 5% energy spread, and spreading out less than ½ milliradian (i.e. ½ millimeter per meter of travel). Low energy spread and high beam collimation are important for applications of accelerators as coherent x-ray sources or particle colliders. This advance was made possible by exploiting unique properties of the Texas Petawatt Laser, a powerful laser at the University of Texas at Austin that produces pulses of 150 femtoseconds (1 femtosecond is 10-15 seconds) in duration and 150 Joules in energy (equivalent to the muzzle energy of a small pistol bullet). This duration was well matched to the natural electron density oscillation period of plasma of 1/100 atmospheric density, enabling efficient excitation of a plasma wake, while this energy was sufficient to drive a high-amplitude wake of the right shape to produce an energetic, collimated electron beam. Continuing research is aimed at increasing electron energy even further, increasing the number of electrons captured and accelerated, and developing applications of the compact, multi-GeV accelerator as a coherent, hard x-ray source for materials science, biomedical imaging and homeland security applications. The second major advance under this project was to develop new methods of visualizing the laser-driven plasma wake structures that underlie laser-plasma accelerators. Visualizing these structures is essential to understanding, optimizing and scaling laser-plasma accelerators. Yet prior to work under this project, computer simulations based on estimated initial conditions were the sole source of detailed knowledge of the complex, evolving internal structure of laser-driven plasma wakes. In this project we developed and demonstrated a suite of optical visualization methods based on well-known methods such as holography, streak cameras, and coherence tomography, but adapted to the ultrafast, light-speed, microscopic world of laser-driven plasma wakes. Our methods output images of laser-driven plasma structures in a single laser shot. We first reported snapshots of low-amplitude laser wakes in Nature Physics in 2006. We subsequently reported images of high-amplitude laser-driven plasma “bubbles”, which are important for producing electron beams with low energy spread, in Physical Review Letters in 2010. More recently, we have figured out how to image laser-driven structures that change shape while propagating in a single laser shot. The latter techniques, which use the methods of computerized tomography, were demonstrated on test objects – e.g. laser-driven filaments in air and glass – and reported in Optics Letters in 2013 and Nature Communications in 2014. Their output is a multi-frame movie rather than a snapshot. Continuing research is aimed at applying these tomographic methods directly to evolving laser-driven plasma accelerator structures in our laboratory, then, once perfected, to exporting them to plasma-based accelerator laboratories around the world as standard in-line metrology instruments.

43 PARTICLE ACCELERATORS↗

FEM design and simulation of a short, 10 MV, S-band Linac with Monte Carlo dose simulations

Purpose: Current commercial 10 MV Linac waveguides are 1.5 m. The authors’ current 6 MV linear accelerator–magnetic resonance imager (Linac–MR) system fits in typical radiotherapy vaults. To allow 10 MV treatments with the Linac–MR and still fit within typical vaults, the authors design a 10 MV Linac with an accelerator waveguide of the same length (27.5 cm) as current 6 MV Linacs. Methods: The first design stage is to design a cavity such that a specific experimental measurement for breakdown is applicable to the cavity. This is accomplished through the use of finite element method (FEM) simulations to match published shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength from an electric breakdown study. A full waveguide is then designed and tuned in FEM simulations based on this cavity design. Electron trajectories are computed through the resulting radio frequency fields, and the waveguide geometry is modified by shifting the first coupling cavity in order to optimize the electron beam properties until the energy spread and mean energy closely match values published for an emulated 10 MV Linac. Finally, Monte Carlo dose simulations are used to compare the resulting photon beam depth dose profile and penumbra with that produced by the emulated 10 MV Linac. Results: The shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength are all matched to within 0.1%. A first coupling cavity shift of 1.45 mm produces an energy spectrum width of 0.347 MeV, very close to the published value for the emulated 10 MV of 0.315 MeV, and a mean energy of 10.53 MeV, nearly identical to the published 10.5 MeV for the emulated 10 MV Linac. The depth dose profile produced by their new Linac is within 1% of that produced by the emulated 10 MV spectrum for all depths greater than 1.5 cm. The penumbra produced is 11% narrower, as measured from 80% to 20% of the central axis dose. Conclusions: The authors have successfully designed and simulated an S-band waveguide of length of 27.5 cm capable of producing a 10 MV photon beam. This waveguide operates well within the breakdown threshold determined for the cavity geometry used. The designed Linac produces depth dose profiles similar to those of the emulated 10 MV Linac (waveguide-length of 1.5 m) but yields a narrower penumbra.

07 ISOTOPE AND RADIATION SOURCES↗

Milky Way and M31 rotation curves: Λ CDM versus MOND

We analyze the existing rotation-curve data of the Milky Way and M31 galaxies that extend to very large distances and low accelerations. We find a systematic downward trend in the weak acceleration (large distances) segment of the radial acceleration. A similar downward trend has been noticed in the ΛCDM EAGLE simulation [Phys. Rev. D 96, 124016 (2017).], while the deviation from the generic MOND prediction would need to be ascribed to an external field effect, or possibly a post facto selected acceleration function μ(x).

79 ASTRONOMY AND ASTROPHYSICS↗

Local radiation hydrodynamic simulations of massive star envelopes at the iron opacity peak

We perform three-dimensional radiation hydrodynamic simulations of the structure and dynamics of the radiation-dominated envelopes of massive stars at the location of the iron opacity peak. One-dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, τ{sub 0}, to τ{sub c}=c/c{sub g}, where c{sub g}≈50 km s{sup −1} is the isothermal sound speed in the gas alone. When τ{sub 0}≫τ{sub c}, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation-dominated regime. When τ{sub 0}≪τ{sub c}, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed c{sub g}, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Monte Carlo simulation of the dose response of a novel 2D silicon diode array for use in hybrid MRI–LINAC systems

Purpose: MRI-guided radiation therapy systems (MRIgRT) are being developed to improve online imaging during treatment delivery. At present, the operation of single point dosimeters and an ionization chamber array have been characterized in such systems. This work investigates a novel 2D diode array, named “magic plate,” for both single point calibration and 2D positional performance, the latter being a key element of modern radiotherapy techniques that will be delivered by these systems. Methods: GEANT4 Monte Carlo methods have been employed to study the dose response of a silicon diode array to 6 MV photon beams, in the presence of in-line and perpendicularly aligned uniform magnetic fields. The array consists of 121 silicon diodes (dimensions 1.5 × 1.5 × 0.38 mm{sup 3}) embedded in kapton substrate with 1 cm pitch, spanning a 10 × 10 cm{sup 2} area in total. A geometrically identical, water equivalent volume was simulated concurrently for comparison. The dose response of the silicon diode array was assessed for various photon beam field shapes and sizes, including an IMRT field, at 1 T. The dose response was further investigated at larger magnetic field strengths (1.5 and 3 T) for a 4 × 4 cm{sup 2} photon field size. Results: The magic plate diode array shows excellent correspondence (< ± 1%) to water dose in the in-line orientation, for all beam arrangements and magnetic field strengths investigated. The perpendicular orientation, however, exhibits a dose shift with respect to water at the high-dose-gradient beam edge of jaw-defined fields [maximum (4.3 ± 0.8)% over-response, maximum (1.8 ± 0.8)% under-response on opposing side for 1 T, uncertainty 1σ]. The trend is not evident in areas with in-field dose gradients typical of IMRT dose maps. Conclusions: A novel 121 pixel silicon diode array detector has been characterized by Monte Carlo simulation for its performance inside magnetic fields representative of current prototype and proposed MRI–linear accelerator systems. In the in-line orientation, the silicon dose is directly proportional to the water dose. In the perpendicular orientation, there is a shift in dose response relative to water in the highest dose gradient regions, at the edge of jaw-defined and single-segment MLC fields. The trend was not observed in-field for an IMRT beam. The array is expected to be a valuable tool in MRIgRT dosimetry.

07 ISOTOPE AND RADIATION SOURCES↗

Development and performance of a HemeLB GPU code for human-scale blood flow simulation

In recent years, it has become increasingly common for high performance computers (HPC) to possess some level of heterogeneous architecture - typically in the form of GPU accelerators. In some machines these are isolated within a dedicated partition, whilst in others they are integral to all compute nodes - often with multiple GPUs per node - and provide the majority of a machine's compute performance. In light of this trend, it is becoming essential that codes deployed on HPC are updated to execute on accelerator hardware. Here, in this paper, we introduce a GPU implementation of the 3D blood flow simulation code HemeLB that has been developed using CUDA C++. We demonstrate how taking advantage of NVIDIA GPU hardware can achieve significant performance improvements compared to the equivalent CPU only code on which it has been built whilst retaining the excellent strong scaling characteristics that have been repeatedly demonstrated by the CPU version. With HPC positioned on the brink of the exascale era, we use HemeLB as a motivation to provide a discussion on some of the challenges that many users will face when deploying their own applications on upcoming exascale machines.

59 BASIC BIOLOGICAL SCIENCES↗

Synchronous acceleration with tapered dielectric-lined waveguides

Here, we present a general concept to accelerate non-relativistic charged particles. Our concept employs an adiabatically-tapered dielectric-lined waveguide which supports accelerating phase velocities for synchronous acceleration. We propose an ansatz for the transient field equations, show it satisfies Maxwell's equations under an adiabatic approximation and find excellent agreement with a finite-difference time-domain computer simulation. The fields were implemented into the particle-tracking program {\sc astra} and we present beam dynamics results for an accelerating field with a 1-mm-wavelength and peak electric field of 100~MV/m. The numerical simulations indicate that a $\sim 200$-keV electron beam can be accelerated to an energy of $\sim10$~MeV over $\sim 10$~cm. The novel scheme is also found to form electron beams with parameters of interest to a wide range of applications including, e.g., future advanced accelerators, and ultra-fast electron diffraction.

43 PARTICLE ACCELERATORS↗

Shock–Cloud Interaction in the Solar Corona

Flare-associated coronal shock waves sometimes interact with solar prominences, leading to large-amplitude prominence oscillations (LAPOs). Such prominence activation gives us a unique opportunity to track the time evolution of shock–cloud interaction in cosmic plasmas. Although the dynamics of interstellar shock–cloud interaction has been extensively studied, coronal shock–solar prominence interaction is rarely studied in the context of shock–cloud interaction. Associated with the X5.4 class solar flare that occurred on 2012 March 7, a globally propagated coronal shock wave interacted with a polar prominence, leading to LAPO. In this paper, we studied bulk acceleration and excitation of the internal flow of the shocked prominence using three-dimensional magnetohydrodynamic (MHD) simulations. We studied eight MHD simulation runs, each with different mass density structure of the prominence, and one hydrodynamic simulation run, and we compared the result. In order to compare the observed motion of activated prominence with the corresponding simulation, we also studied prominence activation by injection of a triangular-shaped coronal shock. We found that the prominence is first accelerated mainly by magnetic tension force as well as direct transmission of the shock, and later decelerated mainly by magnetic tension force. The internal flow, on the other hand, is excited during the shock front sweeps through the prominence and damps almost exponentially. We construct a phenomenological model of bulk momentum transfer from the shock to the prominence, which agreed quantitatively with all the simulation results. Based on the phenomenological prominence activation model, we diagnosed physical parameters of the coronal shock wave. The estimated energy of the coronal shock is several percent of the total energy released during the X5.4 flare.

79 ASTRONOMY AND ASTROPHYSICS↗

Water Contamination Impacts on Biodiesel Antioxidants and Storage Stability

Biodiesel (fatty acid methyl esters) can oxidize in storage to form acids and gums that negatively impact engine performance and durability. Antioxidant additives are used to increase biodiesel storage stability, and previous studies that evaluated the effectiveness of antioxidants demonstrated that more polar antioxidants tend to be the most effective (provide the largest improvement in stability per unit of antioxidant added). However, polar antioxidants have significant water solubility, and diesel fuel storage tanks are commonly contaminated with water (forming a layer of water under the fuel). This study investigates whether nonpolar antioxidants, which are less effective in dry environments, might be more effective under wet conditions simulating real-world storage. Biodiesel blends treated with polar and nonpolar antioxidants were subjected to accelerated aging using the ASTM D4625 protocol (storage at 43 °C, open to air, for up to 24 weeks) both with and without added water. Fuels treated with polar/higher-effectiveness compounds and stored in contact with water (simulating water in a storage tank) or high humidity showed accelerated loss of stability compared to dry storage. The same fuel treated with a nonpolar antioxidant and stored in the same conditions did not exhibit accelerated stability loss and thus had higher storage stability over the long term despite treatment with an initially less effective additive. Analysis of the fuels during aging showed that this loss of stability was not due to oxidation but rather extraction of the polar antioxidant into the water layer. Antioxidant additives that are incompatible with wet or humid storage conditions were found to cause faster-than-anticipated loss of stability, which was preventable with use of nonpolar additives.

09 BIOMASS FUELS↗

A review of neutronics and thermal hydraulics–based screening methods applied to accelerated nuclear fuel qualification

This paper reviews the state-of-the-art engineering approach for using thermal hydraulic (TH) and neutronics modeling and simulation (M&S) tools to perform rapid screening studies of novel nuclear fuel concepts within the context of accelerated fuel qualification. Global research efforts have introduced nuclear fuel and material concepts that mark a significant departure from traditional reactor materials. The number of new technologies being considered for development for light water reactors and advanced reactor types has created the need for an accelerated fuel qualification procedure. A key component of this procedure is the rapid identification of the most promising fuel concepts using computational screening studies. Advanced TH and neutronic M&S tools should be leveraged to efficiently determine whether the reactor performance and safety characteristics of a given concept warrant additional studies or whether the concept requires modification or elimination. This paper reviews best practices for performing these TH and neutronics screening studies at various stages during a fuel concept's progression through the qualification procedure. The motivation behind standardizing this approach is to minimize time and resources spent on qualification activities for fuel concepts that could be quickly refined or eliminated from consideration based on their reactor physics and TH characteristics. Adoption of this screening procedure—which focuses primarily on nuclear fuels but may be applicable to other reactor materials—will also help accelerate new material qualification by generating boundary conditions crucial to fuel performance evaluations and highlighting needed areas of separate effects experimentation. This article reviews the motivation behind the introduction of novel nuclear fuel concepts, provides incentive for utilizing TH- and neutronics-based screening studies, describes the screening approach and methodology, and includes discussion on how to interpret screening results to provide recommendations for the continued development of a given concept.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Methods, systems, and apparatuses for calculating global fluence for neutron and photon monte carlo transport using expected value estimators

Global fluence estimators may be calculated on accelerators and processors for neutron and photon Monte Carlo transport. Monte Carlo random walk simulation may be performed on the processors and the calculation of a Volumetric-Ray-Casting (VRC) estimator may be offloaded to the accelerators. The VRC estimator may modify an expected-value estimator to extend a pseudo-particle ray along the direction of the emitted particle from source and collision event through not only the event volume, but also through all volumes that describe the problem geometry. Additionally, many pseudo-particle rays may be sampled per event, rather than just a single pseudo-particle ray per event, in order to provide more complete angular coverage.

Sweezy, Jeremy Ed↗

Spontaneous ion beam formation in the laboratory, space, and simulation

We present experimental evidence for the spontaneous formation of multiple double layers within a single divergent magnetic field structure. Downstream of the divergent magnetic field, multiple accelerated ion populations are observed. The similarity of the accelerated ion populations observed in these laboratory experiments to ion populations observed in the magnetosphere and in numerical simulations suggests that the observation of a complex ion velocity distribution alone is insufficient to distinguish between simple plasma expansion and magnetic reconnection. Further, the effective temperature of the aggregate ion population is significantly larger than the temperatures of the individual ion population components, suggesting that insufficiently resolved measurements could misidentify multiple beam creation as ion heating. Ions accelerated in randomly oriented electric fields that mimic heating would have an ion heating rate dependent on the ion charge and mass that is qualitatively consistent with recent experimental observations of ion heating during magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-T-507: Modeling Cerenkov Emissions From Medical Linear Accelerators: A Monte Carlo Study

Purpose: Cerenkov emissions are a natural byproduct of MV radiotherapy but are typically ignored as inconsequential. However, Cerenkov photons may be useful for activation of drugs such as psoralen. Here, we investigate Cerenkov radiation from common radiotherapy beams using Monte Carlo simulations. Methods: GAMOS, a GEANT4-based framework for Monte Carlo simulations, was used to model 6 and 18MV photon beams from a Varian medical linac. Simulations were run to track Cerenkov production from these beams when irradiating a 50cm radius sphere of water. Electron contamination was neglected. 2 million primary photon histories were run for each energy, and values scored included integral dose and total track length of Cerenkov photons between 100 and 400 nm wavelength. By lowering process energy thresholds, simulations included low energy Bremsstrahlung photons to ensure comprehensive evaluation of UV production in the medium. Results: For the same number of primary photons, UV Cerenkov production for 18MV was greater than 6MV by a factor of 3.72 as determined by total track length. The total integral dose was a factor of 2.31 greater for the 18MV beam. Bremsstrahlung photons were a negligibly small component of photons in the wavelength range of interest, comprising 0.02% of such photons. Conclusion: Cerenkov emissions in water are 1.6x greater for 18MV than 6MV for the same integral dose. Future work will expand the analysis to include optical properties of tissues, and to investigate strategies to maximize Cerenkov emission per unit dose for MV radiotherapy.

60 APPLIED LIFE SCIENCES↗

Refine your search to explore more results.