Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculations”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

PyCrystalField : software for calculation, analysis and fitting of crystal electric field Hamiltonians

PyCrystalField is a Python software package for calculating single-ion crystal electric field (CEF) Hamiltonians. This software can calculate a CEF Hamiltonian ab initio from a point charge model for any transition or rare earth ion in either the J basis or the LS basis, perform symmetry analysis to identify nonzero CEF parameters, calculate the energy spectrum and observables such as neutron spectrum and magnetization, and fit CEF Hamiltonians to any experimental data. Here, the theory, implementation and examples of its use are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A ray-tracing algorithm for ab-initio calculation of thermal load in undulator-based synchrotron beamlines

The OASYS suite and its powerful integration features are used to implement a ray-tracing algorithm to accurately calculate the thermal load in any component of an undulator-based synchrotron beamline. This is achieved by sampling and converting the SRW source of a given energy into a Shadow source and using the latter code to ray trace the full beamline. The accuracy of the algorithm is proved by reconstructing the full undulator radiation distribution through an aperture and comparing the result with direct calculaton of the total power using SRW. The algorithm is particularly suited to analyze cases with complex beamline layouts and optical elements, such as crystals, multilayers, and compound refractive lenses. Examples of its use to calculate the power load on elements of two of the feature beamlines at the Advanced Photon Source Upgrade Project and a comparison of the results with analytical calculations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Measurements of SNAC2 area dosimeters placed in different configurations around the PROSPERO reactor and comparison with TRIPOLI-4 calculations

In the context of criticality accident alarm system tests, several experiments were carried out in 2013 on the PROSPERO reactor to study the response to neutron and gamma of different devices and dosimeters, particularly on the SNAC2 dosimeter. This article presents the results of this criticality dosimeter in different configurations, and compares the experimental measurements with the results of calculation performed with the TRIPOLI-4 Monte-Carlo Neutral Particles transport code. PROSPERO is a metallic critical assembly managed by the Criticality, Neutron Science and Measurement Department located at the French CEA Research Center of Valduc. The core, surrounded by a reflector of depleted uranium, is composed of 2 horizontal cylindrical blocks made of a highly enriched uranium alloy which can be placed in contact, and of 4 depleted uranium control rods which allow the reactor to be driven. This reactor, placed in a cell 10 m x 8 m x 6 m high, with 1.4-meter-thick concrete walls, is used as a fast neutron spectrum source and is operated at stable power level in delayed critical state, which can vary from 3 mW to 3 kW. PROSPERO is extensively used for electronic hardening or to study the effect of the neutrons on various materials. The SNAC2 criticality dosimeter is a zone dosimeter allowing the off line measurement of criticality accident neutron doses. This dosimeter consists of the pile up of seven activation foils embedded into a 23 mm diameter x 21 mm height cadmium container. The activation measurement of each foil, using a gamma spectroscopy technique, gives information about the neutron reaction rates. The SNAC2 software allows the spectrum unfolding from these values, taking into account the hypothesis of a particular spectrum shape, in three components: a Maxwell spectrum component for the thermal range, a 1/E component for the epithermal range, and a Watt spectrum component for the high energy range. Moreover, from the neutron spectrum, the SNAC software can calculate the neutron fluence integrated by the dosimeter and the neutron dose. During the 3 weeks measurement campaign many radioprotection devices were used. To modify the spectrum seen by these devices, several shields of various thicknesses made of concrete or polyethylene, with or without cadmium covers, were placed in the PROSPERO cell. These devices allow the study of criticality accident spectra in several environments: from metal to pseudo liquid. The fluxes measured by the SNAC2 devices were compared with TRIPOLI-4 calculations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Evaluating the Incident Energy of Arcs in Photovoltaic DC Systems: Comparison Between Calculated and Experimental Data

Solar Photovoltaic (PV) systems have permeated the energy generation world at a very high rate, some of the safety codes and standards are still lagging in accurately assessing the hazards and risks associated with PV array arcing energies. Safety professionals and maintenance workers using NFPA 70E have utilized the Doan, Stokes & Oppenlander or Enrique models, meant to determine arc energies in DC power systems using the maximum power method. These methods may lead to an overestimate of energy available in PV systems during a fault. Since PV modules/arrays are non-linear, current limited DC devices, some of these calculation methods may not accurately predict fault energy. This paper will validate current arc energy models for PV systems by comparing experimental and calculated data. Additionally, this data will help modify the current NFPA 70E models related to smaller solar arrays. Understanding where the real safety threshold for DC arc flash in PV systems exists will help maintenance and safety professionals better prepare for a variety of work related activities. This paper will analyze real arc data taken for PV systems <1000VDC and <60amps and compare this to the calculated incident energy models, to include 70E. Thus, using these comparisons, it may be possible to reduce the safety hazard severity and thus relax the PPE requirements for installation and maintenance crews.

14 SOLAR ENERGY↗
DOE OSTI2014

SU-E-I-22: Dependence On Calibration Phantom and Field Area of the Conversion Factor Used to Calculate Skin Dose During Neuro-Interventional Fluoroscopic Procedures

Purpose: To determine the appropriate calibration factor to use when calculating skin dose with our real-time dose-tracking system (DTS) during neuro-interventional fluoroscopic procedures by evaluating the difference in backscatter from different phantoms and as a function of entrance-skin field area. Methods: We developed a dose-tracking system to calculate and graphically display the cumulative skin-dose distribution in real time. To calibrate the DTS for neuro-interventional procedures, a phantom is needed that closely approximates the scattering properties of the head. We compared the x-ray backscatter from eight phantoms: 20-cm-thick solid water, 16-cm diameter water-filled container, 16-cm CTDI phantom, modified-ANSI head phantom, 20-cm-thick PMMA, Kyoto-Kagaku PBU- 50 head, Phantom-Labs SK-150 head, and RSD RS-240T head. The phantoms were placed on the patient table with the entrance surface at 15 cm tube-side from the isocenter of a Toshiba Infinix C-arm, and the entrance-skin exposure was measured with a calibrated 6-cc PTW ionization chamber. The measurement included primary radiation, backscatter from the phantom and forward scatter from the table and pad. The variation in entrance-skin exposure was also measured as a function of the skin-entrance area for a 30x30 cm by 20-cm-thick PMMA phantom and the SK-150 head phantom using four different added beam filters. Results: The entranceskin exposure values measured for eight different phantoms differed by up to 12%, while the ratio of entrance exposure of all phantoms relative to solid water showed less than 3% variation with kVp. The change in entrance-skin exposure with entrance-skin area was found to differ for the SK-150 head compared to the 20-cm PMMA phantom and the variation with field area was dependent on the added beam filtration. Conclusion: To accurately calculate skin dose for neuro-interventional procedures with the DTS, the phantom for calibration should be carefully chosen since different phantoms can contribute different backscatter for identical exposure parameters. Research supported in part by Toshiba Medical Systems and NIH Grants R43FD0158401, R44FD0158402 and R01EB002873.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-J-30: Benchmark Image-Based TCP Calculation for Evaluation of PTV Margins for Lung SBRT Patients

Purpose: Tumor control probability (TCP) calculated with accumulated radiation doses may help design appropriate treatment margins. Image registration errors, however, may compromise the calculated TCP. The purpose of this study is to develop benchmark CT images to quantify registration-induced errors in the accumulated doses and their corresponding TCP. Methods: 4DCT images were registered from end-inhale (EI) to end-exhale (EE) using a “demons” algorithm. The demons DVFs were corrected by an FEM model to get realistic deformation fields. The FEM DVFs were used to warp the EI images to create the FEM-simulated images. The two images combined with the FEM DVF formed a benchmark model. Maximum intensity projection (MIP) images, created from the EI and simulated images, were used to develop IMRT plans. Two plans with 3 and 5 mm margins were developed for each patient. With these plans, radiation doses were recalculated on the simulated images and warped back to the EI images using the FEM DVFs to get the accumulated doses. The Elastix software was used to register the FEM-simulated images to the EI images. TCPs calculated with the Elastix-accumulated doses were compared with those generated by the FEM to get the TCP error of the Elastix registrations. Results: For six lung patients, the mean Elastix registration error ranged from 0.93 to 1.98 mm. Their relative dose errors in PTV were between 0.28% and 6.8% for 3mm margin plans, and between 0.29% and 6.3% for 5mm-margin plans. As the PTV margin reduced from 5 to 3 mm, the mean TCP error of the Elastix-reconstructed doses increased from 2.0% to 2.9%, and the mean NTCP errors decreased from 1.2% to 1.1%. Conclusion: Patient-specific benchmark images can be used to evaluate the impact of registration errors on the computed TCPs, and may help select appropriate PTV margins for lung SBRT patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-101: Retroactive Calculation of TLD and Film Dose in Anthropomorphic Phantom as Assessment of Updated TPS Performance

Purpose: To demonstrate a quick and comprehensive method verifying the accuracy of the updated dose model by recalculating dose distribution in an anthropomorphic phantom with a new version of the TPS and comparing the results to measured values. Methods: CT images and IMRT plan of an RPC anthropomorphic head phantom, previously calculated by Pinnacle 9.0, was re-computed using Pinnacle 9.2 and 9.6. The dosimeters within the phantom include four TLD capsules representing a primary PTV, two TLD capsules representing a secondary PTV, and two TLD capsules representing an organ at risk. Also included were three sheets of Gafchromic film. Performance of the updated TPS version was assessed by recalculating point doses and dose profiles corresponding to TLD and film position respectively and then comparing the results to reported values by the RPC. Results: Comparing calculated doses to reported measured doses from the RPC yielded an average disagreement of 1.48%, 2.04% and 2.10% for versions 9.0, 9.2, 9.6 respectively. Computed doses points all meet the RPC's passing criteria with the exception of the point representing the superior organ at risk in version 9.6. However, qualitative analysis of the recalculated dose profiles showed improved agreement with those of the RPC, especially in the penumbra region. Conclusion: This work has demonstrated the calculation results of Pinnacle 9.2 and 9.6 vs 9.0 version. Additionally, this study illustrates a method for the user to gain confidence upgrade to a newer version of the treatment planning system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-175: Proton Dose Calculation On Scatter-Corrected CBCT Image: Feasibility Study for Adaptive Proton Therapy

Purpose: To demonstrate feasibility of proton dose calculation on scattercorrected CBCT images for the purpose of adaptive proton therapy. Methods: Two CBCT image sets were acquired from a prostate cancer patient and a thorax phantom using an on-board imaging system of an Elekta infinity linear accelerator. 2-D scatter maps were estimated using a previously introduced CT-based technique, and were subtracted from each raw projection image. A CBCT image set was then reconstructed with an open source reconstruction toolkit (RTK). Conversion from the CBCT number to HU was performed by soft tissue-based shifting with reference to the plan CT. Passively scattered proton plans were simulated on the plan CT and corrected/uncorrected CBCT images using the XiO treatment planning system. For quantitative evaluation, water equivalent path length (WEPL) was compared in those treatment plans. Results: The scatter correction method significantly improved image quality and HU accuracy in the prostate case where large scatter artifacts were obvious. However, the correction technique showed limited effects on the thorax case that was associated with fewer scatter artifacts. Mean absolute WEPL errors from the plans with the uncorrected and corrected images were 1.3 mm and 5.1 mm in the thorax case and 13.5 mm and 3.1 mm in the prostate case. The prostate plan dose distribution of the corrected image demonstrated better agreement with the reference one than that of the uncorrected image. Conclusion: A priori CT-based CBCT scatter correction can reduce the proton dose calculation error when large scatter artifacts are involved. If scatter artifacts are low, an uncorrected CBCT image is also promising for proton dose calculation when it is calibrated with the soft-tissue based shifting.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-275: Impact of the Intra and Inter Observer Variability in the Delineation of Parotid Glands On the Dose Calculation During Head and Neck Helical Tomotherapy

Purpose: The intra- and inter-observer variability in delineation of the parotids on the kilo-voltage computed tomography (kVCT) and mega-voltage computed tomography (MVCT) were examined to establish their impact on the dose calculation during adaptive head and neck helical tomotherapy (HT). Methods: Three observers delineated left and right parotids for ten randomly selected patients with oropharynx cancer treated on HT. The pre-treatment kVCT and the MVCT from the first fraction of irradiation were selected to delineation. The delineation procedure was repeated three times by each observer. The parotids were delineated according to the institutional protocol. The analyses included intra-observer reproducibility and inter-structure, -observer and -modality variability of the volume and dose. Results: The differences between the left and right parotid outlines were not statistically significant (p>0.3). The reproducibility of the delineation was confirmed for each observer on the kVCT (p>0.2) and on the MVCT (p>0.1). The inter-observer variability of the outlines was significant (p<0.001) as well as the inter-modality variability (p<0.006). The parotids delineated on the MVCT were 10% smaller than on the kVCT. The inter-observer variability of the parotids delineation did not affect the average dose (p=0.096 on the kVCT and p=0.176 on the MVCT). The dose calculated on the MVCT was higher by 3.3% than dose from the kVCT (p=0.009). Conclusion: Usage of the institutional protocols for the parotids delineation reduces intra-observer variability and increases reproducibility of the outlines. These protocols do not eliminate delineation differences between the observers, but these differences are not clinically significant and do not affect average doses in the parotids. The volumes of the parotids delineated on the MVCT are smaller than on the kVCT, which affects the differences in the calculated doses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-53: Benchmarking a Monte Carlo Model for Patient Plane Leakage Calculations of Low Energy 6MV Unique Linacs

Purpose: To validate the physics parameters of a Monte Carlo model for patient plane leakage calculations on the 6MV Unique linac by comparing the simulations against IEC patient plane leakage measurements. The benchmarked model can further be used for shielding design optimization, to predict leakage in the proximity of intended treatment fields, reduce the system weight and cost, and improve components reliability. Methods: The treatment head geometry of the Unique linac was simulated in Geant4 (v9.4.p02 with “Opt3” standard electromagnetic physics list) based on CAD drawings of all collimation and shielding components projected from the target to the area within 2m from isocenter. A 4×4m2 scorer was inserted 1m from the target in the patient plane and multiple phase space files were recorded by performing a 40-node computing cluster simulation on the EC2 cloud. The photon energy fluence was calculated relative to the value at isocenter for a 10×10cm2 field using 10×10mm2 bins. Tungsten blocks were parked accordingly to represent MLC120. The secondary particle contamination to patient plane was eliminated by “killing” those particles prior to the primary collimator entrance using a “kill-plane”, which represented the upper head shielding components not being modeled. Both IEC patient-plane leakage and X/Y-jaws transmission were simulated. Results: The contribution of photons to energy fluence was 0.064% on average, in excellent agreement with the experimental data available at 0.5, 1.0, and 1.5m from isocenter, characterized by an average leakage of 0.045% and a maximum leakage of 0.085%. X- and Y-jaws transmissions of 0.43% and 0.44% were found in good agreement with measurements of 0.48% and 0.43%, respectively. Conclusion: A Geant4 model based on energy fluence calculations for the 6MV Unique linac was created and validated using IEC patient plane leakage measurements. The “kill-plane” has effectively eliminated electron contamination to patient plane in these simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

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↗
DOE OSTI2014

SU-E-T-162: Evaluation of Dose Calculation of RayStation Planning System in Heterogeneous Media

Purpose: To investigate the clinical reliability of heterogeneity-based dose algorithm using RayStation treatment planning system v.4.0. Methods: The collapsed cone dose calculations in RayStation (RaySearch, Sweden) were compared with the measurements (ion chamber and EBT2 film) and with an in-house Monte Carlo algorithm. A heterogeneous multi-layer phantom and CT images of 4 lung cancer patients were used here. The phantom, composed of multiple solid water slabs and Styrofoams, was irradiated with 6MV beams perpendicular to the layers. The MLC-defined field sizes were 5×5, 10×10, 15×15 and 20×20cm{sup 2}. The chamber was positioned at center of central solid water layer, and the films were placed at interfaces of solid water and Styrofoam. The RayStation dose and Monte Carlo dose were compared by performing absolute gamma analysis (3mm/3%): 1D gamma for PDD in the phantom and 3D gamma for patient volumes receiving dose above 10% of maximum dose. Results: The point dose differences between RayStation and ion chamber measurement were smaller than 1% for all of the field sizes. Between RayStation and film measurement, 5×5cm2 field had the maximum differences : <4mm for the penumbra and <0.3mm for the field width at all Styrofoam-and-solid-water interfaces. The absolute gamma analysis showed good agreement between RayStation and Monte Carlo. For PDD along beam axis in the phantom, the 1D gamma was 95.4, 98.6, 99.6 and 99.3% for field size 5×5, 10×10, 15×15 and 20×202 respectively. For dose comparison using patient CT images, 3D gamma was > 95% for all the patients. Conclusion: With respect to ion chamber/film measurement and Monte Carlo calculation, the collapsed cone algorithm in RayStation computed reasonable dose in both phantom and patient cases. Heterogeneity-based dose calculation of RayStation is clinically acceptable in heterogeneous media.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-191: First Principle Calculation of Quantum Yield in Photodynamic Therapy

Purpose: We present a first-principle method to calculate the spin transfer efficiency in oxygen induced by any photon fields especially in MeV energy range. The optical pumping is mediated through photosensitizers, e.g., porphyrin and/or ensemble of quantum dots. Methods: Under normal conditions, oxygen molecules are in the relatively non-reactive triplet state. In the presence of certain photosensitizer compounds such as porphyrins, electromagnetic radiation of specific wavelengths can excite oxygen to highly reactive singlet state. With selective uptake of photosensitizers by certain malignant cells, photon irradiation of phosensitized tumors can lead to selective killing of cancer cells. This is the basis of photodynamic therapy (PDT). Despite several attempts, PDT has not been clinically successful except in limited superficial cancers. Many parameters such as photon energy, conjugation with quantum dots etc. can be potentially combined with PDT in order to extend the role of PDT in cancer management. The key quantity for this optimization is the spin transfer efficiency in oxygen by any photon field. The first principle calculation model presented here, is an attempt to fill this need. We employ stochastic density matrix description of the quantum jumps and the rate equation methods in quantum optics based on Markov/Poisson processes and calculate time evolution of the population of the optically pumped singlet oxygen. Results: The results demonstrate the feasibility of our model in showing the dependence of the optical yield in generating spin-singlet oxygen on the experimental conditions. The adjustable variables can be tuned to maximize the population of the singlet oxygen hence the efficacy of the photodynamic therapy. Conclusion: The present model can be employed to fit and analyze the experimental data and possibly to assist researchers in optimizing the experimental conditions in photodynamic therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-330: To Analyze the Calculation Error of Live Dose-Volume Indices Applying 4D-CT in Radiotherapy for PTVs Within the Liver Completely

Purpose: To study the variation rule of normal liver dose-volume indices calculation for the liver malignancy patients whose plan target volumes were in the liver completely in all breath phases. Methods: Ten patients who accepted radiotherapy for malignant tumor were selected in our study. All patients underwent 4D-CT simulation and 3D-CT simulation in free breathing(FB). 4D-CT was sorted into 10 series CT images according to breath phase, named CT0, CT10 to CT90, respectively. And GTVs were contoured on different CT series, and the individual target volume(ITV) was obtained by merging 10 GTVs from 4D-CT. The PTVs were obtained from ITV applying margins. The PTVs were not beyond the boundary of liver in all breath phase observed by dynamic 4D-CT. The radiotherapy plans were designed and irradiation dose was calculated on 3D-CT images, and the livers were contoured on different series CT images and mapped to 3D-CT images applying rigid registration. To compare the dose-volume difference of livers based on distinct CT images. Results: (1)The liver volumes were similar on 4D-CT and 3D-CT images(CTFB 1485±500cm{sup 3}, CT0 1413±377cm{sup 3}, CT10 1409±396cm{sup 3}, CT20 1419±418cm{sup 3},CT30 1485±500cm{sup 3}, CT40 1438±392cm{sup 3}, CT50 1437±404cm{sup 3}, CT60 1439±409cm{sup 3}, CT70 1408±393cm{sup 3}, CT80 1384±397cm{sup 3}, CT90 1398±397cm{sup 3}; F=0.064,p=1.00) (2) The PTVs volume were 30.17±14.62cm{sup 3};(3) The mean dose and V5 to V10 of liver were similar among 4D-CT different series CT images(p>0.05), and the indices varied less than ±4% refer to liver on CT50. Conclusion: The calculation affection of liver dose-volume indices induced by breath motion were not significant for the PTV within liver completely as estimation before. And more objective prediction indices for radiation induced l.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-416: VMAT Dose Calculations Using Cone Beam CT Images: A Preliminary Study

Purpose: Cone beam CT (CBCT) images have been used routinely for patient positioning throughout the treatment course. However, use of CBCT for dose calculation is still investigational. The purpose of this study is to assess the utility of CBCT images for Volumetric Modulated Arc Therapy (VMAT) plan dose calculation. Methods: A CATPHAN 504 phantom (The Phantom Laboratory, Salem, NY) was used to compare the dosimetric and geometric accuracy between conventional CT and CBCT (in both full and half fan modes). Hounsfield units (HU) profiles at different density areas were evaluated. A C shape target that surrounds a central avoidance structure was created and a VMAT plan was generated on the CT images and copied to the CBCT phantom images. Patient studies included three brain patients, and one head and neck (H'N) patient. VMAT plans generated on the patients treatment planning CT was applied to CBCT images obtained during the first treatment. Isodose distributions and dosevolume- histograms (DVHs) were compared. Results: For the phantom study, the HU difference between CT and CBCT is within 100 (maximum 96 HU for Teflon CBCT images in full fan mode). The impact of these differences on the calculated dose distributions was clinically insignificant. In both phantom and patient studies, target DVHs based on CBCT images were in excellent agreement with those based on planning CT images. Mean, Median, near minimum (D98%), and near maximum (D2%) doses agreed within 0-2.5%. A slightly larger discrepancy is observed in the patient studies compared to that seen in the phantom study, (0-1% vs. 0 - 2.5%). Conclusion: CBCT images can be used to accurately predict dosimetric results, without any HU correction. It is feasible to use CBCT to evaluate the actual dose delivered at each fraction. The dosimetric consequences resulting from tumor response and patient geometry changes could be monitored.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-439: An Improved Formula of Scatter-To-Primary Ratio for Photon Dose Calculation

Purpose: Scatter-to-primary ratio (SPR) is an important dosimetric quantity that describes the contribution from the scatter photons in an external photon beam. The purpose of this study is to develop an improved analytical formula to describe SPR as a function of circular field size (r) and depth (d) using Monte Carlo (MC) simulation. Methods: MC simulation was performed for Mohan photon spectra (Co-60, 4, 6, 10, 15, 23 MV) using EGSNRC code. Point-spread scatter dose kernels in water are generated. The scatter-to-primary ratio (SPR) is also calculated using MC simulation as a function of field size for circular field size with radius r and depth d. The doses from forward scatter and backscatter photons are calculated using a convolution of the point-spread scatter dose kernel and by accounting for scatter photons contributing to dose before (z' d) reaching the depth of interest, d, where z' is the location of scatter photons, respectively. The depth dependence of the ratio of the forward scatter and backscatter doses is determined as a function of depth and field size. Results: We are able to improve the existing 3-parameter (a, w, d0) empirical formula for SPR by introducing depth dependence for one of the parameter d0, which becomes 0 for deeper depths. The depth dependence of d0 can be directly calculated as a ratio of backscatter-to-forward scatter doses for otherwise the same field and depth. With the improved empirical formula, we can fit SPR for all megavoltage photon beams to within 2%. Existing 3-parameter formula cannot fit SPR data for Co-60 to better than 3.1%. Conclusion: An improved empirical formula is developed to fit SPR for all megavoltage photon energies to within 2%.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-491: Importance of Energy Dependent Protons Per MU Calibration Factors in IMPT Dose Calculations Using Monte Carlo Technique

Purpose: In intensity modulated proton therapy (IMPT), energy dependent, protons per monitor unit (MU) calibration factors are important parameters that determine absolute dose values from energy deposition data obtained from Monte Carlo (MC) simulations. Purpose of this study was to assess the sensitivity of MC-computed absolute dose distributions to the protons/MU calibration factors in IMPT. Methods: A “verification plan” (i.e., treatment beams applied individually to water phantom) of a head and neck patient plan was calculated using MC technique. The patient plan had three beams; one posterior-anterior (PA); two anterior oblique. Dose prescription was 66 Gy in 30 fractions. Of the total MUs, 58% was delivered in PA beam, 25% and 17% in other two. Energy deposition data obtained from the MC simulation were converted to Gy using energy dependent protons/MU calibrations factors obtained from two methods. First method is based on experimental measurements and MC simulations. Second is based on hand calculations, based on how many ion pairs were produced per proton in the dose monitor and how many ion pairs is equal to 1 MU (vendor recommended method). Dose distributions obtained from method one was compared with those from method two. Results: Average difference of 8% in protons/MU calibration factors between method one and two converted into 27 % difference in absolute dose values for PA beam; although dose distributions preserved the shape of 3D dose distribution qualitatively, they were different quantitatively. For two oblique beams, significant difference in absolute dose was not observed. Conclusion: Results demonstrate that protons/MU calibration factors can have a significant impact on absolute dose values in IMPT depending on the fraction of MUs delivered. When number of MUs increases the effect due to the calibration factors amplify. In determining protons/MU calibration factors, experimental method should be preferred in MC dose calculations. Research supported by National Cancer Institute grant P01CA021239.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-F-19A-10: Recalculation and Reporting Clinical HDR 192-Ir Head and Neck Dose Distributions Using Model Based Dose Calculation

Purpose: To retrospectively re-calculate dose distributions for selected head and neck cancer patients, earlier treated with HDR 192Ir brachytherapy, using Monte Carlo (MC) simulations and compare results to distributions from the planning system derived using TG43 formalism. To study differences between dose to medium (as obtained with the MC code) and dose to water in medium as obtained through (1) ratios of stopping powers and (2) ratios of mass energy absorption coefficients between water and medium. Methods: The MC code Algebra was used to calculate dose distributions according to earlier actual treatment plans using anonymized plan data and CT images in DICOM format. Ratios of stopping power and mass energy absorption coefficients for water with various media obtained from 192-Ir spectra were used in toggling between dose to water and dose to media. Results: Differences between initial planned TG43 dose distributions and the doses to media calculated by MC are insignificant in the target volume. Differences are moderate (within 4–5 % at distances of 3–4 cm) but increase with distance and are most notable in bone and at the patient surface. Differences between dose to water and dose to medium are within 1-2% when using mass energy absorption coefficients to toggle between the two quantities but increase to above 10% for bone using stopping power ratios. Conclusion: MC predicts target doses for head and neck cancer patients in close agreement with TG43. MC yields improved dose estimations outside the target where a larger fraction of dose is from scattered photons. It is important with awareness and a clear reporting of absorbed dose values in using model based algorithms. Differences in bone media can exceed 10% depending on how dose to water in medium is defined.

60 APPLIED LIFE SCIENCES↗