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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Fundamental radiological and geometric performance of two types of proton beam modulated discrete scanning systems

Purpose: The purpose of this investigation was to compare and contrast the measured fundamental properties of two new types of modulated proton scanning systems. This provides a basis for clinical expectations based on the scanned beam quality and a benchmark for computational models. Because the relatively small beam and fast scanning gave challenges to the characterization, a secondary purpose was to develop and apply new approaches where necessary to do so.Methods: The following performances of the proton scanning systems were investigated: beamlet alignment, static in-air beamlet size and shape, scanned in-air penumbra, scanned fluence map accuracy, geometric alignment of scanning system to isocenter, maximum field size, lateral and longitudinal field uniformity of a 1 l cubic uniform field, output stability over time, gantry angle invariance, monitoring system linearity, and reproducibility. A range of detectors was used: film, ionization chambers, lateral multielement and longitudinal multilayer ionization chambers, and a scintillation screen combined with a digital video camera. Characterization of the scanned fluence maps was performed with a software analysis tool.Results: The resulting measurements and analysis indicated that the two types of delivery systems performed within specification for those aspects investigated. The significant differences were observed between the two types of scanning systems where one type exhibits a smaller spot size and associated penumbra than the other. The differential is minimum at maximum energy and increases inversely with decreasing energy. Additionally, the large spot system showed an increase in dose precision to a static target with layer rescanning whereas the small spot system did not.Conclusions: The measured results from the two types of modulated scanning types of system were consistent with their designs under the conditions tested. The most significant difference between the types of system was their proton spot size and associated resolution, factors of magnetic optics, and vacuum length. The need and benefit of mutielement detectors and high-resolution sensors was also shown. The use of a fluence map analytical software tool was particularly effective in characterizing the dynamic proton energy-layer scanning.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2010

Calculation of Dose Deposition in 3D Voxels by Heavy Ions

The biological response to high-LET radiation is very different from low-LET radiation, and can be partly attributed to the energy deposition by the radiation. Several experiments, notably detection of gamma-H2AX foci by immunofluorescence, has revealed important differences in the nature and in the spatial distribution of double-strand breaks (DSB) induced by low- and high-LET radiations. Many calculations, most of which are based on amorphous track models with radial dose, have been combined with chromosome models to calculate the number and distribution of DSB within nuclei and chromosome aberrations. In this work, the Monte-Carlo track structure simulation code RITRACKS have been used to calculate directly the energy deposition in voxels (3D pixels). A cubic volume of 5 micrometers of side was irradiated by 1) 450 (1)H+ ions of 300 MeV (LET is approximately 0.3 keV/micrometer) and 2) by 1 (56)Fe26+ ion of 1 GeV/amu (LET is approximately 150 keV/micrometer). In both cases, the dose deposited in the volume is approximately 1 Gy. All energy deposition events are recorded and dose is calculated in voxels of 20 micrometers of side. The voxels are then visualized in 3D by using a color scale to represent the intensity of the dose in a voxel. This simple approach has revealed several important points which may help understand experimental observations. In both simulations, voxels which receive low dose are the most numerous, and those corresponding to electron track ends received a dose which is in the higher range. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. The distribution of the voxels shows major differences for the (56)Fe26+ ion. The track structure can still be seen, and voxels with much higher dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and may be responsible for DSB that are more difficult to repair. By applying a threshold on the dose visualization, voxels corresponding to electron track ends are evidenced and the spatial distribution of voxels is very similar to the distribution of DSB observed in gamma H2AX experiments, even if no chromosomes have been included in the simulation. Furthermore, this work has shown that a significant dose is deposited in voxels corresponding to electron track ends. Since some delta-rays from iron ion can travel several millimeters, they may also be of radiobiological importance.

Plante, Ianik↗
DOE OSTI2016

Investigating ion recombination effects in a liquid-filled ionization chamber array used for IMRT QA measurements

Purpose: PTW’s Octavius 1000 SRS array performs IMRT quality assurance (QA) measurements with liquid-filled ionization chambers (LICs) to allow closer detector spacing and higher resolution, compared to air-filled QA devices. However, reduced ion mobility in LICs relative to air leads to increased ion recombination effects and reduced collection efficiencies that are dependent on Linac pulse frequency and pulse dose. These pulse parameters are variable during an IMRT delivery, which affects QA results. In this study, (1) 1000 SRS collection efficiencies were measured as a function of pulse frequency and pulse dose, (2) two methods were developed to correct changes in collection efficiencies during IMRT QA measurements, and the effects of these corrections on QA pass rates were compared. Methods: To obtain collection efficiencies, the OCTAVIUS 1000 SRS was used to measure open fields of varying pulse frequency, pulse dose, and beam energy with results normalized to air-filled chamber measurements. Changes in ratios of 1000 SRS to chamber measured dose were attributed to changing collection efficiencies, which were then correlated to pulse parameters using regression analysis. The usefulness of the derived corrections was then evaluated using 6 MV and 10FFF SBRT RapidArc plans delivered to the OCTAVIUS 4D system using a TrueBeam (Varian Medical Systems) linear accelerator equipped with a high definition multileaf collimator. For the first correction, MATLAB software was developed that calculates pulse frequency and pulse dose for each detector, using measurement and DICOM RT Plan files. Pulse information is converted to collection efficiency, and measurements are corrected by multiplying detector dose by ratios of calibration to measured collection efficiencies. For the second correction the MU/min in the daily 1000 SRS calibration was chosen to match the average MU/min of the volumetric modulated arc therapy plan. Effects of the two corrections on QA results were examined by performing 3D gamma analysis comparing predicted to measured dose, with and without corrections. Results: Collection efficiencies correlated linearly to pulse dose, while correlations with pulse frequency were less defined, generally increasing as pulse frequency decreased. After complex MATLAB corrections, average 3D gamma pass rates improved by [0.07%,0.40%,1.17%] for 6 MV and [0.29%,1.40%,4.57%] for 10FFF using [3%/3 mm,2%/2 mm,1%/1 mm] criteria. Maximum changes in gamma pass rates were [0.43%,1.63%,3.05%] for 6 MV and [1.00%,4.80%,11.2%] for 10FFF using [3%/3 mm,2%/2 mm,1%/1 mm] criteria. On average, pass rates of simple daily calibration corrections were within 1% of complex MATLAB corrections. Conclusions: OCTAVIUS 1000 SRS ion recombination effects have little effect on 6 MV measurements. However, the effect could potentially be clinically significant for higher pulse dose unflattened beams when using tighter gamma tolerances, especially when small aperture sizes are used, as is common for SRS/SBRT. In addition, ion recombination effects are strongly correlated to changing MU/min, therefore MU/min used in daily 1000 SRS calibrations should be matched to the expected average MU/min of the IMRT plan.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Validation of a Monte Carlo model for the MLC-Elekta Precise with MCNP6 - Paper 76

Cancer treatment using radiotherapy techniques needs an accurate prediction of the dose delivered to the patient at each position. Treatment planning systems (TPS) allow the optimization of the dose delivered to the tumor volume while minimizing the dose absorbed by the surrounding healthy tissue. Commercial treatment planning systems use empirical deterministic algorithms that do not correctly estimate the dose in the surroundings of heterogeneities. This lack of accuracy is taken into account as bigger security margins in the treatment, reducing the potential effect of the treatment. However, Monte Carlo (MC) simulation techniques have proved to predict the dose distribution more accurately and also representing experimental data better. Nevertheless, the clinical application of these techniques has been limited by the long calculation times they require. As high performance computers decrease in cost MC simulations become more feasible to be used in TPS. The present work improves a simplified model of the MLC-Elekta Precise linear accelerator (linac or unit from now on) previously developed and correctly validated by the group. A multileaf collimator (MLC) has been modelled according to the manufacturer specifications and added to the previous model to allow higher conformity in the studies. The full geometry has been entered to the MCNP6 Monte Carlo code to simulate the transport of particles through the linac. The geometry also includes a cubic water phantom for dosimetry studies and validation with experimental data. The experimental data employed for the dose distribution comparison has been measured at the linac from the linac Hospital Universitari Clinic de Valencia (Spain). Additionally, particles are stored in a phase space surface before their entrance in the water to be analyzed and for performance optimization. Given the results presented, we conclude that this first approach to validate a model of the Elekta Precise linac from the Hospital Cl nic Universitari de Valancia represents well the depth-dose distribution, but fails at representing the dose profiles. It has been shown that they predict too much penumbra, probably caused by the filtration between the leaves. In order to solve this problems, further work is needed to improve the simulated model. This work will mainly focus on the multileaf collimator. Also, the codes developed by the group to analyze phase space files produced by MCNP6 have proven highly useful. As work progress and the codes are debugged, they will be published so that the community can take advantage of them.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Computer program for the sensitivity calculation of a CR-39 detector in a diffusion chamber for radon measurements

Computer software for calculation of the sensitivity of a CR-39 detector closed in a diffusion chamber to radon is described in this work. The software consists of two programs, both written in the standard Fortran 90 programming language. The physical background and a numerical example are given. Presented software is intended for numerous researches in radon measurement community. Previously published computer programs TRACK-TEST.F90 and TRACK-VISION.F90 [D. Nikezic and K. N. Yu, Comput. Phys. Commun. 174, 160 (2006); D. Nikezic and K. N. Yu, Comput. Phys. Commun. 178, 591 (2008)] are used here as subroutines to calculate the track parameters and to determine whether the track is visible or not, based on the incident angle, impact energy, etching conditions, gray level, and visibility criterion. The results obtained by the software, using five different V functions, were compared with the experimental data found in the literature. Application of two functions in this software reproduced experimental data very well, while other three gave lower sensitivity than experiment.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Feature-space assessment of electrical impedance tomography coregistered with computed tomography in detecting multiple contrast targets

Purpose: Fusion of electrical impedance tomography (EIT) with computed tomography (CT) can be useful as a clinical tool for providing additional physiological information about tissues, but requires suitable fusion algorithms and validation procedures. This work explores the feasibility of fusing EIT and CT images using an algorithm for coregistration. The imaging performance is validated through feature space assessment on phantom contrast targets. Methods: EIT data were acquired by scanning a phantom using a circuit, configured for injecting current through 16 electrodes, placed around the phantom. A conductivity image of the phantom was obtained from the data using electrical impedance and diffuse optical tomography reconstruction software (EIDORS). A CT image of the phantom was also acquired. The EIT and CT images were fused using a region of interest (ROI) coregistration fusion algorithm. Phantom imaging experiments were carried out on objects of different contrasts, sizes, and positions. The conductive medium of the phantoms was made of a tissue-mimicking bolus material that is routinely used in clinical radiation therapy settings. To validate the imaging performance in detecting different contrasts, the ROI of the phantom was filled with distilled water and normal saline. Spatially separated cylindrical objects of different sizes were used for validating the imaging performance in multiple target detection. Analyses of the CT, EIT and the EIT/CT phantom images were carried out based on the variations of contrast, correlation, energy, and homogeneity, using a gray level co-occurrence matrix (GLCM). A reference image of the phantom was simulated using EIDORS, and the performances of the CT and EIT imaging systems were evaluated and compared against the performance of the EIT/CT system using various feature metrics, detectability, and structural similarity index measures. Results: In detecting distilled and normal saline water in bolus medium, EIT as a stand-alone imaging system showed contrast discrimination of 47%, while the CT imaging system showed a discrimination of only 1.5%. The structural similarity index measure showed a drop of 24% with EIT imaging compared to CT imaging. The average detectability measure for CT imaging was found to be 2.375 ± 0.19 before fusion. After complementing with EIT information, the detectability measure increased to 11.06 ± 2.04. Based on the feature metrics, the functional imaging quality of CT and EIT were found to be 2.29% and 86%, respectively, before fusion. Structural imaging quality was found to be 66% for CT and 16% for EIT. After fusion, functional imaging quality improved in CT imaging from 2.29% to 42% and the structural imaging quality of EIT imaging changed from 16% to 66%. The improvement in image quality was also observed in detecting objects of different sizes. Conclusions: The authors found a significant improvement in the contrast detectability performance of CT imaging when complemented with functional imaging information from EIT. Along with the feature assessment metrics, the concept of complementing CT with EIT imaging can lead to an EIT/CT imaging modality which might fully utilize the functional imaging abilities of EIT imaging, thereby enhancing the quality of care in the areas of cancer diagnosis and radiotherapy treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Microdosimetry of proton and carbon ions

Purpose: To investigate microdosimetry properties of 160 MeV/u protons and 290 MeV/u{sup 12}C ion beams in small volumes of diameters 10–100 nm. Methods: Energy distributions of primary particles and nuclear fragments in the beams were calculated from simulations with the general purpose code SHIELD-HIT, while energy depositions by monoenergetic ions in nanometer volumes were obtained from the event-by-event Monte Carlo track structure ion code PITS99 coupled with the electron track structure code KURBUC. Results: The results are presented for frequencies of energy depositions in cylindrical targets of diameters 10–100 nm, dose distributionsyd(y) in lineal energy y, and dose-mean lineal energies y{sup ¯}{sub D}. For monoenergetic ions, the y{sup ¯}{sub D} was found to increase with an increasing target size for high-linear energy transfer (LET) ions, but decrease with an increasing target size for low-LET ions. Compared to the depth dose profile of the ion beams, the maximum of the y{sup ¯}{sub D} depth profile for the 160 MeV proton beam was located at ∼0.5 cm behind the Bragg peak maximum, while the y{sup ¯}{sub D} peak of the 290 MeV/u {sup 12}C beam coincided well with the peak of the absorbed dose profile. Differences between the y{sup ¯}{sub D} and dose-averaged linear energy transfer (LET{sub D}) were large in the proton beam for both target volumes studied, and in the {sup 12}C beam for the 10 nm diameter cylindrical volumes. The y{sup ¯}{sub D} determined for 100 nm diameter cylindrical volumes in the {sup 12}C beam was approximately equal to the LET{sub D}. The contributions from secondary particles to the y{sup ¯}{sub D} of the beams are presented, including the contributions from secondary protons in the proton beam and from fragments with atomic number Z = 1–6 in the {sup 12}C beam. Conclusions: The present investigation provides an insight into differences in energy depositions in subcellular-size volumes when irradiated by proton and carbon ion beams. The results are useful for characterizing ion beams of practical importance for biophysical modeling of radiation-induced DNA damage response and repair in the depth profiles of protons and carbon ions used in radiotherapy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

How does C-VIEW image quality compare with conventional 2D FFDM?

Purpose: The FDA approved the use of digital breast tomosynthesis (DBT) in 2011 as an adjunct to 2D full field digital mammography (FFDM) with the constraint that all DBT acquisitions must be paired with a 2D image to assure adequate interpretative information is provided. Recently manufacturers have developed methods to provide a synthesized 2D image generated from the DBT data with the hope of sparing patients the radiation exposure from the FFDM acquisition. While this much needed alternative effectively reduces the total radiation burden, differences in image quality must also be considered. The goal of this study was to compare the intrinsic image quality of synthesized 2D C-VIEW and 2D FFDM images in terms of resolution, contrast, and noise. Methods: Two phantoms were utilized in this study: the American College of Radiology mammography accreditation phantom (ACR phantom) and a novel 3D printed anthropomorphic breast phantom. Both phantoms were imaged using a Hologic Selenia Dimensions 3D system. Analysis of the ACR phantom includes both visual inspection and objective automated analysis using in-house software. Analysis of the 3D anthropomorphic phantom includes visual assessment of resolution and Fourier analysis of the noise. Results: Using ACR-defined scoring criteria for the ACR phantom, the FFDM images scored statistically higher than C-VIEW according to both the average observer and automated scores. In addition, between 50% and 70% of C-VIEW images failed to meet the nominal minimum ACR accreditation requirements—primarily due to fiber breaks. Software analysis demonstrated that C-VIEW provided enhanced visualization of medium and large microcalcification objects; however, the benefits diminished for smaller high contrast objects and all low contrast objects. Visual analysis of the anthropomorphic phantom showed a measureable loss of resolution in the C-VIEW image (11 lp/mm FFDM, 5 lp/mm C-VIEW) and loss in detection of small microcalcification objects. Spectral analysis of the anthropomorphic phantom showed higher total noise magnitude in the FFDM image compared with C-VIEW. Whereas the FFDM image contained approximately white noise texture, the C-VIEW image exhibited marked noise reduction at midfrequency and high frequency with far less noise suppression at low frequencies resulting in a mottled noise appearance. Conclusions: Their analysis demonstrates many instances where the C-VIEW image quality differs from FFDM. Compared to FFDM, C-VIEW offers a better depiction of objects of certain size and contrast, but provides poorer overall resolution and noise properties. Based on these findings, the utilization of C-VIEW images in the clinical setting requires careful consideration, especially if considering the discontinuation of FFDM imaging. Not explicitly explored in this study is how the combination of DBT + C-VIEW performs relative to DBT + FFDM or FFDM alone.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Modeling of Pyroprocessing Hot Cell for Process Monitoring Evaluation

In this report, we examine studies undertaken to determine the applicability of applying measurements with a neutron counter within a pyroprocessing hot cell. The studies focus on the modeling of the hot cell and neutron counter in Monte Carlo N- Particle Code (MCNP) to determine expected radiation effects and applicability of neutron counters. The neutron counter of focus was the High Dose Neutron Detector (HDND) developed at Los Alamos National Laboratory (LANL). We studied the HDND to determine what would be the detector response under normal background conditions with operating process equipment as well as the background’s potential effect on counting applications involving a U/TRU ingot produced by the pyroprocessing electrorefiner. To perform this task, we were required to develop models of the pyroprocessing process cell and unit operations and generate source terms for various different times in the process. The data for the pyroprocessing materials came from SNL’s Safeguards and Security Process Model (SSPM) as part of MPACT’s advanced integration mission. This report will discuss the development of the pyroprocessing process cell and unit operations models in MCNP as well as various simulations to determine the applicability of the HDND for process monitoring efforts in the hot cell.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2017

Materials Degradation in the Jovian Radiation Environment

The radiation environment of Jupiter represents a significant hazard for Europa Lander deorbit stage components, and presents a significant potential mission risk. The radiolytic degradation of ammonium perchlorate (AP) oxidizer in solid propellants may affect its properties and performance. The Monte Carlo code MONSOL was used for modeling of laboratory experiments on the electron irradiation of propellant samples. An approach for flattening dose profiles along the depth of irradiated samples is proposed. Depth-dose distributions produced by Jovian electrons in multi-layer slabs of materials are calculated. It is found that the absorbed dose in a particular slab is significantly affected by backscattered electrons and photons from neighboring slabs. The dose and radiolytic decomposition of AP crystals are investigated and radiation-induced chemical yields and weight percent of radical products are reported.

Miloshevsky, Gennady↗
DOE OSTI2011

Facility for Advanced Accelerator Experimental Tests (FACET) at SLAC and its Radiological Considerations

Facility for Advanced Accelerator Experimental Tests (FACET) in SLAC will be used to study plasma wakefield acceleration. FLUKA Monte Carlo code was used to design a maze wall to separate FACET project and LCLS project to allow persons working in FACET side during LCLS operation. Also FLUKA Monte Carlo code was used to design the shielding for FACET dump to get optimum design for shielding both prompt and residual doses, as well as reducing environmental impact. FACET will be an experimental facility that provides short, intense pulses of electrons and positrons to excite plasma wakefields and study a variety of critical issues associated with plasma wakefield acceleration [1]. This paper describes the FACET beam parameters, the lay-out and its radiological issues.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Monte Carlo simulation for the transport beamline

In the framework of the ELIMED project, Monte Carlo (MC) simulations are widely used to study the physical transport of charged particles generated by laser-target interactions and to preliminarily evaluate fluence and dose distributions. An energy selection system and the experimental setup for the TARANIS laser facility in Belfast (UK) have been already simulated with the GEANT4 (GEometry ANd Tracking) MC toolkit. Preliminary results are reported here. Future developments are planned to implement a MC based 3D treatment planning in order to optimize shots number and dose delivery.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Radiological Modeling for Determination of Derived Concentration Levels of an Area with Uranium Residual Material - 13533

As a result of a pilot project developed at the old Spanish 'Junta de Energia Nuclear' to extract uranium from ores, tailings materials were generated. Most of these residual materials were sent back to different uranium mines, but a small amount of it was mixed with conventional building materials and deposited near the old plant until the surrounding ground was flattened. The affected land is included in an area under institutional control and used as recreational area. At the time of processing, uranium isotopes were separated but other radionuclides of the uranium decay series as Th-230, Ra-226 and daughters remain in the residue. Recently, the analyses of samples taken at different ground's depths confirmed their presence. This paper presents the methodology used to calculate the derived concentration level to ensure that the reference dose level of 0.1 mSv y-1 used as radiological criteria. In this study, a radiological impact assessment was performed modeling the area as recreational scenario. The modelization study was carried out with the code RESRAD considering as exposure pathways, external irradiation, inadvertent ingestion of soil, inhalation of resuspended particles, and inhalation of radon (Rn-222). As result was concluded that, if the concentration of Ra-226 in the first 15 cm of soil is lower than, 0.34 Bq g{sup -1}, the dose would not exceed the reference dose. Applying this value as a derived concentration level and comparing with the results of measurements on the ground, some areas with a concentration of activity slightly higher than latter were found. In these zones the remediation proposal has been to cover with a layer of 15 cm of clean material. This action represents a reduction of 85% of the dose and ensures compliance with the reference dose. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

SU-E-T-150: Brachytherapy QA Employing a High Resolution Liquid Filled Ionisation Chamber Array: Initial Experience and Limitations

Purpose: Verifying a complex 3D brachytherapy dose distribution by measurement is non-trivial. Ideally a photon detector array should be independent of energy and angle, have high spatial resolution and be robust for routine clinical use. Methods: An iridium-192 source was used. A PMMA jig was constructed comprising an outer slab and a central insert with eight milled channels for 1.33mm (outer diameter) steel needles, see figure. All calculations were performed using an empty CT study reconstructing eight virtual needles (QA-CT), using the v2 source model (Elekta AG, Sweden). A high resolution liquid filled ionisation chamber array SRS1000, together with Verisoft software v6.0 (PTW Freiburg, Germany), was used to perform measurements of plans of increasing complexity to evaluate its suitability for device- and patient-specific QA. The dimension of backscatter material was investigated. The patient plan dwell time distribution was entered manually into the QA-CT and the dose distribution was calculated. Results: Our measurements indicate that the array is independent of energy and angle. The resulting measured dose values are linearly interpolated to 2025 values. Shifts of 1mm of the entire needle are readily detectable. Individual dwell position shifts (2.5mm) are also readily measurable. Moreover a dwell time increase of 1 second both in the edge and central region are detectable. Conclusion: The high resolution SRS1000 array is a powerful instrument for brachytherapy QA enabling 977 simultaneous measurements to be performed. Our measurements suggest 60mm of RW3 backscatter material upstream and downstream are sufficient. Local percentage difference analysis is useful for device based QA, normalized relative percentage difference is arguably better for patient specific QA. Automated transfer of patient plan dwell time distribution to the QA plan is required to enable a comprehensive patient QA study to be performed. Moreover the described measurement technique can be used for commissioning new brachytherapy applicators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-A-19A-10: Fast Four Dimensional Monte Carlo Dose Computations for Proton Therapy of Lung Cancer

Purpose: To develop and validate a fast and accurate four dimensional (4D) Monte Carlo (MC) dose computation system for proton therapy of lung cancer and other thoracic and abdominal malignancies in which the delivered dose distributions can be affected by respiratory motion of the patient. Methods: A 4D computer tomography (CT) scan for a lung cancer patient treated with protons in our clinic was used to create a time dependent patient model using our in-house, MCNPX-based Monte Carlo system (“MC{sup 2}”). The beam line configurations for two passively scattered proton beams used in the actual treatment were extracted from the clinical treatment plan and a set of input files was created automatically using MC{sup 2}. A full MC simulation of the beam line was computed using MCNPX and a set of phase space files for each beam was collected at the distal surface of the range compensator. The particles from these phase space files were transported through the 10 voxelized patient models corresponding to the 10 phases of the breathing cycle in the 4DCT, using MCNPX and an accelerated (fast) MC code called “FDC”, developed by us and which is based on the track repeating algorithm. The accuracy of the fast algorithm was assessed by comparing the two time dependent dose distributions. Results: The error of less than 1% in 100% of the voxels in all phases of the breathing cycle was achieved using this method with a speedup of more than 1000 times. Conclusion: The proposed method, which uses full MC to simulate the beam line and the accelerated MC code FDC for the time consuming particle transport inside the complex, time dependent, geometry of the patient shows excellent accuracy together with an extraordinary speed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-I-49: Vendor-Independent, Model-Based Iterative Reconstruction On a Rotating Grid with Coordinate-Descent Optimization for CT Imaging Investigations

Purpose: Iterative reconstruction methods show promise for improving image quality and lowering the dose in helical CT. We aim to develop a novel model-based reconstruction method that offers potential for dose reduction with reasonable computation speed and storage requirements for vendor-independent reconstruction from clinical data on a normal desktop computer. Methods: In 2012, Xu proposed reconstructing on rotating slices to exploit helical symmetry and reduce the storage requirements for the CT system matrix. Inspired by this concept, we have developed a novel reconstruction method incorporating the stored-system-matrix approach together with iterative coordinate-descent (ICD) optimization. A penalized-least-squares objective function with a quadratic penalty term is solved analytically voxel-by-voxel, sequentially iterating along the axial direction first, followed by the transaxial direction. 8 in-plane (transaxial) neighbors are used for the ICD algorithm. The forward problem is modeled via a unique approach that combines the principle of Joseph’s method with trilinear B-spline interpolation to enable accurate reconstruction with low storage requirements. Iterations are accelerated with multi-CPU OpenMP libraries. For preliminary evaluations, we reconstructed (1) a simulated 3D ellipse phantom and (2) an ACR accreditation phantom dataset exported from a clinical scanner (Definition AS, Siemens Healthcare). Image quality was evaluated in the resolution module. Results: Image quality was excellent for the ellipse phantom. For the ACR phantom, image quality was comparable to clinical reconstructions and reconstructions using open-source FreeCT-wFBP software. Also, we did not observe any deleterious impact associated with the utilization of rotating slices. The system matrix storage requirement was only 4.5GB, and reconstruction time was 50 seconds per iteration. Conclusion: Our reconstruction method shows potential for furthering research in low-dose helical CT, in particular as part of our ongoing development of an acquisition/reconstruction pipeline for generating images under a wide range of conditions. Our algorithm will be made available open-source as “FreeCT-ICD”. NIH U01 CA181156; Disclosures (McNitt-Gray): Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Decommissioning for the containment vessel of HTR-PM: Monte Carlo simulation and long-time predictions

A high temperature gas cooled reactor-pebble bed module (HTR-PM) is currently under design and construction in China. Calculations of induced radioactivity in the concrete containment vessel of HTR-PM were performed in order to predict the decommissioning problems. The main activation products and the ambient dose equivalent rate corresponding to 40 years of operations followed by 30 days/1 year/10 years of decay, were calculated by using the Monte Carlo particle transport code FLUKA. Finally, according to the Chinese National Standard and on the basis of the knowledge gained from this study, the classifications of the radioactive solid wastes were discussed, and some suggestions for the working duration of the employees were given. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

SU-E-J-166: Sensitivity of Clinically Relevant Dosimetric Parameters to Contouring Uncertainty During Post Implant Dosimetry of Prostate Permanent Seed Implants

Purpose: There is a strong evidence relating post-implant dosimetry for permanent seed prostate brachytherpy to local control rates. The delineation of the prostate on CT images, however, represents a challenge as it is difficult to confidently identify the prostate borders from soft tissue surrounding it. This study aims at quantifying the sensitivity of clinically relevant dosimetric parameters to prostate contouring uncertainty. Methods: The post-implant CT images and plans for a cohort of 43 patients, who have received I–125 permanent prostate seed implant in our centre, were exported to MIM Symphony LDR brachytherapy treatment planning system (MIM Software Inc., Cleveland, OH). The prostate contours in post-implant CT images were expanded/contracted uniformly for margins of ±1.00mm, ±2.00mm, ±3.00mm, ±4.00mm and ±5.00mm (±0.01mm). The values for V100 and D90 were extracted from Dose Volume Histograms for each contour and compared. Results: The mean value of V100 and D90 was obtained as 92.3±8.4% and 108.4±12.3% respectively (Rx=145Gy). V100 was reduced by −3.2±1.5%, −7.2±3.0%, −12.8±4.0%, −19.0±4.8%, − 25.5±5.4% for expanded contours of prostate with margins of +1mm, +2mm, +3mm, +4mm, and +5mm, respectively, while it was increased by 1.6±1.2%, 2.4±2.4%, 2.7±3.2%, 2.9±4.2%, 2.9±5.1% for the contracted contours. D90 was reduced by −6.9±3.5%, −14.5±6.1%, −23.8±7.1%, − 33.6±8.5%, −40.6±8.7% and increased by 4.1±2.6%, 6.1±5.0%, 7.2±5.7%, 8.1±7.3% and 8.1±7.3% for the same set of contours. Conclusion: Systematic expansion errors of more than 1mm may likely render a plan sub-optimal. Conversely contraction errors may Result in labeling a plan likely as optimal. The use of MRI images to contour the prostate should results in better delineation of prostate organ which increases the predictive value of post-op plans. Since observers tend to overestimate the prostate volume on CT, compared with MRI, the impact of the contouring uncertainty on V100 and D90 fortunately, has a conservative effect of underestimating the prostate coverage.

60 APPLIED LIFE SCIENCES↗