Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dose 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 343 records · Page 19

Gamma Dose Measurements in ZPPR-15

Gamma doses were measured in the Zero Power Plutonium Reactor (ZPPR) to provide data for validation of calculations of gamma heating in important reactor components. Stainless steel-encapsulated LiF thermoluminescent dosimeters (TLDs) were irradiated in two configurations in ZPPR-15B and one configuration in ZPPR-15D to provide gamma heating data for the Integral Fast Reactor design in the mid-1980s. The ZPPR-15 experimental records were used to develop detailed as-built Monte Carlo models for the gamma dose measurements in ZPPR-15B and ZPPR-15D. These Monte Carlo models can be converted to deterministic models for calculations with deterministic codes. The reported uncertainty in the published gamma dose measurements is just the statistical uncertainty in the TLD readout. In addition to this statistical uncertainty in the TLD readout, several corrections and other factors contribute to the total uncertainty for the gamma dose measurements. The total uncertainty for the gamma dose measurements can be significantly larger than the published statistical uncertainty for the TLD readout. All known contributions to the total uncertainty were analyzed, and a total uncertainty was derived for the ZPPR-15 gamma dose measurements.

36 MATERIALS SCIENCE↗

Application of Shift Ex-Core Calculations to a Detailed MSR Model

This report documents the collaborative research conducted between Argonne National Laboratory (ANL) and Oak Ridge National Laboratory (ORNL) on applying the Monte Carlo neutronics code Shift to calculate the radiation dose rates for a detailed molten salt reactor (MSR) model based on the Molten Salt Breeder Reactor (MSBR). In this work scope, the ANL and ORNL teams worked together and applied Shift to calculate the radiological environments of this detailed MSR facility modeled with explicit geometries. The radiological conditions within the MSR facility were modeled for when the reactor is at two different operation modes: normal full power and drained state. The FW-CADIS hybrid method in Shift was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Dose rate maps showed that in the current MSBR numerical model, potential pathways exist for neutrons and gammas to stream through the 8-foot concrete shield to reach the top of the MSBR reactor cell. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the Shift simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates, which suggests that a steel liner in the tank model must be added for calculating the dose rates in the drain cell.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SCALE Non-LWR Models for NRC Volume 5

This dataset contains input and result files of computational simulations with the SCALE code system. The simulations cover radionuclide inventory generation, criticality calculations, and dose rate/shielding analyses of various advanced reactors. Users wanting to reproduce results from this dataset are required to obtain a license to the SCALE code system for which details on the distribution can be found here: https://www.ornl.gov/scale/releases

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

AXIOM Unfold 0.7.0, Users Manual

The AXIOM-Unfold application is a computational code for performing spectral unfolds along with uncertainty quantification of the photon spectrum. While this code was principally designed for spectral unfolds on the Saturn source, it is also relevant to other radiation sources such as Pithon. This code is a component of the AXIOM project which was undertaken in order to measure the time-resolved spectrum of the Saturn source; to support this, the AXIOM-Unfold code is able to process time-dependent dose measurements in order to obtain a time-resolved spectrum. This manual contains a full description of the algorithms used by the method. The code features are fully documented along with several worked examples.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Testicular Doses in Image-Guided Radiotherapy of Prostate Cancer

Purpose: To investigate testicular doses contributed by kilovoltage cone-beam computed tomography (kVCBCT) during image-guided radiotherapy (IGRT) of prostate cancer. Methods and Materials: An EGS4 Monte Carlo code was used to calculate three-dimensional dose distributions from kVCBCT on 3 prostate cancer patients. Absorbed doses to various organs were compared between intensity-modulated radiotherapy (IMRT) treatments and kVCBCT scans. The impact of CBCT scanning mode, kilovoltage peak energy (kVp), and CBCT field span on dose deposition to testes and other organs was investigated. Results: In comparison with one 10-MV IMRT treatment, a 125-kV half-fan CBCT scan delivered 3.4, 3.8, 4.1, and 5.7 cGy to the prostate, rectum, bladder, and femoral heads, respectively, accounting for 1.7%, 3.2%, 3.2%, and 8.4% of megavoltage photon dose contributions. However, the testes received 2.9 cGy from the same CBCT scan, a threefold increase as compared with 0.7 cGy received during IMRT. With the same kVp, full-fan mode deposited much less dose to organs than half-fan mode, ranging from 9% less for prostate to 69% less for testes, except for rectum, where full-fan mode delivered 34% more dose. As photon beam energy increased from 60 to 125 kV, kVCBCT-contributed doses increased exponentially for all organs, irrespective of scanning mode. Reducing CBCT field span from 30 to 10 cm in the superior-inferior direction cut testicular doses from 5.7 to 0.2 cGy in half-fan mode and from 1.5 to 0.1 cGy in full-fan mode. Conclusions: Compared with IMRT, kVCBCT-contributed doses to the prostate, rectum, bladder, and femoral heads are clinically insignificant, whereas dose to the testes is threefold more. Full-fan CBCT usually deposits much less dose to organs (except for rectum) than half-fan mode in prostate patients. Kilovoltage CBCT-contributed doses increase exponentially with photon beam energy. Reducing CBCT field significantly cuts doses to testes and other organs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-F-J-99: Dose Accumulation and Evaluation in Lung SBRT Among All Phases of Respiration

Purpose: To calculate the total planning dose on lung tumors (GTV) by accumulating the dose received in all respiration phases. Methods: A patient 4D planning CT (phase-binned, from a Siemens Somatom CT) was used to locate the GTV of a lung tumor in all respiratory phases with Pinnacle (v9.10). GTV contours defined in all phases were projected to the reference phase, where the ITV was defined. Centroids were calculated for all the GTV projections. No deformation or rotation was taken into account. The only GTV contour as defined in the reference phase was voxelized to track each voxel individually. We accumulated the absorbed dose in different phases on each voxel. A 3DCRT and a VMAT plan were designed on the reference phase fulfilling the ITV dosimetric requirements, using the 10MV FFF photon model from an Elekta Versa linac. ITV-to-PTV margins were set to 5mm. In-house developed MATLAB code was used for tumor voxeling and dose accumulation, assuming that the dose distribution planned in the reference phase behaved as a “dose-cloud” during patient breathing. Results: We tested the method on a patient 4DCT set of images exhibiting limited tumor motion (<5mm). For the 3DCRT plan, D95 was calculated for the GTV with motion and for the ITV, showing an agreement of 0.04%. For the VMAT plan, we calculated the D95 for every phase as if the GTV in that phase had received the whole treatment. Differences in D95 for all phases are within 1%, and estimate the potential interplay effect during delivery. Conclusion: A method for dose accumulation and assessment was developed that can compare GTV motion with ITV dosage, and estimate the potential interplay effect for VMAT plans. Work in progress includes the incorporation of deformable image registration and 4D CBCT dose calculation for dose reconstruction and assessment during treatment.

60 APPLIED LIFE SCIENCES↗

Spontaneous wrinkling of soft matter by energetic deposition of Cr and Au

Wrinkling of stiff thin films deposited on compliant substrates is an effect that has been broadly investigated. However, wrinkling consequent to metal ion implantation has been less studied. In the work described here, we have explored the sub-micron wrinkling phenomena that spontaneously occur when metal ions (Au and Cr) are implanted with energy of a few tens of electron volts (49 eV for Au and 72 eV for Cr) into a compliant material (PDMS). This very low energy ion implantation was performed using a Filtered Cathodic Vacuum Arc technique, a process often referred to as energetic deposition or energetic condensation. For comparison, Au and Cr depositions with similar doses were also done using a sputtering technique (with lower particle energy of approximately 2 eV), and no wrinkle formation was then observed. In this way, we can discuss the role of ion energy in wrinkle formation. Depth profiles of the implanted material were calculated using the Tridyn computer simulation code for each metal, for several implantation doses. UV-vis absorption spectroscopy analysis confirmed the presence of metal nanoparticles. Atomic Force Microscopy imaging with spectral processing was used to compare the wrinkle morphology for each case investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU‐C‐105‐05: Reference Dosimetry of High‐Energy Electron Beams with a Farmer‐Type Ionization Chamber

Purpose: To investigate gradient effects and provide Monte Carlo calculated beam quality conversion factors to characterize the Farmer‐type NE2571 ion chamber for high‐energy reference dosimetry of clinical electron beams. Methods: The EGSnrc code system is used to calculate the absorbed dose to water and to the gas in a fully modeled NE2571 chamber as a function of depth in a water phantom. Electron beams incident on the surface of the phantom are modeled using realistic BEAMnrc accelerator simulations and electron beam spectra. Beam quality conversion factors are determined using calculated doses to water and to air in the chamber in high‐energy electron beams and in a cobalt‐60 reference field. Calculated water‐to‐air stopping power ratios are employed for investigation of the overall ion chamber perturbation factor. Results: An upstream shift of 0.3–0.4 multiplied by the chamber radius, r-cav, both minimizes the variation of the overall ion chamber perturbation factor with depth and reduces the difference between the beam quality specifier (R{sub 5} {sub 0}) calculated using ion chamber simulations and that obtained with simulations of dose‐to‐water in the phantom. Beam quality conversion factors are obtained at the reference depth and gradient effects are optimized using a shift of 0.2r-cav. The photon‐electron conversion factor, k-ecal, amounts to 0.906 when gradient effects are minimized using the shift established here and 0.903 if no shift of the data is used. Systematic uncertainties in beam quality conversion factors are investigated and amount to between 0.4 to 1.1% depending on assumptions used. Conclusion: The calculations obtained in this work characterize the use of an NE2571 ion chamber for reference dosimetry of high‐energy electron beams. These results will be useful as the AAPM continues to review their reference dosimetry protocols.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MO-F-CAMPUS-I-04: Characterization of Fan Beam Coded Aperture Coherent Scatter Spectral Imaging Methods for Differentiation of Normal and Neoplastic Breast Structures

Purpose: This study intends to characterize the spectral and spatial resolution limits of various fan beam geometries for differentiation of normal and neoplastic breast structures via coded aperture coherent scatter spectral imaging techniques. In previous studies, pencil beam raster scanning methods using coherent scatter computed tomography and selected volume tomography have yielded excellent results for tumor discrimination. However, these methods don’t readily conform to clinical constraints; primarily prolonged scan times and excessive dose to the patient. Here, we refine a fan beam coded aperture coherent scatter imaging system to characterize the tradeoffs between dose, scan time and image quality for breast tumor discrimination. Methods: An X-ray tube (125kVp, 400mAs) illuminated the sample with collimated fan beams of varying widths (3mm to 25mm). Scatter data was collected via two linear-array energy-sensitive detectors oriented parallel and perpendicular to the beam plane. An iterative reconstruction algorithm yields images of the sample’s spatial distribution and respective spectral data for each location. To model in-vivo tumor analysis, surgically resected breast tumor samples were used in conjunction with lard, which has a form factor comparable to adipose (fat). Results: Quantitative analysis with current setup geometry indicated optimal performance for beams up to 10mm wide, with wider beams producing poorer spatial resolution. Scan time for a fixed volume was reduced by a factor of 6 when scanned with a 10mm fan beam compared to a 1.5mm pencil beam. Conclusion: The study demonstrates the utility of fan beam coherent scatter spectral imaging for differentiation of normal and neoplastic breast tissues has successfully reduced dose and scan times whilst sufficiently preserving spectral and spatial resolution. Future work to alter the coded aperture and detector geometries could potentially allow the use of even wider fans, thereby making coded aperture coherent scatter imaging a clinically viable method for breast cancer detection. United States Department of Homeland Security; Duke University Medical Center - Department of Radiology; Carl E Ravin Advanced Imaging Laboratories; Duke University Medical Physics Graduate Program.

60 APPLIED LIFE SCIENCES↗

Status of the Slide Rule Update - Phase 4 Plutonium Configurations – Delayed Gamma [Slides]

The fourth phase of Slide-Rule update gives a calculation scheme for DFG doses and its application to five plutonium material, as well as comparisons between radiation transport and depletion codes such as MCNP, SCALE, COG and FISPACT. Overall, DFG dose rates calculated by each laboratory led to consistent results. Extra effort is under way to identify the cause of the remaining differences, especially for case 2, by comparing precisely the gamma source, and particularly nuclides inventories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Determining the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5

Purpose: Recent advances in the imaging of {sup 90}Y using positron emission tomography (PET) and improved uncertainty in the branching ratio for the internal pair production component of {sup 90}Y decay allow for a more accurate determination of the activity distribution of {sup 90}Y microspheres within a patient. This improved activity distribution can be convolved with the dose kernel of {sup 90}Y to calculate the dose distribution within a patient. This work investigates the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5 and compares the results of these two transport codes. Methods: Monte Carlo simulations were performed with egsnrc and mcnp5 to calculate the dose rate at multiple radial distances around various {sup 90}Y sources. Point source simulations were completed with mcnp5 to determine the optimal electron transport settings for this work. After determining the optimal settings, point source simulations were completed using egsnrc (user code edknrc) and mcnp5 in water and liver [as defined by the International Commission on Radiation Units and Measurements (ICRU) Report 44]. The results were compared to ICRU Report 72 reference data. Point source simulations were also completed in water with a density of 1.06 g{center_dot}cm{sup -3} to evaluate the effect of the density of the surrounding material. Glass and resin microsphere simulations were performed with average and maximum diameter and density values (based on values given in the literature) in water and in liver. The results were compared to point source simulation results using the same transport code and in the same surrounding material. All simulations had statistical uncertainties less than 1%. Results: The optimal transport settings in mcnp5 for this work included using the energy-and step-specific algorithm (DBCN 17J 2) and ESTEP set to 10. These settings were used for all subsequent simulations with mcnp5. The point source simulations in water for both egsnrc and mcnp5 were found to agree within 2% of the ICRU 72 reference data over the investigated range. Point source simulations in liver had large differences relative to ICRU 72, approaching -60% near the maximum range of {sup 90}Y. These differences are mostly attributed to the difference in density between water (1.0 g{center_dot}cm{sup -3}) and liver (1.06 g{center_dot}cm{sup -3}). Glass and resin microsphere simulations showed a slight decrease in the dose rate near the maximum range of {sup 90}Y relative to the point source simulations. The largest relative differences were approximately -4.2% and -2.8% for the glass and resin microspheres, respectively. Agreement between the egsnrc and mcnp5 simulations results was generally good. Conclusions: The presence of the microsphere material causes slight differences in the {sup 90}Y dose kernel compared to those calculated with point sources. Large differences were seen between simulations in water and those in liver. For the most accurate calculation of the dose distribution, the density of the patient's liver should be accounted for in the calculation of the dose kernel. Lastly, due to the need to determine the optimal transport settings with mcnp5, electron transport with this code should be used with caution.

61 RADIATION PROTECTION AND DOSIMETRY↗

Methodology comparison for gamma-heating calculations in material-testing reactors

The Jules Horowitz Reactor (JHR) is a Material-Testing Reactor (MTR) under construction in the south of France at CEA Cadarache (French Alternative Energies and Atomic Energy Commission). It will typically host about 20 simultaneous irradiation experiments in the core and in the beryllium reflector. These experiments will help us better understand the complex phenomena occurring during the accelerated ageing of materials and the irradiation of nuclear fuels. Gamma heating, i.e. photon energy deposition, is mainly responsible for temperature rise in non-fuelled zones of nuclear reactors, including JHR internal structures and irradiation devices. As temperature is a key parameter for physical models describing the behavior of material, accurate control of temperature, and hence gamma heating, is required in irradiation devices and samples in order to perform an advanced suitable analysis of future experimental results. From a broader point of view, JHR global attractiveness as a MTR depends on its ability to monitor experimental parameters with high accuracy, including gamma heating. Strict control of temperature levels is also necessary in terms of safety. As JHR structures are warmed up by gamma heating, they must be appropriately cooled down to prevent creep deformation or melting. Cooling-power sizing is based on calculated levels of gamma heating in the JHR. Due to these safety concerns, accurate calculation of gamma heating with well-controlled bias and associated uncertainty as low as possible is all the more important. There are two main kinds of calculation bias: bias coming from nuclear data on the one hand and bias coming from physical approximations assumed by computer codes and by general calculation route on the other hand. The former must be determined by comparison between calculation and experimental data; the latter by calculation comparisons between codes and between methodologies. In this presentation, we focus on this latter kind of bias. Nuclear heating is represented by the physical quantity called absorbed dose (energy deposition induced by particle-matter interactions, divided by mass). Its calculation with Monte Carlo codes is possible but computationally expensive as it requires transport simulation of charged particles, along with neutrons and photons. For that reason, the calculation of another physical quantity, called KERMA, is often preferred, as KERMA calculation with Monte Carlo codes only requires transport of neutral particles. However, KERMA is only an estimator of the absorbed dose and many conditions must be fulfilled for KERMA to be equal to absorbed dose, including so-called condition of electronic equilibrium. Also, Monte Carlo computations of absorbed dose still present some physical approximations, even though there is only a limited number of them. Some of these approximations are linked to the way how Monte Carlo codes apprehend the transport simulation of charged particles and the productive and destructive interactions between photons, electrons and positrons. There exists a huge variety of electromagnetic shower models which tackle this topic. Differences in the implementation of these models can lead to discrepancies in calculated values of absorbed dose between different Monte Carlo codes. The magnitude of order of such potential discrepancies should be quantified for JHR gamma-heating calculations. We consequently present a two-pronged plan. In a first phase, we intend to perform compared absorbed dose / KERMA Monte Carlo calculations in the JHR. This way, we will study the presence or absence of electronic equilibrium in the different JHR structures and experimental devices and we will give recommendations for the choice of KERMA or absorbed dose when calculating gamma heating in the JHR. In a second phase, we intend to perform compared TRIPOLI4 / MCNP absorbed dose calculations in a simplified JHR-representative geometry. For this comparison, we will use the same nuclear data library for both codes (the European library JEFF3.1.1 and photon library EPDL97) so as to isolate the effects from electromagnetic shower models on absorbed dose calculation. This way, we hope to get insightful feedback on these models and their implementation in Monte Carlo codes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Fast patient-specific Monte Carlo brachytherapy dose calculations via the correlated sampling variance reduction technique

Purpose: To demonstrate potential of correlated sampling Monte Carlo (CMC) simulation to improve the calculation efficiency for permanent seed brachytherapy (PSB) implants without loss of accuracy. Methods: CMC was implemented within an in-house MC code family (PTRAN) and used to compute 3D dose distributions for two patient cases: a clinical PSB postimplant prostate CT imaging study and a simulated post lumpectomy breast PSB implant planned on a screening dedicated breast cone-beam CT patient exam. CMC tallies the dose difference, {Delta}D, between highly correlated histories in homogeneous and heterogeneous geometries. The heterogeneous geometry histories were derived from photon collisions sampled in a geometrically identical but purely homogeneous medium geometry, by altering their particle weights to correct for bias. The prostate case consisted of 78 Model-6711 {sup 125}I seeds. The breast case consisted of 87 Model-200 {sup 103}Pd seeds embedded around a simulated lumpectomy cavity. Systematic and random errors in CMC were unfolded using low-uncertainty uncorrelated MC (UMC) as the benchmark. CMC efficiency gains, relative to UMC, were computed for all voxels, and the mean was classified in regions that received minimum doses greater than 20%, 50%, and 90% of D{sub 90}, as well as for various anatomical regions. Results: Systematic errors in CMC relative to UMC were less than 0.6% for 99% of the voxels and 0.04% for 100% of the voxels for the prostate and breast cases, respectively. For a 1 x 1 x 1 mm{sup 3} dose grid, efficiency gains were realized in all structures with 38.1- and 59.8-fold average gains within the prostate and breast clinical target volumes (CTVs), respectively. Greater than 99% of the voxels within the prostate and breast CTVs experienced an efficiency gain. Additionally, it was shown that efficiency losses were confined to low dose regions while the largest gains were located where little difference exists between the homogeneous and heterogeneous doses. On an AMD 1090T processor, computing times of 38 and 21 sec were required to achieve an average statistical uncertainty of 2% within the prostate (1 x 1 x 1 mm{sup 3}) and breast (0.67 x 0.67 x 0.8 mm{sup 3}) CTVs, respectively. Conclusions: CMC supports an additional average 38-60 fold improvement in average efficiency relative to conventional uncorrelated MC techniques, although some voxels experience no gain or even efficiency losses. However, for the two investigated case studies, the maximum variance within clinically significant structures was always reduced (on average by a factor of 6) in the therapeutic dose range generally. CMC takes only seconds to produce an accurate, high-resolution, low-uncertainly dose distribution for the low-energy PSB implants investigated in this study.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-E-T-58: Calculation of Dose Distribution of Accuboost Brachytherapy in Deformable Polyvinil Alcohol Breast Phantom Using Biomechanical Modeling and Monte Carlo Simulation

Purpose: the accuboost is the most modern method of breast brachytherapy that is a boost method in compressed tissue by a mammography unit. the dose distribution in uncompressed tissue, as compressed tissue is important that should be characterized. Methods: In this study, the mechanical behavior of breast in mammography loading, the displacement of breast tissue and the dose distribution in compressed and uncompressed tissue, are investigated. Dosimetry was performed by two dosimeter methods of Monte Carlo simulations using MCNP5 code and thermoluminescence dosimeters. For Monte Carlo simulations, the dose values in cubical lattice were calculated using tally F6. The displacement of the breast elements was simulated by Finite element model and calculated using ABAQUS software, from which the 3D dose distribution in uncompressed tissue was determined. The geometry of the model is constructed from MR images of 6 volunteers. Experimental dosimetery was performed by placing the thermoluminescence dosimeters into the polyvinyl alcohol breast equivalent phantom and on the proximal edge of compression plates to the chest. Results: The results indicate that using the cone applicators would deliver more than 95% of dose to the depth of 5 to 17mm, while round applicator will increase the skin dose. Nodal displacement, in presence of gravity and 60N forces, i.e. in mammography compression, was determined with 43% contraction in the loading direction and 37% expansion in orthogonal orientation. Finally, in comparison of the acquired from thermoluminescence dosimeters with MCNP5, they are consistent with each other in breast phantom and in chest's skin with average different percentage of 13.7±5.7 and 7.7±2.3, respectively. Conclusion: The major advantage of this kind of dosimetry is the ability of 3D dose calculation by FE Modeling. Finally, polyvinyl alcohol is a reliable material as a breast tissue equivalent dosimetric phantom that provides the ability of TLD dosimetry for validation.

60 APPLIED LIFE SCIENCES↗

PCP METHODOLOGY FOR DETERMINING DOSE RATES FOR SMALL GRAM QUANTITIES IN SHIPPING PACKAGINGS

The Small Gram Quantity (SGQ) concept is based on the understanding that small amounts of hazardous materials, in this case radioactive materials, are significantly less hazardous than large amounts of the same materials. This study describes a methodology designed to estimate an SGQ for several neutron and gamma emitting isotopes that can be shipped in a package compliant with 10 CFR Part 71 external radiation level limits regulations. These regulations require packaging for the shipment of radioactive materials perform, under both normal and accident conditions, the essential functions of material containment, subcriticality, and maintain external radiation levels within regulatory limits. 10 CFR 71.33(b)(1)(2)&(3) state radioactive and fissile materials must be identified and their maximum quantity, chemical and physical forms be included in an application. Furthermore, the U.S. Federal Regulations require application contain an evaluation demonstrating the package (i.e., the packaging and its contents) satisfies the external radiation standards for all packages (10 CFR 71.31(2), 71.35(a), & 71.47). By placing the contents in a He leak-tight containment vessel, and limiting the mass to ensure subcriticality, the first two essential functions are readily met. Some isotopes emit sufficiently strong photon radiation that small amounts of material can yield a large external dose rate. Quantifying of the dose rate for a proposed content is a challenging issue for the SGQ approach. It is essential to quantify external radiation levels from several common gamma and neutron sources that can be safely placed in a specific packaging, to ensure compliance with federal regulations. The Packaging Certification Program (PCP) Methodology for Determining Dose Rate for Small Gram Quantities in Shipping Packagings described in this report provides bounding mass limits for a set of proposed SGQ isotopes. Methodology calculations were performed to estimate external radiation levels for the 9977 shipping package using the MCNP radiation transport code to develop a set of response multipliers (Green's functions) for 'dose per particle' for each neutron and photon spectral group. The source spectrum for each isotope generated using the ORIGEN-S and RASTA computer codes was folded with the response multipliers to generate the dose rate per gram of each isotope in the 9977 shipping package and its associated shielded containers. The maximum amount of a single isotope that could be shipped within the regulatory limits contained in 10 CFR 71.47 for dose rate at the surface of the package is determined. If a package contains a mixture of isotopes, the acceptability for shipment can be determined by a sum of fractions approach. Furthermore, the results of this analysis can be easily extended to additional radioisotopes by simply evaluating the neutron and/or photon spectra of those isotopes and folding the spectral data with the Green's functions provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Comparison of film measurements and Monte Carlo simulations of dose delivered with very high-energy electron beams in a polystyrene phantom

Purpose: To measure radiation dose in a water-equivalent medium from very high-energy electron (VHEE) beams and make comparisons to Monte Carlo (MC) simulation results. Methods: Dose in a polystyrene phantom delivered by an experimental VHEE beam line was measured with Gafchromic films for three 50 MeV and two 70 MeV Gaussian beams of 4.0–6.9 mm FWHM and compared to corresponding MC-simulated dose distributions. MC dose in the polystyrene phantom was calculated with the EGSnrc/BEAMnrc and DOSXYZnrc codes based on the experimental setup. Additionally, the effect of 2% beam energy measurement uncertainty and possible non-zero beam angular spread on MC dose distributions was evaluated. Results: MC simulated percentage depth dose (PDD) curves agreed with measurements within 4% for all beam sizes at both 50 and 70 MeV VHEE beams. Central axis PDD at 8 cm depth ranged from 14% to 19% for the 5.4–6.9 mm 50 MeV beams and it ranged from 14% to 18% for the 4.0–4.5 mm 70 MeV beams. MC simulated relative beam profiles of regularly shaped Gaussian beams evaluated at depths of 0.64 to 7.46 cm agreed with measurements to within 5%. A 2% beam energy uncertainty and 0.286° beam angular spread corresponded to a maximum 3.0% and 3.8% difference in depth dose curves of the 50 and 70 MeV electron beams, respectively. Absolute dose differences between MC simulations and film measurements of regularly shaped Gaussian beams were between 10% and 42%. Conclusions: The authors demonstrate that relative dose distributions for VHEE beams of 50–70 MeV can be measured with Gafchromic films and modeled with Monte Carlo simulations to an accuracy of 5%. The reported absolute dose differences likely caused by imperfect beam steering and subsequent charge loss revealed the importance of accurate VHEE beam control and diagnostics.

07 ISOTOPE AND RADIATION SOURCES↗

SU-F-T-679: Monte Carlo Simulation On Surface Dose in Preclinical Irradiation Using Monoenergetic Photon Beams

Purpose: This study investigated the surface dose variation in preclinical irradiation using small animal, when monoenergetic photon beams with energy range from 50 keV to 1.25 MeV were used. Methods: Inhomogeneous, homogeneous and bone-tissue homogeneous mouse phantom based on the same CT image set were used. The homogeneous and bone-tissue homogeneous phantom were created with the relative electron density of all and only bone voxels of the mouse overridden to one, respectively. Monte Carlo simulation based on the EGSnrc-based code was used to calculate the surface dose, when the phantoms were irradiated by a 360° photon arc with energies ranging from 50 keV to 1.25 MeV. The mean surface doses of the three phantoms were calculated. In addition, the surface doses from partial arcs, 45°–315°, 125°–225°, 45°–125° and 225°–315° covering the anterior, posterior, right lateral and left lateral region of the mouse were determined using different photon beam energies. Results: When the prescribed dose at the isocenter of the mouse was 2 Gy, the maximum mean surface doses, found at the 50-keV photon beams, were 0.358 Gy, 0.363 Gy and 0.350 Gy for the inhomogeneous, homogeneous and bone-tissue homogeneous mouse phantom, respectively. The mean surface dose of the mouse was found decreasing with an increase of the photon beam energy. For surface dose in different orientations, the lateral regions of the mouse were receiving lower dose than the anterior and posterior regions. This may be due to the increase of beam attenuation along the horizontal (left-right) axis than the vertical (anterior-posterior) in the mouse. Conclusion: It is concluded that consideration of phantom inhomogeneity in the dose calculation resulted in a lower mean surface dose of the mouse. The mean surface dose also decreased with an increase of photon beam energy in the kilovoltage range.

60 APPLIED LIFE SCIENCES↗

MO-E-18C-04: Advanced Computer Simulation and Visualization Tools for Enhanced Understanding of Core Medical Physics Concepts

Purpose: Most medical physics programs emphasize proficiency in routine clinical calculations and QA. The formulaic aspect of these calculations and prescriptive nature of measurement protocols obviate the need to frequently apply basic physical principles, which, therefore, gradually decay away from memory. E.g. few students appreciate the role of electron transport in photon dose, making it difficult to understand key concepts such as dose buildup, electronic disequilibrium effects and Bragg-Gray theory. These conceptual deficiencies manifest when the physicist encounters a new system, requiring knowledge beyond routine activities. Methods: Two interactive computer simulation tools are developed to facilitate deeper learning of physical principles. One is a Monte Carlo code written with a strong educational aspect. The code can “label” regions and interactions to highlight specific aspects of the physics, e.g., certain regions can be designated as “starters” or “crossers,” and any interaction type can be turned on and off. Full 3D tracks with specific portions highlighted further enhance the visualization of radiation transport problems. The second code calculates and displays trajectories of a collection electrons under arbitrary space/time dependent Lorentz force using relativistic kinematics. Results: Using the Monte Carlo code, the student can interactively study photon and electron transport through visualization of dose components, particle tracks, and interaction types. The code can, for instance, be used to study kerma-dose relationship, explore electronic disequilibrium near interfaces, or visualize kernels by using interaction forcing. The electromagnetic simulator enables the student to explore accelerating mechanisms and particle optics in devices such as cyclotrons and linacs. Conclusion: The proposed tools are designed to enhance understanding of abstract concepts by highlighting various aspects of the physics. The simulations serve as virtual experiments that give deeper and long lasting understanding of core principles. The student can then make sound judgements in novel situations encountered beyond routine clinical activities.

60 APPLIED LIFE SCIENCES↗