Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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,765 records · Page 98

Radiation Transport Tools for Space Applications: A Review

This slide presentation contains a brief discussion of nuclear transport codes widely used in the space radiation community for shielding and scientific analyses. Seven radiation transport codes that are addressed. The two general methods (i.e., Monte Carlo Method, and the Deterministic Method) are briefly reviewed.

Cosmic Ray Effects on Micro-Electronics (CREME96)↗

Use of SCALE Continuous-Energy Monte Carlo Tools for Eigenvalue Sensitivity Coefficient Calculations

The TSUNAMI code within the SCALE code system makes use of eigenvalue sensitivity coefficients for an extensive number of criticality safety applications, such as quantifying the data-induced uncertainty in the eigenvalue of critical systems, assessing the neutronic similarity between different critical systems, and guiding nuclear data adjustment studies. The need to model geometrically complex systems with improved fidelity and the desire to extend TSUNAMI analysis to advanced applications has motivated the development of a methodology for calculating sensitivity coefficients in continuous-energy (CE) Monte Carlo applications. The CLUTCH and Iterated Fission Probability (IFP) eigenvalue sensitivity methods were recently implemented in the CE KENO framework to generate the capability for TSUNAMI-3D to perform eigenvalue sensitivity calculations in continuous-energy applications. This work explores the improvements in accuracy that can be gained in eigenvalue and eigenvalue sensitivity calculations through the use of the SCALE CE KENO and CE TSUNAMI continuous-energy Monte Carlo tools as compared to multigroup tools. The CE KENO and CE TSUNAMI tools were used to analyze two difficult models of critical benchmarks, and produced eigenvalue and eigenvalue sensitivity coefficient results that showed a marked improvement in accuracy. The CLUTCH sensitivity method in particular excelled in terms of efficiency and computational memory requirements.

97 MATHEMATICS AND COMPUTING↗

Monte-Carlo simulation of noise in hard X-ray Transmission Crystal Spectrometers: Identification of contributors to the background noise and shielding optimization

Transmission crystal spectrometers (TCS) are used on many laser facilities to record hard X-ray spectra. During experiments, signal recorded on imaging plates is often degraded by a background noise. Monte-Carlo simulations made with the code GEANT4 show that this background noise is mainly generated by diffusion of MeV electrons and very hard X-rays. An experiment, carried out at LULI2000, confirmed that the use of magnets in front of the diagnostic, that bent the electron trajectories, reduces significantly this background. The new spectrometer SPECTIX (Spectromètre PETAL à Cristal en TransmIssion X), built for the LMJ/PETAL facility, will include this optimized shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Material Assessment for ITER's Collective Thomson Scattering first mirror

The International Thermonuclear Energy Reactor (ITER) Collective Thomson Scattering (CTS) system is a diagnostic instrument that measures plasma density and velocity through Thomson scattering of microwave radiation. Some of the key components of the CTS are quasi-optical mirrors that are used to produce astigmatic beam patterns, which have impact on the strength and spatial resolution of the diagnostic signal. The mirrors are exposed to neutron radiation, which may alter the quality of the signal received. In this work, three different materials (molybdenum (Mo), stainless steel 316 (SS-316) and tungsten (W)) are considered for the first mirror of the CTS. The objective is to access which of the material studied are best suited for this mirror, considering different neutron radiation loads simulated scenarios defined by ITER, based on the resultant stresses and temperature distributions. For it, the neutron irradiation, and subsequent heat-load on the mirrors are simulated using the Monte Carlo N-Particle (MCNP) code. Based on the MCNP heat-load results, transient thermal-structural Finite Element Analysis (FEA) of the mirror over a 400 s discharge, with and without cooling on the rear side, are conducted using in commercial FEA software ANSYS. Results show that of the tested materials Mo and W are the most suitable material for this application. Even though, this study does not yet consider the variation of the material properties with temperature, it presents a quick initial satisfactory assessment that may be considered in subsequent and more complex analysis. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Monte Carlo simulation of photon buildup factors for shielding materials in radiotherapy x-ray facilities

Purpose: This paper presents the results of a series of calculations to determine buildup factors for ordinary concrete, baryte concrete, lead, steel, and iron in broad beam geometry for photons energies from 0.125 to 25.125 MeV at 0.250 MeV intervals.Methods: Monte Carlo N-particle radiation transport computer code has been used to determine the buildup factors for the studied shielding materials.Results: The computation of the primary broad beams using buildup factors data was done for nine published megavoltage photon beam spectra ranging from 4 to 25 MV in nominal energies, representing linacs made by the three major manufacturers. The first tenth value layer and the equilibrium tenth value layer are calculated from the broad beam transmission for these nine primary megavoltage photon beam spectra.Conclusions: The results, compared with published data, show the ability of these buildup factor data to predict shielding transmission curves for the primary radiation beam. Therefore, the buildup factor data can be combined with primary, scatter, and leakage x-ray spectra to perform computation of broad beam transmission for barriers in radiotherapy shielding x-ray facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗

Activation estimates of the RHIC beam dump

In anticipation of the expiration of the RHIC ASE in December 2025 and the need to demonstrate compliance with nuclear facility rule (10CFR830), this study quantifies the radionuclide inventory and respective activity at the beam dumps by simulating a period of five years of operation using conservative assumptions. Simulations were performed with the Monte Carlo particle transport and interaction code FLUKA and its dosimetry results were compared with measured dose rate data, showing good agreement. The beam dump was deemed to constitute the main hotspot of residual radiation and also representative of a worst case scenario in terms of induced activation. The activation levels calculated indicate that the RHIC tunnel will be well below the thresholds of a DOE’s category-3 nuclear facility classification.

43 PARTICLE ACCELERATORS↗

MO-AB-BRA-02: A Novel Scatter Imaging Modality for Real-Time Image Guidance During Lung SBRT

Purpose: A novel scatter imaging modality is developed and its feasibility for image-guided radiation therapy (IGRT) during stereotactic body radiation therapy (SBRT) for lung cancer patients is assessed using analytic and Monte Carlo models as well as experimental testing. Methods: During treatment, incident radiation interacts and scatters from within the patient. The presented methodology forms an image of patient anatomy from the scattered radiation for real-time localization of the treatment target. A radiographic flat panel-based pinhole camera provides spatial information regarding the origin of detected scattered radiation. An analytical model is developed, which provides a mathematical formalism for describing the scatter imaging system. Experimental scatter images are acquired by irradiating an object using a Varian TrueBeam accelerator. The differentiation between tissue types is investigated by imaging simple objects of known compositions (water, lung, and cortical bone equivalent). A lung tumor phantom, simulating materials and geometry encountered during lung SBRT treatments, is fabricated and imaged to investigate image quality for various quantities of delivered radiation. Monte Carlo N-Particle (MCNP) code is used for validation and testing by simulating scatter image formation using the experimental pinhole camera setup. Results: Analytical calculations, MCNP simulations, and experimental results when imaging the water, lung, and cortical bone equivalent objects show close agreement, thus validating the proposed models and demonstrating that scatter imaging differentiates these materials well. Lung tumor phantom images have sufficient contrast-to-noise ratio (CNR) to clearly distinguish tumor from surrounding lung tissue. CNR=4.1 and CNR=29.1 for 10MU and 5000MU images (equivalent to 0.5 and 250 second images), respectively. Conclusion: Lung SBRT provides favorable treatment outcomes, but depends on accurate target localization. A comprehensive approach, employing multiple simulation techniques and experiments, is taken to demonstrate the feasibility of a novel scatter imaging modality for the necessary real-time image guidance.

60 APPLIED LIFE SCIENCES↗

GPU-BASED MONTE CARLO DUST RADIATIVE TRANSFER SCHEME APPLIED TO ACTIVE GALACTIC NUCLEI

A three-dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing-time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons. Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one- and two-dimensional benchmarks are reproduced at high precision. The parallelization capability of various MC algorithms is analyzed and included in our treatment. We utilize the Lucy algorithm for the optical thin case where the Poisson noise is high, the iteration-free Bjorkman and Wood method to reduce the calculation time, and the Fleck and Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGNs) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 {mu}m silicate features in absorption or emission is discussed. The SED of the radio-loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far-infrared emission.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-412: Evaluation of Tungsten-Based Functional Paper for Attenuation Device in Intraoperative Radiotherapy for Breast Cancer

Purpose: Intraoperative radiotherapy (IORT) with an electron beam is one of the accelerated partial breast irradiation methods that have recently been used in early-stage breast cancer. A protective acrylic resin-copper disk is inserted between the breast tissue and the pectoralis muscle to minimize the dose received by the posterior structures. However, a problem with this protective disk is that the surgical incision must be larger than the field size because the disk is manufactured from stiff and unyielding materials. The purpose of this study was to assess the applicability of a new tungsten-based functional paper (TFP) as an alternative to the existing protective disk in IORT. Methods: The newly introduced TFP (Toppan Printing Co., Ltd., Tokyo, JP) is anticipated to become a useful device that is lead-free, light, flexible, and easily processed. The radiation shielding performance of TFP was verified by experimental measurements and Monte Carlo (MC) simulations using PHITS code. The doses transmitted through the protective disk or TFP were measured on a Mobetron mobile accelerator. The same geometries were then reproduced, and the dose distributions were simulated by the MC method. Results: The percentages of transmitted dose relative to the absence of the existing protective disk were lower than 2% in both the measurements and MC simulations. In the experimental measurements, the percentages of transmitted dose for a 9 MeV electron beam were 48.1, 2.3, and 0.6% with TFP thicknesses of 1.9, 3.7, and 7.4 mm, respectively. The percentages for a 12 MeV were 76.0, 49.3, 20.0, and 5.5% with TFP thicknesses of 1.9, 3.7, 7.4, and 14.8 mm, respectively. The results of the MC simulation showed a slight dose increase at the incident surface of the TFP caused by backscattered radiation. Conclusion: The results indicate that a small-incision procedure may be possible by the use of TFP.

60 APPLIED LIFE SCIENCES↗

SU-F-T-577: Comparison of Small Field Dosimetry Measurements in Fields Shaped with Conical Applicators On Two Different Accelerating Systems

Purpose: To investigate small field dosimetry measurements and associated uncertainties when conical applicators are used to shape treatment fields from two different accelerating systems. Methods: Output factor measurements are made in water in beams from the CyberKnife radiosurgery system, which uses conical applicators to shape fields from a (flattening filter-free) 6 MV beam, and in a 6 MV beam from the Elekta Precise linear accelerator (with flattening filter) with BrainLab external conical applicators fitted to shape the field. The measurements use various detectors: (i) an Exradin A16 ion chamber, (ii) two Exradin W1 plastic scintillation detectors, (iii) a Sun Nuclear Edge diode, and (iv) two PTW microDiamond synthetic diamond detectors. Profiles are used for accurate detector positioning and to specify field size (FWHM). Output factor measurements are corrected with detector specific correction factors taken from the literature where available and/or from Monte Carlo simulations using the EGSnrc code system. Results: Differences in measurements of up to 1.7% are observed with a given detector type in the same beam (i.e., intra-detector variability). Corrected results from different detectors in the same beam (inter-detector differences) show deviations up to 3 %. Combining data for all detectors and comparing results from the two accelerators results in a 5.9% maximum difference for the smallest field sizes (FWHM=5.2–5.6 mm), well outside the combined uncertainties (∼1% for the smallest beams) and/or differences among detectors. This suggests that the FWHM of a measured profile is not a good specifier to compare results from different small fields with the same nominal energy. Conclusion: Large differences in results for both intra-detector variability and inter-detector differences suggest potentially high uncertainties in detector-specific correction factors. Differences between the results measured in circular fields from different accelerating systems provide insight into sources of variability in small field dosimetric measurements reported in the literature.

07 ISOTOPE AND RADIATION SOURCES↗

Optimal Temperature Grid For Accurate Doppler Kernel Reconstruction

For computational efficiency, both Monte Carlo and deterministic neutron transport codes often rely on pre-tabulated cross section sets and simple interpolation schemes to represent the Doppler broadening of the cross sections at a needed temperature. Accurate accounting of Doppler effects is of paramount importance for reactor analysis since it impacts the power distribution and the reactor response during transients. Temperature interpolation methods for Doppler broadening have traditionally been physically motivated engineering curve-fits, selected upon best-performance criteria on important nuclides. Kernel reconstruction methods have recently been introduced to perform fast and efficient Doppler broadening interpolation as the result of an optimization problem,, . Out of all the possible kernels analyzed, the temperature-optimized free Doppler kernel reconstruction method was shown to significantly outperform traditional Doppler broadening interpolation methods. Detailed derivations will be part of an upcoming journal paper. In this article, key formulae and results for performing free Doppler kernel reconstruction are summarized, and a study of the optimal temperature grid (T{sub j}) selection is presented. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A New Light Weight Structural Material for Nuclear Structures

Radiation shielding materials are commonly used in nuclear facilities to attenuate the background ionization radiations to a minimum level for creating a safer workplace, meeting regulatory requirements and maintaining high quality performance. The conventional radiation shielding materials have a number of drawbacks: heavy concrete contains a high amount of elements that are not desirable for an effective shielding such as oxygen, silicon, and calcium; a well known limitation of lead is its low machinability and toxicity, which is causing a major environmental concern. Therefore, an effective and environmentally friendly shielding material with increased attenuation and low mass density is desirable. Close-cell composite metal foams (CMFs) and open-cell Al foam with fillers are light-weight candidate materials that we have studied in this project. Close-cell CMFs possess several suitable properties that are unattainable by conventional radiation shielding materials such as low density and high strength for structural applications, high surface area to volume ratio for excellent thermal isolation with an extraordinary energy absorption capability. Open-cell foam is made up of a network of interconnected solid struts, which allows gas or fluid media to pass through it. This unique structure provided a further motive to investigate its application as radiation shields by infiltrating original empty pores with high hydrogen or boron compounds, which are well known for their excellent neutron shielding capability. The resulting open-cell foam with fillers will not only exhibit light weight and high specific surface area, but also possess excellent radiation shielding capability and good processability. In this study, all the foams were investigated for their radiation shielding efficiency in terms of X-ray, gamma ray and neutron. X-ray transmission measurements were carried out on a high-resolution microcomputed tomography (microCT) system. Gamma-emitting sources: 3.0mCi 60 Co, 1.8mCi 137 Cs, 13.5mCi 241 Am, and 5.0mCi 133 Ba were used for gamma-ray attenuation analysis. The evaluations of neutron transmission measurements were conducted at the Neutron Powder Diffractometer beam facility at North Carolina State University. The experimental results were verified theoretically through XCOM and Monte Carlo Z-particle Transport Code (MCNP). A mechanical investigation was performed by means of quasi-static compressive testing. Thermal characterizations were carried out through effective thermal conductivity and thermal expansion analyses in terms of high temperature guarded-comparative-longitudinal heat flow technique and thermomechanical analyzer (TMA), respectively. The experimental results were compared with analytical results obtained from, respectively, Brailsford and Major’s model and modified Turner’s model for verification. Flame test was performed in accordance with United States Nuclear Regulatory Commission (USNRC) standard. CMF sample and a 304L stainless steel control sample were subjected to a fully engulfing fire with an average flame temperature of 800°C for a period of 30 minutes. Finite Element Analysis was conducted to secure the credibility of the experimental results. This research indicates the potential of utilizing the light-weight close-cell CMFs and open-cell Al foam with fillers as shielding material replacing current heavy structures with additional advantage of high-energy absorption and excellent thermal characteristics.

36 MATERIALS SCIENCE↗

DEAD, UNDEAD, AND ZOMBIE ZONES IN PROTOSTELLAR DISKS AS A FUNCTION OF STELLAR MASS

We investigate the viability of the magnetorotational instability (MRI) in X-ray ionized viscous accretion disks around both solar-type stars and very low mass stars. In particular, we determine the disk regions where the MRI can be shut off either by Ohmic resistivity (the so-called dead and undead zones) or by ambipolar diffusion (a region we term the zombie zone). We consider two stellar masses: M {sub *} = 0.7 M {sub Sun} and 0.1 M {sub Sun }. In each case, we assume that: the disk surface density profile is that of a scaled Minimum Mass Solar Nebula, with M {sub disk}/M {sub *} = 0.01 as suggested by current data; disk ionization is driven primarily by stellar X-rays, complemented by cosmic rays and radionuclides; and the stellar X-ray luminosity scales with bolometric luminosity as L{sub X} /L {sub *} Almost-Equal-To 10{sup -3.5}, as observed. Ionization rates are calculated with the MOCCASIN Monte Carlo X-ray transport code, and ionization balance determined using a simplified chemical network, including well-mixed 0.1 {mu}m grains at various levels of depletion. We find that (1) ambipolar diffusion is the primary factor controlling MRI activity in disks around both solar-type and very low mass classical T Tauri stars. Assuming that the MRI yields the maximum possible field strength at each radius, we further find that: (2) the MRI-active layer constitutes only {approx}5%-10% of the total disk mass; (3) the accretion rate ( M-dot ) varies radially in both magnitude and sign (inward or outward), implying time-variable accretion as well as the creation of disk gaps and overdensities, with consequences for planet formation and migration; (4) achieving the empirical accretion rates in solar-type and very low mass stars requires a depletion of well-mixed small grains (via grain growth and/or settling) by a factor of 10-1000 relative to the standard dust-to-gas mass ratio of 10{sup -2}; and (5) the current non-detection of polarized emission from field-aligned grains in the outer disk regions is consistent with active MRI at those radii.

79 ASTRONOMY AND ASTROPHYSICS↗

On the suitability of ultrathin detectors for absorbed dose assessment in the presence of high-density heterogeneities

Purpose: The aim of this study was to evaluate the suitability of several detectors for the determination of absorbed dose in bone. Methods: Three types of ultrathin LiF-based thermoluminescent dosimeters (TLDs)—two LiF:Mg,Cu,P-based (MCP-Ns and TLD-2000F) and a{sup 7}Li-enriched LiF:Mg,Ti-based (MTS-7s)—as well as EBT2 Gafchromic films were used to measure percentage depth-dose distributions (PDDs) in a water-equivalent phantom with a bone-equivalent heterogeneity for 6 and 18 MV and a set of field sizes ranging from 5×5 cm{sup 2} to 20×20 cm{sup 2}. MCP-Ns, TLD-2000F, MTS-7s, and EBT2 have active layers of 50, 20, 50, and 30 μm, respectively. Monte Carlo (MC) dose calculations (PENELOPE code) were used as the reference and helped to understand the experimental results and to evaluate the potential perturbation of the fluence in bone caused by the presence of the detectors. The energy dependence and linearity of the TLDs’ response was evaluated. Results: TLDs exhibited flat energy responses (within 2.5%) and linearity with dose (within 1.1%) within the range of interest for the selected beams. The results revealed that all considered detectors perturb the electron fluence with respect to the energy inside the bone-equivalent material. MCP-Ns and MTS-7s underestimated the absorbed dose in bone by 4%–5%. EBT2 exhibited comparable accuracy to MTS-7s and MCP-Ns. TLD-2000F was able to determine the dose within 2% accuracy. No dependence on the beam energy or field size was observed. The MC calculations showed that a50 μm thick detector can provide reliable dose estimations in bone regardless of whether it is made of LiF, water or EBT’s active layer material. Conclusions: TLD-2000F was found to be suitable for providing reliable absorbed dose measurements in the presence of bone for high-energy x-ray beams.

60 APPLIED LIFE SCIENCES↗

Polyethylene density uncertainty in the simulation of neutron detectors

Simulations with the Monte Carlo N-Particle® (MCNP®) code are frequently used to study neutron coincidence counters for international safeguards. Reducing sources of uncertainty enables more accurate simulations, which reduces the need for expensive or impossible measurements. Accurately specifying the density of polyethylene is crucial for accurate simulations. The density of polyethylene ranges from 0.88 to 0.97 g/cm 3 and high-density polyethylene (HDPE) ranges from about 0.944 to 0.965 g/cm 3 . The density of a specific detector should be measured exactly for the best simulations. Here, the density of the blue/white Active Well Coincidence Counter (AWCC) bulk polyethylene was measured to be 0.9594 +- 0.0013 g/cm 3 at 1-sigma confidence. The density was calculated from repeated measurements of the various polyethylene dimensions for volume, and weighing the polyethylene for the mass. The AWCC neutron detection efficiency was simulated for each polyethylene density. The relative range of efficiencies (3σ) was 3.4% for unspecified polyethylene, 0.9% for HDPE, and 0.3% for the measured density. The impact of modeling small features such as screws was also studied and was found to be a negligible 0.2%. Measuring polyethylene density can reduce its role in uncertainty from one of the largest to one of the smallest contributors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Monte Carlo simulation of TrueBeam flattening-filter-free beams using Varian phase-space files: Comparison with experimental data

Purpose: Phase-space files for Monte Carlo simulation of the Varian TrueBeam beams have been made available by Varian. The aim of this study is to evaluate the accuracy of the distributed phase-space files for flattening filter free (FFF) beams, against experimental measurements from ten TrueBeam Linacs. Methods: The phase-space files have been used as input in PRIMO, a recently released Monte Carlo program based on thePENELOPE code. Simulations of 6 and 10 MV FFF were computed in a virtual water phantom for field sizes 3 × 3, 6 × 6, and 10 × 10 cm{sup 2} using 1 × 1 × 1 mm{sup 3} voxels and for 20 × 20 and 40 × 40 cm{sup 2} with 2 × 2 × 2 mm{sup 3} voxels. The particles contained in the initial phase-space files were transported downstream to a plane just above the phantom surface, where a subsequent phase-space file was tallied. Particles were transported downstream this second phase-space file to the water phantom. Experimental data consisted of depth doses and profiles at five different depths acquired at SSD = 100 cm (seven datasets) and SSD = 90 cm (three datasets). Simulations and experimental data were compared in terms of dose difference. Gamma analysis was also performed using 1%, 1 mm and 2%, 2 mm criteria of dose-difference and distance-to-agreement, respectively. Additionally, the parameters characterizing the dose profiles of unflattened beams were evaluated for both measurements and simulations. Results: Analysis of depth dose curves showed that dose differences increased with increasing field size and depth; this effect might be partly motivated due to an underestimation of the primary beam energy used to compute the phase-space files. Average dose differences reached 1% for the largest field size. Lateral profiles presented dose differences well within 1% for fields up to 20 × 20 cm{sup 2}, while the discrepancy increased toward 2% in the 40 × 40 cm{sup 2} cases. Gamma analysis resulted in an agreement of 100% when a 2%, 2 mm criterion was used, with the only exception of the 40 × 40 cm{sup 2} field (∼95% agreement). With the more stringent criteria of 1%, 1 mm, the agreement reduced to almost 95% for field sizes up to 10 × 10 cm{sup 2}, worse for larger fields. Unflatness and slope FFF-specific parameters are in line with the possible energy underestimation of the simulated results relative to experimental data. Conclusions: The agreement between Monte Carlo simulations and experimental data proved that the evaluated Varian phase-space files for FFF beams from TrueBeam can be used as radiation sources for accurate Monte Carlo dose estimation, especially for field sizes up to 10 × 10 cm{sup 2}, that is the range of field sizes mostly used in combination to the FFF, high dose rate beams.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dosimetric evaluation of new approaches in GRID therapy using nonconventional radiation sources

Purpose: Spatial fractionation of the dose has proven to be a promising approach to increase the tolerance of healthy tissue, which is the main limitation of radiotherapy. A good example of that is GRID therapy, which has been successfully used in the management of large tumors with low toxicity. The aim of this work is to explore new avenues using nonconventional sources: GRID therapy by using kilovoltage (synchrotron) x-rays, the use of very high-energy electrons, and proton GRID therapy. They share in common the use of the smallest possible grid sizes in order to exploit the dose–volume effects. Methods: Monte Carlo simulations (PENELOPE/PENEASY and GEANT4/GATE codes) were used as a method to study dose distributions resulting from irradiations in different configurations of the three proposed techniques. As figure of merit, percentage (peak and valley) depth dose curves, penumbras, and central peak-to-valley dose ratios (PVDR) were evaluated. As shown in previous biological experiments, high PVDR values are requested for healthy tissue sparing. A superior tumor control may benefit from a lower PVDR. Results: High PVDR values were obtained in the healthy tissue for the three cases studied. When low energy photons are used, the treatment of deep-seated tumors can still be performed with submillimetric grid sizes. Superior PVDR values were reached with the other two approaches in the first centimeters along the beam path. The use of protons has the advantage of delivering a uniform dose distribution in the tumor, while healthy tissue benefits from the spatial fractionation of the dose. In the three evaluated techniques, there is a net reduction in penumbra with respect to radiosurgery. Conclusions: The high PVDR values in the healthy tissue and the use of small grid sizes in the three presented approaches might constitute a promising alternative to treat tumors with such spatially fractionated radiotherapy techniques. The dosimetric results presented here support the interest of performing radiobiology experiments in order to evaluate these new avenues.

07 ISOTOPE AND RADIATION SOURCES↗

Absolute x-ray dosimetry on a synchrotron medical beam line with a graphite calorimeter

Purpose: The absolute dose rate of the Imaging and Medical Beamline (IMBL) on the Australian Synchrotron was measured with a graphite calorimeter. The calorimetry results were compared to measurements from the existing free-air chamber, to provide a robust determination of the absolute dose in the synchrotron beam and provide confidence in the first implementation of a graphite calorimeter on a synchrotron medical beam line. Methods: The graphite calorimeter has a core which rises in temperature when irradiated by the beam. A collimated x-ray beam from the synchrotron with well-defined edges was used to partially irradiate the core. Two filtration sets were used, one corresponding to an average beam energy of about 80 keV, with dose rate about 50 Gy/s, and the second filtration set corresponding to average beam energy of 90 keV, with dose rate about 20 Gy/s. The temperature rise from this beam was measured by a calibrated thermistor embedded in the core which was then converted to absorbed dose to graphite by multiplying the rise in temperature by the specific heat capacity for graphite and the ratio of cross-sectional areas of the core and beam. Conversion of the measured absorbed dose to graphite to absorbed dose to water was achieved using Monte Carlo calculations with the EGSnrc code. The air kerma measurements from the free-air chamber were converted to absorbed dose to water using the AAPM TG-61 protocol. Results: Absolute measurements of the IMBL dose rate were made using the graphite calorimeter and compared to measurements with the free-air chamber. The measurements were at three different depths in graphite and two different filtrations. The calorimetry measurements at depths in graphite show agreement within 1% with free-air chamber measurements, when converted to absorbed dose to water. The calorimetry at the surface and free-air chamber results show agreement of order 3% when converted to absorbed dose to water. The combined standard uncertainty is 3.9%. Conclusions: The good agreement of the graphite calorimeter and free-air chamber results indicates that both devices are performing as expected. Further investigations at higher dose rates than 50 Gy/s are planned. At higher dose rates, recombination effects for the free-air chamber are much higher and expected to lead to much larger uncertainties. Since the graphite calorimeter does not have problems associated with dose rate, it is an appropriate primary standard detector for the synchrotron IMBL x rays and is the more accurate dosimeter for the higher dose rates expected in radiotherapy applications.

61 RADIATION PROTECTION AND DOSIMETRY↗