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 253 records · Page 14

Recent Improvement of Measurement Instrumentation to Supervise Nuclear Operations and to Contribute Input Data to 3D Simulation Code - 13289

The CEA has developed many compact characterization tools to follow sensitive operations in a nuclear environment. Usually, these devices are made to carry out radiological inventories, to prepare nuclear interventions or to supervise some special operations. These in situ measurement techniques mainly take place at different stages of clean-up operations and decommissioning projects, but they are also in use to supervise sensitive operations when the nuclear plant is still operating. In addition to this, such tools are often associated with robots to access very highly radioactive areas, and thus can be used in accident situations. Last but not least, the radiological data collected can be entered in 3D calculation codes used to simulate the doses absorbed by workers in real time during operations in a nuclear environment. Faced with these ever-greater needs, nuclear measurement instrumentation always has to involve on-going improvement processes. Firstly, this paper will describe the latest developments and results obtained in both gamma and alpha imaging techniques. The gamma camera has been used by the CEA since the 1990's and several changes have made this device more sensitive, more compact and more competitive for nuclear plant operations. It is used to quickly identify hot spots, locating irradiating sources from 50 keV to 1500 keV. Several examples from a wide field of applications will be presented, together with the very latest developments. The alpha camera is a new camera used to see invisible alpha contamination on several kinds of surfaces. The latest results obtained allow real time supervision of a glove box cleaning operation (for {sup 241}Am contamination). The detection principle as well as the main trials and results obtained will be presented. Secondly, this paper will focus on in situ gamma spectrometry methods developed by the CEA with compact gamma spectrometry probes (CdZnTe, LaBr{sub 3}, NaI, etc.). The radiological data collected is used to quantify the activity of hot spots and can also then be entered in 3D models of nuclear plants to simulate intervention scenarios. Recent developments and results will be presented regarding this. Finally, thanks to a large amount of feedback, the interest of using complementary measurements will be discussed. In fact, the recent use of 3D simulation codes requires very accurate knowledge of nuclear plant radiological data. The use of coupled devices such as imaging devices, (gamma and alpha cameras), gamma spectrometry, dose rate mapping, collimated / un-collimated measurements and many other physical values gives an approach to the radiological knowledge of a process or plant with the lowest possible uncertainty. In line with this, the paper will conclude with the future developments and trials that could be assessed in that field of application. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with GEANT 4

Purpose: This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and GEANT4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field. Methods: A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm{sup 2}. Dose scoring was performed using a 3D grid comprising 1 mm{sup 3} voxels. The production thresholds were equivalent for both codes. Results were analyzed based upon a voxel by voxel dose difference between the two codes and also using a volumetric gamma analysis. Results: Comparisons were drawn from central axis depth doses, cross beam profiles, and isodose contours. Both in the presence and absence of a 1.5 T static magnetic field the relative differences in doses scored along the beam central axis were less than 1% for the homogeneous water phantom and all results matched within a maximum of ±2% for heterogeneous phantoms. Volumetric gamma analysis indicated that more than 99% of the examined volume passed gamma criteria of 2%—2 mm (dose difference and distance to agreement, respectively). These criteria were chosen because the minimum primary statistical uncertainty in dose scoring voxels was 0.5%. The presence of the magnetic field affects the dose at the interface depending upon the density of the material on either sides of the interface. This effect varies with the field size. For example, at the water-lung interface a 33.94% increase in dose was observed (relative to the D{sub max}), by both GPUMCD and GEANT4 for the field size of 2 × 2 cm{sup 2} (compared to no B-field case), which increased to 47.83% for the field size of 5 × 5 cm{sup 2} in the presence of the magnetic field. Similarly, at the lung-water interface, the dose decreased by 19.21% (relative to D{sub max}) for a field size of 2 × 2 cm{sup 2} and by 30.01% for 5 × 5 cm{sup 2} field size. For more complex combinations of materials the dose deposition also becomes more complex. Conclusions: The GPUMCD algorithm showed good agreement against GEANT4 both in the presence and absence of a 1.5 T external magnetic field. The application of 1.5 T magnetic field significantly alters the dose at the interfaces by either increasing or decreasing the dose depending upon the density of the material on either side of the interfaces.

61 RADIATION PROTECTION AND DOSIMETRY↗

ALTERNATIVES OF MACCS2 IN LANL DISPERSION ANALYSIS FOR ONSITE AND OFFSITE DOSES

In modeling atmospheric dispersion to determine accidental release of radiological material, one of the common statistical analysis tools used at Los Alamos National Laboratory (LANL) is MELCOR Accident Consequence Code System, Version 2 (MACCS2). MACCS2, however, has some limitations and shortfalls for both onsite and offsite applications. Alternative computer codes, which could provide more realistic calculations, are being investigated for use at LANL. In the Yucca Mountain Project (YMP), the suitability of MACCS2 for the calculation of onsite worker doses was a concern; therefore, ARCON96 was chosen to replace MACCS2. YMP's use of ARCON96 provided results which clearly demonstrated the program's merit for onsite worker safety analyses in a wide range of complex configurations and scenarios. For offsite public exposures, the conservatism of MACCS2 on the treatment of turbulence phenomena at LANL is examined in this paper. The results show a factor of at least two conservatism in calculated public doses. The new EPA air quality model, AERMOD, which implements advanced meteorological turbulence calculations, is a good candidate for LANL applications to provide more confidence in the accuracy of offsite public dose projections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The alanine detector in BNCT dosimetry: Dose response in thermal and epithermal neutron fields

Purpose: The response of alanine solid state dosimeters to ionizing radiation strongly depends on particle type and energy. Due to nuclear interactions, neutron fields usually also consist of secondary particles such as photons and protons of diverse energies. Various experiments have been carried out in three different neutron beams to explore the alanine dose response behavior and to validate model predictions. Additionally, application in medical neutron fields for boron neutron capture therapy is discussed. Methods: Alanine detectors have been irradiated in the thermal neutron field of the research reactor TRIGA Mainz, Germany, in five experimental conditions, generating different secondary particle spectra. Further irradiations have been made in the epithermal neutron beams at the research reactors FiR 1 in Helsinki, Finland, and Tsing Hua open pool reactor in HsinChu, Taiwan ROC. Readout has been performed with electron spin resonance spectrometry with reference to an absorbed dose standard in a {sup 60}Co gamma ray beam. Absorbed doses and dose components have been calculated using the Monte Carlo codes FLUKA and MCNP. The relative effectiveness (RE), linking absorbed dose and detector response, has been calculated using the Hansen and Olsen alanine response model. Results: The measured dose response of the alanine detector in the different experiments has been evaluated and compared to model predictions. Therefore, a relative effectiveness has been calculated for each dose component, accounting for its dependence on particle type and energy. Agreement within 5% between model and measurement has been achieved for most irradiated detectors. Significant differences have been observed in response behavior between thermal and epithermal neutron fields, especially regarding dose composition and depth dose curves. The calculated dose components could be verified with the experimental results in the different primary and secondary particle fields. Conclusions: The alanine detector can be used without difficulty in neutron fields. The response has been understood with the model used which includes the relative effectiveness. Results and the corresponding discussion lead to the conclusion that application in neutron fields for medical purpose is limited by its sensitivity but that it is a useful tool as supplement to other detectors and verification of neutron source descriptions.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-296: Single Field Per Day Vs. Multiple Fields Per Day and the Impact On BED in Proton Therapy Treatment

Purpose: A common practice, in proton therapy, is to deliver a rotating subset of fields from the treatment plan for the daily fractions. This study compares the impact this practice has on the biological effective dose (BED) versus delivering all planned fields daily. Methods: For two scenarios (a phantom with a geometry approximating the anatomy of a prostate treatment with opposing lateral beams, and a clinical 3-field brain treatment), treatment plans were produced in Eclipse (Varian) to simulate delivery of one, two, and three fields per fraction. The RT-Dose file, structure set, and α/β ratios were processed using in-house MATLAB code to return a new RT-Dose file containing the BED (including a proton RBE of 1.1) which was imported back into Eclipse for analysis. Results: For targets and regions of field overlap in the treatment plan, BED is not affected by delivery regimen. In the phantom, BED in the femoral heads showed increased by 20% when a single field was used rather than two fields. In the brain treatment, the minimum BED to the left optic nerve and the pituitary gland increased by 13% and 10% respectively, for a one-field regime compared to three-fields per fraction. Comparing the two-field and threefield regimes, the optic nerve BED was not significantly affected and the minimum pituitary BED was 4% higher for two fields per day. Conclusion: Hypo-fractionation effects, in regions of non-overlap of fields, significantly increase the BED to the involved tissues by as much as 20%. Care should be taken to avoid inadvertently sacrificing plan effectiveness in the interest of reduced treatment time.

60 APPLIED LIFE SCIENCES↗

Dose-to-water conversion for the backscatter-shielded EPID: A frame-based method to correct for EPID energy response to MLC transmitted radiation

Purpose: To develop a frame-by-frame correction for the energy response of amorphous silicon electronic portal imaging devices (a-Si EPIDs) to radiation that has transmitted through the multileaf collimator (MLC) and to integrate this correction into the backscatter shielded EPID (BSS-EPID) dose-to-water conversion model. Methods: Individual EPID frames were acquired using a Varian frame grabber and iTools acquisition software then processed using in-house software developed inMATLAB. For each EPID image frame, the region below the MLC leaves was identified and all pixels in this region were multiplied by a factor of 1.3 to correct for the under-response of the imager to MLC transmitted radiation. The corrected frames were then summed to form a corrected integrated EPID image. This correction was implemented as an initial step in the BSS-EPID dose-to-water conversion model which was then used to compute dose planes in a water phantom for 35 IMRT fields. The calculated dose planes, with and without the proposed MLC transmission correction, were compared to measurements in solid water using a two-dimensional diode array. Results: It was observed that the integration of the MLC transmission correction into the BSS-EPID dose model improved agreement between modeled and measured dose planes. In particular, the MLC correction produced higher pass rates for almost all Head and Neck fields tested, yielding an average pass rate of 99.8% for 2%/2 mm criteria. A two-sample independentt-test and fisher F-test were used to show that the MLC transmission correction resulted in a statistically significant reduction in the mean and the standard deviation of the gamma values, respectively, to give a more accurate and consistent dose-to-water conversion. Conclusions: The frame-by-frame MLC transmission response correction was shown to improve the accuracy and reduce the variability of the BSS-EPID dose-to-water conversion model. The correction may be applied as a preprocessing step in any pretreatment portal dosimetry calculation and has been shown to be beneficial for highly modulated IMRT fields.

60 APPLIED LIFE SCIENCES↗

Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions

Purpose: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a {sup 192}Ir source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. Methods: The PENELOPE2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a {sup 192}Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth–dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. Results: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. Conclusions: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.

60 APPLIED LIFE SCIENCES↗

TSD Dose and X/Q calculation

The purposes of this study are: Calculate the radiological dose due to a potential accident for the Transportation Safety Document (TSD). The calculation is based on the DOE Tool Box code MACCS2 Version 1.13.1 and FGR-13 dose conversion factor. Calculate the dilution factor X/Q which will be used to determine the consequence of a potential toxic chemical release. The calculation is based on the DOE Tool Box code MACCS2 Version 1.13.1. Calculate the Equivalent Source Term Mass Multiplier (ESTMM) relative to the 30-year aged fuel grade Pu for the 30-year aged weapon grade Pu, natural Uranium, enriched uranium, plutonium isotopes, uranium isotopes, and several other isotopes. The dose and X/Q are calculated based on the following two TSD potential accident scenarios for the release of 1 Ci of Pu-239 and 1 g of fuel grade plutonium: a ground level spill without building wake and a ground level fire with 5 MW sensible heat and with 20'W x 30'H building wake. The Equivalent Source Term Mass Multiplier (ESTMM) is calculated for various isotope materials. The dose and X/Q at 90 m are also calculated by extrapolation which shows that the dose and X/Q calculated by MACCS2 based on a spatial interval (0, 180 m) could be very conservative.

61 RADIATION PROTECTION AND DOSIMETRY↗

Determination of peripheral underdosage at the lung-tumor interface using Monte Carlo radiation transport calculations

Prediction of dose distributions in close proximity to interfaces is difficult. In the context of radiotherapy of lung tumors, this may affect the minimum dose received by lesions and is particularly important when prescribing dose to covering isodoses. The objective of this work is to quantify underdosage in key regions around a hypothetical target using Monte Carlo dose calculation methods, and to develop a factor for clinical estimation of such underdosage. A systematic set of calculations are undertaken using 2 Monte Carlo radiation transport codes (EGSnrc and GEANT4). Discrepancies in dose are determined for a number of parameters, including beam energy, tumor size, field size, and distance from chest wall. Calculations were performed for 1-mm{sup 3} regions at proximal, distal, and lateral aspects of a spherical tumor, determined for a 6-MV and a 15-MV photon beam. The simulations indicate regions of tumor underdose at the tumor-lung interface. Results are presented as ratios of the dose at key peripheral regions to the dose at the center of the tumor, a point at which the treatment planning system (TPS) predicts the dose more reliably. Comparison with TPS data (pencil-beam convolution) indicates such underdosage would not have been predicted accurately in the clinic. We define a dose reduction factor (DRF) as the average of the dose in the periphery in the 6 cardinal directions divided by the central dose in the target, the mean of which is 0.97 and 0.95 for a 6-MV and 15-MV beam, respectively. The DRF can assist clinicians in the estimation of the magnitude of potential discrepancies between prescribed and delivered dose distributions as a function of tumor size and location. Calculation for a systematic set of 'generic' tumors allows application to many classes of patient case, and is particularly useful for interpreting clinical trial data.

61 RADIATION PROTECTION AND DOSIMETRY↗

User's Guide for RESRAD-OFFSITE

The RESRAD-OFFSITE code can be used to model the radiological dose or risk to an offsite receptor. This User’s Guide for RESRAD-OFFSITE Version 3.1 is an update of the User’s Guide for RESRAD-OFFSITE Version 2 contained in the Appendix A of the User’s Manual for RESRAD-OFFSITE Version 2 (ANL/EVS/TM/07-1, DOE/HS-0005, NUREG/CR-6937). This user’s guide presents the basic information necessary to use Version 3.1 of the code. It also points to the help file and other documents that provide more detailed information about the inputs, the input forms and features/tools in the code; two of the features (overriding the source term and computing area factors) are discussed in the appendices to this guide. Section 2 describes how to download and install the code and then verify the installation of the code. Section 3 shows ways to navigate through the input screens to simulate various exposure scenarios and to view the results in graphics and text reports. Section 4 has screen shots of each input form in the code and provides basic information about each parameter to increase the user’s understanding of the code. Section 5 outlines the contents of all the text reports and the graphical output. It also describes the commands in the two output viewers. Section 6 deals with the probabilistic and sensitivity analysis tools available in the code. Section 7 details the various ways of obtaining help in the code.

61 RADIATION PROTECTION AND DOSIMETRY↗

Exploring {sup 57}Co as a new isotope for brachytherapy applications

Purpose: The characteristics of the radionuclide {sup 57}Co make it interesting for use as a brachytherapy source. {sup 57}Co combines a possible high specific activity with the emission of relatively low-energy photons and a half-life (272 days) suitable for regular source exchanges in an afterloader. {sup 57}Co decays by electron capture to the stable {sup 57}Fe with emission of 136 and 122 keV photons. Methods: A hypothetical {sup 57}Co source based on the Flexisource brachytherapy encapsulation with the active core set as a pure cobalt cylinder (length 3.5 mm and diameter 0.6 mm) covered with a cylindrical stainless-steel capsule (length 5 mm and thickness 0.125 mm) was simulated using Geant4 Monte Carlo (MC) code version 9.4. The radial dose function, g(r), and anisotropy function F(r,{theta}), for the line source approximation were calculated following the TG-43U1 formalism. The results were compared to well-known {sup 192}Ir and {sup 125}I radionuclides, representing the higher and the lower energy end of brachytherapy, respectively. Results: The mean energy of photons in water, after passing through the core and the encapsulation material was 123 keV. This hypothetical {sup 57}Co source has an increasing g(r) due to multiple scatter of low-energy photons, which results in a more uniform dose distribution than {sup 192}Ir. Conclusions: {sup 57}Co has many advantages compared to {sup 192}Ir due to its low-energy gamma emissions without any electron contamination. {sup 57}Co has an increasing g(r) that results in a more uniform dose distribution than {sup 192}Ir due to its multiple scattered photons. The anisotropy of the {sup 57}Co source is comparable to that of {sup 192}Ir. Furthermore, {sup 57}Co has lower shielding requirements than {sup 192}Ir.

07 ISOTOPE AND RADIATION SOURCES↗

Feasibility of using intermediate x-ray energies for highly conformal extracranial radiotherapy

Purpose: To investigate the feasibility of using intermediate energy 2 MV x-rays for extracranial robotic intensity modulated radiation therapy. Methods: Two megavolts flattening filter free x-rays were simulated using the Monte Carlo code MCNP (v4c). A convolution/superposition dose calculation program was tuned to match the Monte Carlo calculation. The modeled 2 MV x-rays and actual 6 MV flattened x-rays from existing Varian Linacs were used in integrated beam orientation and fluence optimization for a head and neck, a liver, a lung, and a partial breast treatment. A column generation algorithm was used for the intensity modulation and beam orientation optimization. Identical optimization parameters were applied in three different planning modes for each site: 2, 6 MV, and dual energy 2/6 MV. Results: Excellent agreement was observed between the convolution/superposition and the Monte Carlo calculated percent depth dose profiles. For the patient plans, overall, the 2/6 MV x-ray plans had the best dosimetry followed by 2 MV only and 6 MV only plans. Between the two single energy plans, the PTV coverage was equivalent but 2 MV x-rays improved organs-at-risk sparing. For the head and neck case, the 2MV plan reduced lips, mandible, tongue, oral cavity, brain, larynx, left and right parotid gland mean doses by 14%, 8%, 4%, 14%, 24%, 6%, 30% and 16%, respectively. For the liver case, the 2 MV plan reduced the liver and body mean doses by 17% and 18%, respectively. For the lung case, lung V20, V10, and V5 were reduced by 13%, 25%, and 30%, respectively. V10 of heart with 2 MV plan was reduced by 59%. For the partial breast treatment, the 2 MV plan reduced the mean dose to the ipsilateral and contralateral lungs by 27% and 47%, respectively. The mean body dose was reduced by 16%. Conclusions: The authors showed the feasibility of using flattening filter free 2 MV x-rays for extracranial treatments as evidenced by equivalent or superior dosimetry compared to 6 MV plans using the same inverse noncoplanar intensity modulated planning method.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-J-93: Development of Pre-Contoured Human Model Library in DICOM-RT Format for the Epidemiological Study of the Radiotherapy Patients

Purpose: Prior to 3D conformal radiation therapy planning, patient anatomy information was mostly limited to 2D beams-eye-view from the conventional simulator. To analyze the outcomes of such treatments for radiation late effects, 3D computational human models are often used in commercial treatment planning systems (TPSs). However, several underlying difficulties such as time-consuming manual delineation procedures of a large number of structures in the model have always limited its applications. Primary objective of this work was to develop a human model library for the epidemiological study by converting 3D-surface model organs to DICOM-RT format (DICOM-RT structure) using an in-house built software. We converted the ICRP reference human models to DICOM-RT models, which can be readily adopted for various dose calculations. Methods: MATLAB based code were utilized to convert the contour drawings extracted in text-format from the 3D graphic-tool, Rhinoceros into DICOM-RT structure format for 50 different organs of each model using a 16GB dual-core processor. The conversion periods were measured for each DICOM-RT models, and the reconstructed structure volumes were validated against the original 3D-surface models in the TPS. Ten reference hybrid whole-body models (8-pediatric and 2-adults) were automatically processed to create DICOM-RT computational human model library. Results: Mean contour conversion period was found to be 580 (N=2) and 394.5 (N=8) seconds for 50 organs in the adult and pediatric models respectively. A good agreement for large organs (NRMSD <1.0%) and small organs (NRMSD <7.7%) was also observed between the original volumes and corresponding DICOM-RT structure volumes of the organs. Conclusion: The ICRP reference human models were converted into DICOM-RT format to support the epidemiological study using a large cohort of conventional radiotherapy patients. Due to its DICOM-compatibility, the library may be implemented to many other different applications. We also expect to develop the library by including additional models in future.

60 APPLIED LIFE SCIENCES↗

Position Paper on Practicable Performance Criteria for the Removal Efficiency of Volatile Radionuclides

As a result of fuel reprocessing, volatile radionuclides may be released from the facility stack if no processes are put in place to remove them. The radionuclides that are of concern in this document are 3 H, 14 C, 85 Kr, and 129 I. The question we attempted to answer is how efficient must this removal process be for each of these radionuclides? To answer this question, we examined the three regulations that may impact the degree to which these radionuclides must be reduced before process gases can be released from the facility. These regulations are 40 CFR 61 (EPA 2010a), 40 CFR 190(EPA 2010b), and 10 CFR 20 (NRC 2012). These regulations apply to the total radionuclide release and to a particular organ - the thyroid. Because these doses can be divided amongst all the radionuclides in different ways and even within the four radionuclides in question, we provided several cases. We first looked at the inventories for these radionuclides for three fuel types (PWR UOX, PWR MOX, and AHTGR), several burn-up values, and time out of reactor extending to 200 y. We calculated doses to the maximum exposed individual (MEI) with the EPA code CAP-88 (Rosnick 1992). Finally, we looked at two dose cases. Allocating all of the allowable dose to be used by the volatile radionuclides is one case, but, perhaps, unrealistic. In lieu of this, we arbitrarily selected a value of 10% of the allowable dose to be assigned to the volatile radionuclides. We calculated the required decontamination factors (DFs) for both of these cases, including the case for the thyroid dose for which 14 C and 129 I were the main contributors. With respect to 129 I doses, we found that the highest dose was calculated with iodine as a fine particulate. The dose scaled as the fraction of the total 129 I that was particulate. Therefore, we assumed for all of our calculations that 100% of the 129 I was particulate and allow the user of the results given here to scale our calculated doses to their needs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3D inpatient dose reconstruction from the PET-CT imaging of {sup 90}Y microspheres for metastatic cancer to the liver: Feasibility study

Purpose: The introduction of radioembolization with microspheres represents a significant step forward in the treatment of patients with metastatic disease to the liver. This technique uses semiempirical formulae based on body surface area or liver and target volumes to calculate the required total activity for a given patient. However, this treatment modality lacks extremely important information, which is the three-dimensional (3D) dose delivered by microspheres to different organs after their administration. The absence of this information dramatically limits the clinical efficacy of this modality, specifically the predictive power of the treatment. Therefore, the aim of this study is to develop a 3D dose calculation technique that is based on the PET imaging of the infused microspheres.Methods: The Fluka Monte Carlo code was used to calculate the voxel dose kernel for {sup 90}Y source with voxel size equal to that of the PET scan. The measured PET activity distribution was converted to total activity distribution for the subsequent convolution with the voxel dose kernel to obtain the 3D dose distribution. In addition, dose-volume histograms were generated to analyze the dose to the tumor and critical structures.Results: The 3D inpatient dose distribution can be reconstructed from the PET data of a patient scanned after the infusion of microspheres. A total of seven patients have been analyzed so far using the proposed reconstruction method. Four patients underwent treatment with SIR-Spheres for liver metastases from colorectal cancer and three patients were treated with Therasphere for hepatocellular cancer. A total of 14 target tumors were contoured on post-treatment PET-CT scans for dosimetric evaluation. Mean prescription activity was 1.7 GBq (range: 0.58–3.8 GBq). The resulting mean maximum measured dose to targets was 167 Gy (range: 71–311 Gy). Mean minimum dose to 70% of target (D70) was 68 Gy (range: 25–155 Gy). Mean minimum dose to 90% of target (D90) was 53 Gy (range: 13–125 Gy).Conclusions: A three-dimensional inpatient dose reconstruction method has been developed that is based on the PET/CT data of a patient treated with {sup 90}Y microspheres. It allows for a complete description of the absorbed dose by the tumor and critical structures. It represents the first step in building predictive models for treatment outcomes for patients receiving this therapeutic modality as well as it allows for better analysis of patients' dose response and will ultimately improve future treatment administration.

61 RADIATION PROTECTION AND DOSIMETRY↗

Absolute calibration of the Gamma Knife{sup ®} Perfexion™ and delivered dose verification using EPR/alanine dosimetry

Purpose: Elekta Leksell Gamma Knife{sup ®} (LGK) is a radiotherapy beam machine whose features are not compliant with the international calibration protocols for radiotherapy. In this scope, the Laboratoire National Henri Becquerel and the Pitié-Salpêtrière Hospital decided to conceive a new LKG dose calibration method and to compare it with the currently used one. Furthermore, the accuracy of the dose delivered by the LGK machine was checked using an “end-to-end” test. This study also aims to compare doses delivered by the two latest software versions of the Gammaplan treatment planning system (TPS). Methods: The dosimetric method chosen is the electron paramagnetic resonance (EPR) of alanine. Dose rate (calibration) verification was done without TPS using a spherical phantom. Absolute calibration was done with factors calculated by Monte Carlo simulation (MCNP-X). For “end-to-end” test, irradiations in an anthropomorphic head phantom, close to real treatment conditions, are done using the TPS in order to verify the delivered dose. Results: The comparison of the currently used calibration method with the new one revealed a deviation of +0.8% between the dose rates measured by ion chamber and EPR/alanine. For simple fields configuration (less than 16 mm diameter), the “end-to-end” tests showed out average deviations of −1.7% and −0.9% between the measured dose and the calculated dose by Gammaplan v9 and v10, respectively. Conclusions: This paper shows there is a good agreement between the new calibration method and the currently used one. There is also a good agreement between the calculated and delivered doses especially for Gammaplan v10.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

TH-A-19A-04: Latent Uncertainties and Performance of a GPU-Implemented Pre-Calculated Track Monte Carlo Method

Purpose: Assessing the performance and uncertainty of a pre-calculated Monte Carlo (PMC) algorithm for proton and electron transport running on graphics processing units (GPU). While PMC methods have been described in the past, an explicit quantification of the latent uncertainty arising from recycling a limited number of tracks in the pre-generated track bank is missing from the literature. With a proper uncertainty analysis, an optimal pre-generated track bank size can be selected for a desired dose calculation uncertainty. Methods: Particle tracks were pre-generated for electrons and protons using EGSnrc and GEANT4, respectively. The PMC algorithm for track transport was implemented on the CUDA programming framework. GPU-PMC dose distributions were compared to benchmark dose distributions simulated using general-purpose MC codes in the same conditions. A latent uncertainty analysis was performed by comparing GPUPMC dose values to a “ground truth” benchmark while varying the track bank size and primary particle histories. Results: GPU-PMC dose distributions and benchmark doses were within 1% of each other in voxels with dose greater than 50% of Dmax. In proton calculations, a submillimeter distance-to-agreement error was observed at the Bragg Peak. Latent uncertainty followed a Poisson distribution with the number of tracks per energy (TPE) and a track bank of 20,000 TPE produced a latent uncertainty of approximately 1%. Efficiency analysis showed a 937× and 508× gain over a single processor core running DOSXYZnrc for 16 MeV electrons in water and bone, respectively. Conclusion: The GPU-PMC method can calculate dose distributions for electrons and protons to a statistical uncertainty below 1%. The track bank size necessary to achieve an optimal efficiency can be tuned based on the desired uncertainty. Coupled with a model to calculate dose contributions from uncharged particles, GPU-PMC is a candidate for inverse planning of modulated electron radiotherapy and scanned proton beams. This work was supported in part by FRSQ-MSSS (Grant No. 22090), NSERC RG (Grant No. 432290) and CIHR MOP (Grant No. MOP-211360)

07 ISOTOPE AND RADIATION SOURCES↗

Metallic artifact mitigation and organ-constrained tissue assignment for Monte Carlo calculations of permanent implant lung brachytherapy

Purpose: To investigate methods of generating accurate patient-specific computational phantoms for the Monte Carlo calculation of lung brachytherapy patient dose distributions. Methods: Four metallic artifact mitigation methods are applied to six lung brachytherapy patient computed tomography (CT) images: simple threshold replacement (STR) identifies high CT values in the vicinity of the seeds and replaces them with estimated true values; fan beam virtual sinogram replaces artifact-affected values in a virtual sinogram and performs a filtered back-projection to generate a corrected image; 3D median filter replaces voxel values that differ from the median value in a region of interest surrounding the voxel and then applies a second filter to reduce noise; and a combination of fan beam virtual sinogram and STR. Computational phantoms are generated from artifact-corrected and uncorrected images using several tissue assignment schemes: both lung-contour constrained and unconstrained global schemes are considered. Voxel mass densities are assigned based on voxel CT number or using the nominal tissue mass densities. Dose distributions are calculated using the EGSnrc user-code BrachyDose for{sup 125}I, {sup 103}Pd, and {sup 131}Cs seeds and are compared directly as well as through dose volume histograms and dose metrics for target volumes surrounding surgical sutures. Results: Metallic artifact mitigation techniques vary in ability to reduce artifacts while preserving tissue detail. Notably, images corrected with the fan beam virtual sinogram have reduced artifacts but residual artifacts near sources remain requiring additional use of STR; the 3D median filter removes artifacts but simultaneously removes detail in lung and bone. Doses vary considerably between computational phantoms with the largest differences arising from artifact-affected voxels assigned to bone in the vicinity of the seeds. Consequently, when metallic artifact reduction and constrained tissue assignment within lung contours are employed in generated phantoms, this erroneous assignment is reduced, generally resulting in higher doses. Lung-constrained tissue assignment also results in increased doses in regions of interest due to a reduction in the erroneous assignment of adipose to voxels within lung contours. Differences in dose metrics calculated for different computational phantoms are sensitive to radionuclide photon spectra with the largest differences for{sup 103}Pd seeds and smallest but still considerable differences for {sup 131}Cs seeds. Conclusions: Despite producing differences in CT images, dose metrics calculated using the STR, fan beam + STR, and 3D median filter techniques produce similar dose metrics. Results suggest that the accuracy of dose distributions for permanent implant lung brachytherapy is improved by applying lung-constrained tissue assignment schemes to metallic artifact corrected images.

62 RADIOLOGY AND NUCLEAR MEDICINE↗