Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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 973 records · Page 54

Analysis of Composition B Sphere Impact Fragmentation Experiments

Since energetic material fragmentation governs the surface area available for burning and potential escalation to violent reaction modes (e.g., deflagration-to-­detonation transition and XDT modes), it is critical that hazard response models accurately capture impact damage response. The objective of this work is to analyze recent Composition B sphere impact experiments conducted by Naval Air Warfare Center Weapons Division (NAWCWD) that probe fragmentation response for impact speeds from 27 to 127 m/s. In those experiments, fragment size distribution was measured by dry sieving damaged samples. A power law relationship was developed for the normalized cumulative fragment number as a function of the normalized fragment mass, (N/ N$_{tot}$ = A (M$_{bin}$/ M$_{avg}$ ) n , across all impact speeds. The power law fits demonstrated good agreement with the experimental data. We also correlated the power law fitting parameters to the impact speed. The power law prefactor, A, obeyed a polynomial relationship while the power law exponent, n, obeyed a linear relationship in impact speed. While this preliminary analysis shows good agreement with Composition B fragmentation data, it was based on experiments over a limited range of impact speeds and a single sample size. Future modeling-­based studies are critical to improve our understanding of fragmentation response and to more confidently extrapolate to larger length-­scales and higher impact speeds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Area of Interest 1, CO 2 at the Interface: Nature and Dynamics of the Reservoir/Caprock Contact and Implications for Carbon Storage Performance

We examined the influence of geologic features present at the reservoir/caprock interface on the transmission of supercritical CO 2 into and through caprock. We focused on the case of deformation-band faults in reservoir lithologies that intersect the interface and transition to opening-mode fractures in caprock lithologies. Deformation-band faults are exceeding common in potential CO 2 injection units and our fieldwork in Utah indicates that this sort of transition is common. To quantify the impact of these interface features on flow and transport we first described the sedimentology and permeability characteristics of selected sites along the Navajo Sandstone (reservoir lithology) and Carmel Formation (caprock lithology) interface, and along the Slickrock Member (reservoir lithology) and Earthy Member (caprock lithology) of the Entrada Sandstone interface, and used this information to construct conceptual permeability models for numerical analysis. We then examined the impact of these structures on flow using single-phase and multiphase numerical flow models for these study sites. Key findings include: (1) Deformation-band faults strongly compartmentalize the reservoir and largely block cross-fault flow of supercritical CO 2 . (2) Significant flow of CO 2 through the fractures is possible, however, the magnitude is dependent on the small-scale geometry of the contact between the opening-mode fracture and the deformation band fault. (3) Due to the presence of permeable units in the caprock, caprock units are capable of storing significant volumes of CO 2 , particularly when the fracture network does not extend all the way through the caprock. The large-scale distribution of these deformation-bandfault-to-opening-mode-fractures is related to the curvature of the beds, with greater densities of fractures in high curvature regions. We also examined core and outcrops from the Mount Simon Sandstone and Eau Claire Formation reservoir/caprock interface in order to extend our work to a reservoir/caprock pair this is currently being assessed for long-term carbon storage. These analyses indicate that interface features similar to those observed at the Utah sites 3 were not observed. Although not directly related to our main study topic, one byproduct of our investigation is documentation of exceptionally high degrees of heterogeneity in the pore-size distribution of the Mount Simon Sandstone. This suggests that the unit has a greater-than-normal potential for residual trapping of supercritical CO 2 .

54 ENVIRONMENTAL SCIENCES↗

α-Imaging Confirmed Efficient Targeting of CD₄₅-Positive Cells After ²¹¹At-Radioimmunotherapy for Hematopoietic Cell Transplantation

Alpha-radioimmunotherapy (α-RIT) targeting CD45 may substitute for total body irradiation in hematopoietic cell transplantation (HCT) preparative regimens for lymphoma. Our goal was to optimize the anti-CD45 monoclonal antibody (MAb; CA12.10C12) protein dose for astatine-²¹¹(²¹¹At)-RIT, extending the analysis to include intra-organ ²¹¹At activity distribution and α-imaging-based small-scale dosimetry, along with imunohistochemical staining. Methods: Eight normal dogs were injected with either 0.75 (n=5) or 1.00 mg/kg (n=3) of ²¹¹At-B10-CA12.10C12 (11.5–27.6 MBq/kg). Two were euthanized and necropsied 19–22 hours postinjection (p.i.), and six received autologous HCT three days after ²¹¹At-RIT, following lymph node and bone marrow biopsies at 2–4 and/or 19 hours p.i. Blood was sampled to study toxicity and clearance; CD45 targeting was evaluated by flow cytometry. ²¹¹At localization and small scale dosimetry were assessed using two α-imaging : α-camera and iQID. Results: Uptake of ²¹¹At was highest in spleen (0.31–0.61 %IA/g), lymph nodes (0.02–0.16 %IA/g), liver (0.11–0.12 %IA/g), and marrow (0.06–0.08 %IA/g). Lymphocytes in blood and marrow were efficiently targeted using either MAb dose. Lymph nodes remained unsaturated, but displayed targeted ²¹¹At localization in T lymphocyte-rich areas. Absorbed doses to blood, marrow, and lymph nodes were estimated at 3.9, 3.0, and 4.2 Gy/210 MBq, respectively. All transplanted dogs experienced transient hepatic toxicity. Liver enzyme levels were temporarily elevated in 5 of 6 dogs; 1 treated with 1.00 mg MAb/kg developed ascites and was euthanized 136 days after HCT. Conclusion: ²¹¹At-anti-CD45 RIT with 0.75 mg MAb/kg efficiently targeted blood and marrow without severe toxicity. Dosimetry calculations and observed radiation-induced effects indicated that sufficient ²¹¹At-B10-CA12.10C12 localization was achieved for efficient conditioning for HCT.

Radioimmunotherapy, astatine-211, alpha imaging, C↗

THE COVERING FACTOR OF WARM DUST IN WEAK EMISSION-LINE ACTIVE GALACTIC NUCLEI

Weak emission-line active galactic nuclei (WLAGNs) are radio-quiet active galactic nuclei (AGNs) that have nearly featureless optical spectra. We investigate the ultraviolet to mid-infrared spectral energy distributions of 73 WLAGNs (0.4 < z < 3) and find that most of them are similar to normal AGNs. We also calculate the covering factor of warm dust of these 73 WLAGNs. No significant difference is indicated by a KS test between the covering factor of WLAGNs and normal AGNs in the common range of bolometric luminosity. The implication for several models of WLAGNs is discussed. The super-Eddington accretion is unlikely to be the dominant reason for the featureless spectrum of a WLAGN. The present results are still consistent with the evolution scenario, i.e., WLAGNs are in a special stage of AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

Diverse Chemo-Dynamical Properties of Nitrogen-Rich Stars Identified from Low-Resolution Spectra

The second generation of stars in the globular clusters (GCs) of the Milky Way (MW) exhibit unusually high N, Na, or Al, compared to typical Galactic halo stars at similar metallicities. The halo field stars enhanced with such elements are believed to have originated in disrupted GCs or escaped from existing GCs. We identify such stars in the metallicity range –3.0 < [Fe/H] < 0.0 from a sample of ~ 36,800 giant stars observed in the Sloan Digital Sky Survey and Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, and present their dynamical properties. The N-rich population (NRP) and N-normal population (NNP) among our giant sample do not exhibit similarities in either in their metallicity distribution function (MDF) or dynamical properties. We find that, even though the MDF of the NRP looks similar to that of the MW’s GCs in the range of [Fe/H] < –1.0, our analysis of the dynamical properties does not indicate similarities between them in the same metallicity range, implying that the escaped members from existing GCs may account for a small fraction of our N-rich stars, or the orbits of the present GCs have been altered by the dynamical friction of the MW. We also find a significant increase in the fraction of N-rich stars in the halo field in the very metal-poor (VMP; [Fe/H] < –2.0) regime, comprising up to ~ 20% of the fraction of the N-rich stars below [Fe/H] = –2.5, hinting that partially or fully destroyed VMP GCs may have in some degree contributed to the Galactic halo. A more detailed dynamical analysis of the NRP reveals that our sample of N-rich stars do not share a single common origin. Although a substantial fraction of the N-rich stars seem to originate from the GCs formed in situ, more than 60% of them are not associated with those of typical Galactic populations, but probably have extragalactic origins associated with Gaia Sausage/Enceladus, Sequoia, and Sagittarius dwarf galaxies, as well as with presently unrecognized progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗

Isolated nano-tendril bundles on tungsten surfaces exposed to radiofrequency helium plasma

The DIONISOS experiment is used to study the impact of RF helium (He) plasma on the surface morphology of tungsten (W). Modulation of incident He ion energy with a span on the order of 50 eV and below the sputtering threshold at a frequency of 13.56 MHz results in nano-tendril bundles (NTBs) and free-standing W whiskers on the surface at 1020 K. The NTBs are distributed intra-granularly with coverage of less than 10% while reaching up to 30 microns normal to the surface for He ion fluence of 5.8 x 10^25 m^-2 and a flux density of 1.0 x 10^22 m^-2 s^-1. Analysis of the NTB interior and sub-surface structure is provided through focused ion beam cross section.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of the 2017 Solar Eclipse on the Smart Grid

With the increasing interest in using solar energy as a major contributor to the use of renewable generation, and with the focus on using smart grids to optimize the use of electrical energy based on demand and resources from different locations, the need arises to know the moons position in the sky with respect to the sun. When a solar eclipse occurs, the moon disk might totally or partially shade the sun disk, which can affect the irradiance level from the sun disk, consequently affecting a resource on the electric grid. The moons position can then provide smart grid users with information about how potential total or partial solar eclipses might affect different locations on the grid so that other resources on the grid can be directed to where they might be needed when such phenomena occurs. At least five solar eclipses occur yearly at different locations on Earth, they can last 3 hours or more depending on the location, and they can affect smart grid users. On August 21, 2017, a partial and full solar eclipse occurred in many locations in the United States, including at the National Renewable Energy Laboratory in Golden, Colorado. Solar irradiance measurements during the eclipse were compared to the data generated by a model for validation at eight locations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Preprocessing Analysis of the Multi Isotope Process Monitor

The Multi-Isotope Process (MIP) Monitor is a nondestructive monitoring methodology to detect diversion and tampering of nuclear material at a reprocessing facility. The system measures the gamma-emitting isotopes in different stages of the reprocessing cycle through passive gamma spectroscopy. The indicators of isotopic distributions are measured and evaluated in near-real time by comparing measured spectra to expected, normal process spectra through multivariate analysis. The MIP monitor framework is compatible with medium resolution gamma detectors that can be easily deployed in an existing facility. In FY15, the MIP monitor was deployed at the H-Canyon Facility at Savannah River National Laboratory (SRNL) to develop and evaluate analysis algorithms for future in situ deployment. This paper discusses recent results from analysis of the data collected at SRNL. New results examine the preprocessing step of spectral alignment that is associated with the MIP monitor.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Instabilities of Jammed Packings of Frictionless Spheres Under Load

In this work, we consider the contribution to the density of vibrational states and the distribution of energy barrier heights of incipient instabilities in a glass modeled by a jammed packing of spheres. On approaching an instability, the frequency of a normal mode and the height of the energy barrier to cross into a new ground state both vanish. These instabilities produce a contribution to the density of vibrational states that scales as ω 3 at low frequencies ω, and a contribution to the distribution of energy barriers ΔH that scales as ΔH -1/3 at low barrier heights.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On a stochastic approach to a code performance estimation

The main goal of an efficient profiling of software is to minimize the runtime overhead under certain constraints and requirements. The traces built by a profiler during the work, affect the performance of the system itself. One of important aspect of an overhead arises from the randomness of variability in the context in which the application is embedded, e.g., due to possible cache misses, etc. Such uncertainty needs to be taken into account in the design phase. In order to overcome these difficulties we propose to investigate this issue through the analysis of the probability distribution of the difference between profiler’s times for the same code. The approximating model is based on the finite normal mixtures within the framework of the method of moving separation of mixtures. We demonstrate some results for the MATLAB profiler using plotting of 3D surfaces by the function surf. The idea can be used for an estimating of a program efficiency.

97 MATHEMATICS AND COMPUTING↗

Laboratory-developed contact models controlling instability on frictional faults: FORESHOCKS ON FRICTIONAL FAULTS

Laboratory experiments were performed on a polymethyl methacrylate (PMMA)-PMMA frictional interface in a direct shear apparatus in order to gain understanding of fault dynamics leading to gross rupture. Actual asperity sizes and locations along the interface were characterized using a pressure-sensitive film. Slow aseismic slip accumulated nonuniformly along the fault and showed dependency on the applied normal force-increased normal force resulted in higher slip gradients. The slow slip front propagated from the trailing (pushed) edge into a region of more densely distributed asperities at rates between 1 and 9.5 mm/s. Foreshocks were detected and displayed impulsive signals with source radii ranging between 0.21 and 1.09 mm; measurements made using the pressure-sensitive film were between 0.05 and 1.2 mm. The spatiotemporal clustering of foreshocks and their relation to the elastodynamic energy released was dependent on the normal force. In the region where foreshocks occurred, qualitative optical measurements of the asperities along the interface were used to visualize dynamic changes occurring during the slow slip phase. To better understand the nucleation process, a quasi-static asperity finite element (FE) model was developed and focused in the region where foreshocks clustered. The FE model consisted of 172 asperities, located and sized based on pressure-sensitive film measurements. The numerical model provides a plausible explanation as to why foreshocks cluster in space and observed a normal force dependency and lend credence to Ohnaka’s nucleation model.

58 GEOSCIENCES↗

Oister optical and near-infrared observations of type Iax supernova 2012Z

We report observations of the Type Iax supernova (SN Iax) 2012Z at optical and near-infrared (NIR) wavelengths from immediately after the explosion until ∼260 days after the maximum luminosity using the Optical and Infrared Synergetic Telescopes for Education and Research Target-of-Opportunity program and the Subaru Telescope. We found that the NIR light curve evolutions and color evolutions are similar to those of SNe Iax 2005hk and 2008ha. The NIR absolute magnitudes (M{sub J}∼−18.1 mag and M{sub H}∼−18.3 mag) and the rate of decline of the light curve (Δm{sub 15}(B) = 1.6±0.1 mag) are very similar to those of SN 2005hk (M{sub J}∼−17.7 mag, M{sub H}∼−18.0 mag, and Δm{sub 15}(B) ∼ 1.6 mag), yet differ significantly from SNe 2008ha and 2010ae (M{sub J} ∼ −14 to −15 mag and Δm{sub 15}(B) ∼ 2.4–2.7 mag). The estimated rise time is 12.0 ± 3.0 days, which is significantly shorter than that of SN 2005hk or any other SNe Ia. The rapid rise indicates that the {sup 56}Ni distribution may extend into the outer layer or that the effective opacity may be lower than that in normal SNe Ia. The late-phase spectrum exhibits broader emission lines than those of SN 2005hk by a factor of six to eight. Such high velocities of the emission lines indicate that the density profile of the inner ejecta extends more than that of SN 2005hk. We argue that the most favored explosion scenario is a “failed deflagration” model, although the pulsational delayed detonations is not excluded.

79 ASTRONOMY AND ASTROPHYSICS↗

Design and Performance Evaluation of Al{sub 2}O{sub 3}-SiC Composite for Direct-Bonded Copper Substrate

A computational material design approach is applied to propose a novel ceramic material for direct-bonded copper (DBC) substrate with enhanced thermal and structural performance. The material design inherently consists of many competing requirements that require careful decisions regarding key trade-offs in terms of material composition, inclusion size, shape, and distribution to achieve the target properties. The alumina-silicon (Al{sub 2}O{sub 3}-SiC) composite, as compared to commercial alumina, used in DBC is found to be the most suitable design among other candidates with improved thermal and structural properties. In order to study the performance characteristics and the effects of the new ceramic composite with improved properties in terms of structural behavior and fatigue life of the DBC substrate, the normal working and extreme thermal cycling conditions were simulated and analyzed using finite element method. The temperature, strain, and localized stress distribution within the substrate at a steady-state condition were analyzed, and the improved Coffin–Manson law was used to calculate the fatigue life of the substrate under extreme thermal cycling conditions. The proposed Al{sub 2}O{sub 3}-SiC composite is found to be more robust than the commercial alumina as DBC substrates considering the thermal–mechanical performance. The fatigue life cycle of the DBC substrate with the proposed material is predicted to be about two times longer than the commercial alumina DBC ceramic under transient thermal cycling test.

36 MATERIALS SCIENCE↗

Detection of Metastatic Breast and Thyroid Cancer in Lymph Nodes by Desorption Electrospray Ionization Mass Spectrometry Imaging

Ambient ionization mass spectrometry has been widely applied to image lipids and metabolites in primary cancer tissues with the purpose of detecting and understanding metabolic changes associated with cancer development and progression. Here, we report the use of desorption electrospray ionization mass spectrometry (DESI-MS) to image metastatic breast and thyroid cancer in human lymph node tissues. Our results show clear alterations in lipid and metabolite distributions detected in the mass spectra profiles from 42 samples of metastatic thyroid tumors, metastatic breast tumors, and normal lymph node tissues. 2D DESI-MS ion images of selected molecular species allowed discrimination and visualization of specific histologic features within tissue sections, including regions of metastatic cancer, adjacent normal lymph node, and fibrosis or adipose tissues, which strongly correlated with pathologic findings. In thyroid cancer metastasis, increased relative abundances of ceramides and glycerophosphoinisitols were observed. In breast cancer metastasis, increased relative abundances of various fatty acids and specific glycerophospholipids were seen. Trends in the alterations in fatty acyl chain composition of lipid species were also observed through detailed mass spectra evaluation and chemical identification of molecular species. The results obtained demonstrate DESI-MSI as a potential clinical tool for the detection of breast and thyroid cancer metastasis in lymph nodes, although further validation is needed. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

TU-D-209-03: Alignment of the Patient Graphic Model Using Fluoroscopic Images for Skin Dose Mapping

Purpose: The Dose Tracking System (DTS) was developed to provide realtime feedback of skin dose and dose rate during interventional fluoroscopic procedures. A color map on a 3D graphic of the patient represents the cumulative dose distribution on the skin. Automated image correlation algorithms are described which use the fluoroscopic procedure images to align and scale the patient graphic for more accurate dose mapping. Methods: Currently, the DTS employs manual patient graphic selection and alignment. To improve the accuracy of dose mapping and automate the software, various methods are explored to extract information about the beam location and patient morphology from the procedure images. To match patient anatomy with a reference projection image, preprocessing is first used, including edge enhancement, edge detection, and contour detection. Template matching algorithms from OpenCV are then employed to find the location of the beam. Once a match is found, the reference graphic is scaled and rotated to fit the patient, using image registration correlation functions in Matlab. The algorithm runs correlation functions for all points and maps all correlation confidences to a surface map. The highest point of correlation is used for alignment and scaling. The transformation data is saved for later model scaling. Results: Anatomic recognition is used to find matching features between model and image and image registration correlation provides for alignment and scaling at any rotation angle with less than onesecond runtime, and at noise levels in excess of 150% of those found in normal procedures. Conclusion: The algorithm provides the necessary scaling and alignment tools to improve the accuracy of dose distribution mapping on the patient graphic with the DTS. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗

Functional Data Analysis in NTCP Modeling: A New Method to Explore the Radiation Dose-Volume Effects

Purpose/Objective(s): To describe a novel method to explore radiation dose-volume effects. Functional data analysis is used to investigate the information contained in differential dose-volume histograms. The method is applied to the normal tissue complication probability modeling of rectal bleeding (RB) for patients irradiated in the prostatic bed by 3-dimensional conformal radiation therapy. Methods and Materials: Kernel density estimation was used to estimate the individual probability density functions from each of the 141 rectum differential dose-volume histograms. Functional principal component analysis was performed on the estimated probability density functions to explore the variation modes in the dose distribution. The functional principal components were then tested for association with RB using logistic regression adapted to functional covariates (FLR). For comparison, 3 other normal tissue complication probability models were considered: the Lyman-Kutcher-Burman model, logistic model based on standard dosimetric parameters (LM), and logistic model based on multivariate principal component analysis (PCA). Results: The incidence rate of grade ≥2 RB was 14%. V{sub 65Gy} was the most predictive factor for the LM (P=.058). The best fit for the Lyman-Kutcher-Burman model was obtained with n=0.12, m = 0.17, and TD50 = 72.6 Gy. In PCA and FLR, the components that describe the interdependence between the relative volumes exposed at intermediate and high doses were the most correlated to the complication. The FLR parameter function leads to a better understanding of the volume effect by including the treatment specificity in the delivered mechanistic information. For RB grade ≥2, patients with advanced age are significantly at risk (odds ratio, 1.123; 95% confidence interval, 1.03-1.22), and the fits of the LM, PCA, and functional principal component analysis models are significantly improved by including this clinical factor. Conclusion: Functional data analysis provides an attractive method for flexibly estimating the dose-volume effect for normal tissues in external radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Proposed linear energy transfer areal detector for protons using radiochromic film

Radiation therapy depends on predictably and reliably delivering dose to tumors and sparing normal tissues. Protons with kinetic energy of a few hundred MeV can selectively deposit dose to deep seated tumors without an exit dose, unlike x-rays. The better dose distribution is attributed to a phenomenon known as the Bragg peak. The Bragg peak is due to relatively high energy deposition within a given distance or high Linear Energy Transfer (LET). In addition, biological response to radiation depends on the dose, dose rate, and localized energy deposition patterns or LET. At present, the LET can only be measured at a given fixed point and the LET spatial distribution can only be inferred from calculations. The goal of this study is to develop and test a method to measure LET over extended areas. Traditionally, radiochromic films are used to measure dose distribution but not for LET distribution. We report the first use of these films for measuring the spatial distribution of the LET deposited by protons. The radiochromic film sensitivity diminishes for large LET. A mathematical model correlating the film sensitivity and LET is presented to justify relating LET and radiochromic film relative sensitivity. Protons were directed parallel to radiochromic film sandwiched between solid water slabs. This study proposes the scaled-normalized difference (SND) between the Treatment Planning system (TPS) and measured dose as the metric describing the LET. The SND is correlated with a Monte Carlo (MC) calculation of the LET spatial distribution for a large range of SNDs. A polynomial fit between the SND and MC LET is generated for protons having a single range of 20 cm with narrow Bragg peak. Coefficients from these fitted polynomial fits were applied to measured proton dose distributions with a variety of ranges. An identical procedure was applied to the protons deposited from Spread Out Bragg Peak and modulated by 5 cm. Gamma analysis is a method for comparing the calculated LET with the LET measured using radiochromic film at the pixel level over extended areas. Failure rates using gamma analysis are calculated for areas in the dose distribution using parameters of 25% of MC LET and 3 mm. The processed dose distributions find 5%–10% failure rates for the narrow 12.5 and 15 cm proton ranges and 10%–15% for proton ranges of 15, 17.5, and 20 cm and modulated by 5 cm. It is found through gamma analysis that the measured proton energy deposition in radiochromic film and TPS can be used to determine LET. This modified film dosimetry provides an experimental areal LET measurement that can verify MC calculations, support LET point measurements, possibly enhance biologically based proton treatment planning, and determine the polymerization process within the radiochromic film.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-T-135: Assessing the Clinical Impact of Approximations in Analytical Dose Calculations for Proton Therapy

Purpose: To estimate the clinical relevance of approximations made in analytical dose calculation methods (ADCs) used for treatment planning on tumor coverage and tumor control probability (TCP) in proton therapy. Methods: We compared dose distributions planned with ADC to delivered dose distributions (as determined by TOPAS Monte Carlo (MC) simulations). We investigated 10 patients per site for 5 treatment sites (head-and-neck, lung, breast, prostate, liver). We evaluated differences between the two dose distributions analyzing dosimetric indices based on the dose-volume-histograms, the γ-index and the TCP. The normal tissue complication probability (NTCP) was estimated for the bladder and anterior rectum for the prostate patients. Results: We find that the target doses are overestimated by the ADC by 1–2% on average for all patients considered. All dosimetric indices (the mean dose, D95, D50 and D02, the dose values covering 95%, 50% and 2% of the target volume, respectively) are predicted within 5% of the delivered dose. A γ-index with a 3%/3mm criteria had a passing rate for target volumes above 96% for all patients. The TCP predicted by the two algorithms was up to 2%, 2.5%, 6%, 6.5%, and 11% for liver and breast, prostate, head-and-neck and lung patients, respectively. Differences in NTCP for anterior-rectum and bladder for prostate patients were less than 3%. Conclusion: We show that ADC provide adequate dose distributions for most patients, however, they can Result in underdosage of the target by as much as 5%. The TCP was found to be up to 11% lower than predicted. Advanced dose-calculation methods like MC simulations may be necessary in proton therapy to ensure target coverage for heterogeneous patient geometries, in clinical trials comparing proton therapy to conventional radiotherapy to avoid biases due to systematic discrepancies in calculated dose distributions, and, if tighter range margins are considered. Fully funded by NIH grants.

60 APPLIED LIFE SCIENCES↗