Results for “Software dependability” — Page 34

DOE OSTI2021

Validation of Maui PSCAD Model: Motivation, Methodology, and Lessons Learned: Preprint

As power systems become increasingly dependent on inverter-based resources, the dominant physical and control interactions that determine system stability may no longer be accurately captured in positive sequence dynamic simulation tools; electromagnetic transient (EMT) simulations may be needed. This paper describes the process of creating and validating an EMT model of the entire Maui island transmission system in PSCAD software. The model is validated against the utility's PSSE model and against field data for an event that consisted of a single-phase fault followed by a generation trip. The methodology and lessons learned are summarized here.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2015

High-Performance Design Patterns for Modern Fortran

This paper presents ideas for using coordinate-free numerics in modern Fortran to achieve code flexibility in the partial differential equation (PDE) domain. We also show how Fortran, over the last few decades, has changed to become a language well-suited for state-of-the-art software development. Fortran’s new coarray distributed data structure, the language’s class mechanism, and its side-effect-free, pure procedure capability provide the scaffolding on which we implement HPC software. These features empower compilers to organize parallel computations with efficient communication. We present some programming patterns that support asynchronous evaluation of expressions comprised of parallel operations on distributed data. We implemented these patterns using coarrays and the message passing interface (MPI). We compared the codes’ complexity and performance. The MPI code is much more complex and depends on external libraries. The MPI code on Cray hardware using the Cray compiler is 1.5–2 times faster than the coarray code on the same hardware. The Intel compiler implements coarrays atop Intel’s MPI library with the result apparently being 2–2.5 times slower than manually coded MPI despite exhibiting nearly linear scaling efficiency. As compilers mature and further improvements to coarrays comes in Fortran 2015, we expect this performance gap to narrow.

Haveraaen, Magne↗

DOE OSTI2020

Calculation of Moisture-Dependent, Anisotropic Parameters Supporting the Hanford Site's Composite Analysis, Cumulative Impact Evaluation, and Performance Assessments

Variable, moisture (tension) dependent anisotropy (MDA) in unsaturated soils is an effective, large-scale (macroscopic) flow property that results from media heterogeneities at a smaller scale. The purpose of this environmental calculation file (ECF) is to document the variable MDA calculations using a script routine. The accuracy of the script calculations are cross-checked against results using other software. Appendices for this ECF are Appendix A, which documents supporting script use; Appendix B, which provides supporting spreadsheet calculations; and Appendix C, which provides supporting utility calculations. The work is in support of vadose zone flow and transport models used in performance assessment (PA), composite analysis (CA), and cumulative impact evaluation (CIE) at the Hanford Site.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Characterization of extreme precipitation within atmospheric river events over California

Atmospheric rivers (ARs) are large, spatially coherent weather systems with high concentrations of elevated water vapor. These systems often cause severe downpours and flooding over the western coastal United States – and with the availability of more atmospheric moisture in the future under global warming we expect ARs to play an important role as potential causes of extreme precipitation changes. Therefore, we aim to investigate changes in extreme precipitation properties correlated with AR events in a warmer climate, which are large-scale meteorological patterns affecting the weather and climate of California. We have recently developed the TECA (Toolkit for Extreme Climate Analysis) software for automatically identifying and tracking features in climate data sets. Specifically, we can now identify ARs that make landfall on the western coast of North America. Based on this detection procedure, we can investigate the impact of ARs by exploring the spatial extent of AR precipitation using climate model (CMIP5) simulations and characterize spatial patterns of dependence for future projections between AR precipitation extremes under climate change within the statistical framework. Our results show that AR events in the future RCP (Representative Concentration Pathway)8.5 scenario (2076–2100) tend to produce heavier rainfall with higher frequency and longer days than events from the historical run (1981–2005). We also find that the dependence between extreme precipitation events has a shorter spatial range, within localized areas in California, under the high future emissions scenario than under the historical run.

Jeon, S.↗

DOE OSTI2019

Phase field modeling of recrystallization phenomena

In this work, we studied the recrystallization phenomena in copper using phase field modeling technique. EBSD information is successfully translated on to MOOSE and quantitative phase field simulations for Copper are performed on MOOSE. DREAM 3D software is utilized to transform EBSD data into a format which is then used to conduct phase field simulations, with initial grain structure obtained from EBSD data. Several simulations are performed and a comparison is made against various cases: isotropic grain growth without any deformed grains, deformed grains with constant dislocation density, and deformed grains with orientation dependent dislocation density.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Effect of occupant behavior on peak cooling and dehumidification loads in typical and high-efficiency homes

Residential building codes and voluntary labeling programs are continually increasing the energy efficiency requirements of residential buildings. Improving a building's thermal enclosure, installing the ductwork in conditioned space, and improving the building's airtightness results in significant reductions in externally-driven sensible and latent cooling loads. As a building's efficiency is improved, occupant-related internal gains become a larger portion of the building sensible and latent loads. Additionally, internal gains are highly uncertain compared to other load components. In this study, we use a stochastic approach to simulate occupant-related internal gains and compare the internal gains to other sensible and latent heat sources in four house efficiency levels in 10 U.S. climates using whole-building energy simulation software. We compare the expected range in occupant-related internal gains to other building characteristics such as cooling set point, air infiltration rate, and mechanical ventilation rate. We show that in high-efficiency homes, sensible internal gains vary from less than 10% to greater than 40% of the building sensible load under peak total cooling conditions depending on climate and internal gain profile. Likewise, latent internal gains vary from less than 10% to more than 60% of the building latent load under peak total cooling and peak dehumidification conditions depending on climate and internal gain profile.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Generalized framework for likelihood-based field-level inference of growth rate from velocity and density fields

Measuring the growth rate of large-scale structures ( f ) as a function of redshift has the potential to break degeneracies between modified gravity and dark energy models, when combined with expansion-rate probes. Direct estimates of peculiar velocities of galaxies have attracted interest as a means of estimating fσ 8 . In particular, field-level methods can be used to fit the field nuisance parameter along with cosmological parameters simultaneously. This article aims to provide the community with a unified framework for the theoretical modeling of the likelihood-based field-level inference by performing fast field covariance calculations for velocity and density fields. Our purpose is to lay the foundations for a nonlinear extension of the likelihood-based method at the field level. We have developed a generalized framework, implemented in the dedicated software flip to perform a likelihood-based inference of fσ 8 . We derived a new field covariance model, which includes wide-angle corrections. We also included the models previously described in the literature inside our framework. We compared their performance against ours, and we validated our model by comparing it with the two-point statistics of a recent N-body simulation. The tests we performed have allowed us to validate our software and determine the appropriate wavenumber range to integrate our covariance model and its validity in terms of separation. Our framework allows for a wider wavenumber coverage to be used in our calculations than in previous works, which is particularly interesting for nonlinear model extensions. Finally, our generalized framework allows us to efficiently perform a survey geometry-dependent Fisher forecast of the fσ 8 parameter. We show that the Fisher forecast method we developed gives an error bar that is 30% closer to a full likelihood-based estimation than a standard volume Fisher forecast.

Ravoux, Corentin↗

DOE OSTI2024

Reconstructing the reflectivity of liquid surfaces from grazing incidence X-ray off-specular scattering data

The capillary wave model of a liquid surface predicts both the X-ray specular reflection and the diffuse scattering around it. A quantitative method is presented to obtain the X-ray reflectivity (XRR) from a liquid surface through the diffuse scattering data around the specular reflection measured using a grazing incidence X-ray off-specular scattering (GIXOS) geometry at a fixed horizontal offset angle with respect to the plane of incidence. With this approach the entire Q z -dependent reflectivity profile can be obtained at a single, fixed incident angle. This permits a much faster acquisition of the profile than with conventional reflectometry, where the incident angle must be scanned point by point to obtain a Q z -dependent profile. The XRR derived from the GIXOS-measured diffuse scattering, referred to in this paper as pseudo-reflectivity, provides a larger Q z range compared with the reflectivity measured by conventional reflectometry. Transforming the GIXOS-measured diffuse scattering profile to pseudo-XRR opens up the GIXOS method to widely available specular XRR analysis software tools. Here the GIXOS-derived pseudo-XRR is compared with the XRR measured by specular reflectometry from two simple vapor–liquid interfaces at different surface tension, and from a hexadecyltrimethylammonium bromide monolayer on a water surface. For the simple liquids, excellent agreement (beyond 11 orders of magnitude in signal) is found between the two methods, supporting the approach of using GIXOS-measured diffuse scattering to derive reflectivities. Pseudo-XRR obtained at different horizontal offset angles with respect to the plane of incidence yields indistinguishable results, and this supports the robustness of the GIXOS-XRR approach. The pseudo-XRR method can be extended to soft thin films on a liquid surface, and criteria are established for the applicability of the approach.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Multi-Objective Optimization for Size and Resilience of Spiking Neural Networks

Inspired by the connectivity mechanisms in the brain, neuromorphic computing architectures model Spiking Neural Networks (SNNs) in silicon. As such, neuromorphic architectures are designed and developed with the goal of having small, low power chips that can perform control and machine learning tasks. However, the power consumption of the developed hardware can greatly depend on the size of the network that is being evaluated on the chip. Furthermore, the accuracy of a trained SNN that is evaluated on chip can change due to voltage and current variations in the hardware that perturb the learned weights of the network. While efforts are made on the hardware side to minimize those perturbations, a software based strategy to make the deployed networks more resilient can help further alleviate that issue. In this work, we study Spiking Neural Networks in two neuromorphic architecture implementations with the goal of decreasing their size, while at the same time increasing their resiliency to hardware faults. We leverage an evolutionary algorithm to train the SNNs and propose a multiobjective fitness function to optimize the size and resiliency of the SNN. We demonstrate that this strategy leads to well-performing, small-sized networks that are more resilient to hardware faults.

Dimovska, Mihaela↗

DOE OSTI2026

braggedgemodeling

braggedgemodeling (bem) is an open-source Python package for modeling neutron Bragg-edge imaging. It computes the wavelength-dependent total neutron cross-section of a material (coherent and incoherent elastic, coherent and incoherent inelastic scattering, and absorption) from its crystal structure, and implements the March-Dollase texture model and the Jorgensen peak profile, supporting quantitative analysis of energy-resolved neutron imaging data (phase, stress/strain, and texture). Published in the Journal of Open Source Software (2018).

Lin, Jiao [Oak Ridge National Laboratory (ORNL), O↗

DOE OSTI2022

Efficient Implementation of Artificial Neural Networks for Sensor Data Analysis Based on a Genetic Algorithm

The reliability of many industrial processes depends on the sensor system. However, these sensors can be affected by noise, perturbations and failures. Hence, sensor monitoring and diagnosis are fundamental to guarantee the quality of an industrial process. Nowadays, artificial neural networks (ANN) are widely used in sensor signal processing and diagnosis. However, those ANNs usually require many artificial neurons, being difficult to implement in software and hardware due to their high computational costs. This paper presents an optimized implementation of artificial neurons in ANNs for sensor data analysis using a Genetic Algorithm (GA). The objective of GA is to find an adequate segmentation to reduce the activation function approximation error. One of the advantages of the proposed approach is that the cost function used in GA considers the effect of factors such as the ANN architecture or the number of bits used in arithmetic operations. The proposed ANN implementation technique aims to get the best possible approximation for a specific ANN architecture, making easier its implementation in software and hardware. Simulation and experimental results using FPGA (Field Programmable Gate Array) prove the advantages of the proposed approach for implementing sensor data analysis systems based on ANNs.

D estefani, André↗

DOE OSTI2015

Detailed investigation of optoelectronic and microstructural properties of plasma polymerized cyclohexane thin films: Dependence on the radiofrequency power

Optical properties of polymerized cyclohexane films deposited by radiofrequency plasma enhanced chemical vapor deposition technique at different radiofrequency powers onto glass and silicon substrates, are studied and correlated with the microstructure of the films, using a combination of atomic force microscopy, Raman and Fourier Transformer Infrared spectroscopy and optical measurements. The optical constants such as refractive index n, dielectric permittivity ε and extinction k and absorption α coefficients, are extracted from transmission and reflection spectra through the commercial software CODE. These constants lead, by using common theoretical models as Cauchy, Lorentz, Tauc and single effective oscillator, to the determination of the static refractive index n{sub s} and permittivity ε{sub s}, the plasma frequency ω{sub p}, the carrier density to effective mass ratio N/m{sub e}{sup *}, the optical conductivity σ{sub oc}, the optical band gap E{sub g} and the oscillation and dispersion energies E{sub 0} and E{sub d}, respectively. We find that n, ε{sub s}, ω{sub p}, N/m{sub e}{sup *}, E{sub d}, increase with radiofrequency power, while E{sub g} and E{sub 0} decrease in the same range of power. These results are well correlated with those obtained from atomic force microscopy, Raman and infrared measurements. They also indicate that the increase of the radiofrequency power promotes the fragmentation of the precursor and increases the carbon C-sp{sup 2} hybridization proportion, which results in an improvement of the optoelectronic properties of the films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

A comprehensive method of estimating electric fields from vector magnetic field and Doppler measurements

Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic models of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal-toroidal decomposition (PTD) of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the 'PTD-Doppler-FLCT Ideal' (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the FISHPACK software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (ANMHD) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than 1% error in the total Poynting flux and a 10% error in the helicity flux rate at a normal viewing angle (θ = 0) and less than 25% and 10% errors, respectively, at large viewing angles (θ < 60°). We compare our results with other inversion methods at zero viewing angle and find that our method's estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Numerical evaluation of the bispectrum in multiple field inflation—the transport approach with code

We present a complete framework for numerical calculation of the power spectrum and bispectrum in canonical inflation with an arbitrary number of light or heavy fields. Our method includes all relevant effects at tree-level in the loop expansion, including (i) interference between growing and decaying modes near horizon exit; (ii) correlation and coupling between species near horizon exit and on superhorizon scales; (iii) contributions from mass terms; and (iv) all contributions from coupling to gravity. We track the evolution of each correlation function from the vacuum state through horizon exit and the superhorizon regime, with no need to match quantum and classical parts of the calculation; when integrated, our approach corresponds exactly with the tree-level Schwinger or 'in-in' formulation of quantum field theory. In this paper we give the equations necessary to evolve all two- and three-point correlation functions together with suitable initial conditions. The final formalism is suitable to compute the amplitude, shape, and scale dependence of the bispectrum in models with | f {sub NL}| of order unity or less, which are a target for future galaxy surveys such as Euclid, DESI and LSST. As an illustration we apply our framework to a number of examples, obtaining quantitatively accurate predictions for their bispectra for the first time. Two accompanying reports describe publicly-available software packages that implement the method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

SU-E-T-75: Commissioning Optically Stimulated Luminescence Dosimeters for Fast Neutron Therapy

Purpose: Fast neutrons therapy used at the University of Washington is clinically proven to be more effective than photon therapy in treating salivary gland and other cancers. A nanodot optically stimulated luminescence (OSL) system was chosen to be commissioned for patient in vivo dosimetry for neutron therapy. The OSL-based radiation detectors are not susceptible to radiation damage caused by neutrons compared to diodes or MOSFET systems. Methods: An In-Light microStar OSL system was commissioned for in vivo use by radiating Landauer nanodots with neutrons generated from 50.0 MeV protons accelerated onto a beryllium target. The OSLs were calibrated the depth of maximum dose in solid water localized to 150 cm SAD isocenter in a 10.3 cm square field. Linearity was tested over a typical clinical dose fractionation range i.e. 0 to 150 neutron-cGy. Correction factors for transient signal fading, trap depletion, gantry angle, field size, and wedge factor dependencies were also evaluated. The OSLs were photo-bleached between radiations using a tungsten-halogen lamp. Results: Landauer sensitivity factors published for each nanodot are valid for measuring photon and electron doses but do not apply for neutron irradiation. Individually calculated nanodot calibration factors exhibited a 2–5% improvement over calibration factors computed by the microStar InLight software. Transient fading effects had a significant impact on neutron dose reading accuracy compared to photon and electron in vivo dosimetry. Greater accuracy can be achieved by calibrating and reading each dosimeter within 1–2 hours after irradiation. No additional OSL correction factors were needed for field size, gantry angle, or wedge factors in solid water phantom measurements. Conclusion: OSL detectors are a useful for neutron beam in vivo dosimetry verification. Dosimetric accuracy comparable to conventional diode systems can be achieved. Accounting for transient fading effects during the neutron beam calibration is a critical component for achieving comparable accuracy.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2014

MO-F-16A-05: Spot-Scanning Proton Therapy Patient-Specific Quality Assurance: Our Methodology and Results From 295 Treatment Plans

Purpose: To report on the patient-specific quality assurance (PSQA) results for 295 spot-scanning proton therapy treatment plans from the MD Anderson PTC-Houston. We show how the results differed by treatment site and how they were affected by the treatment plan optimization method and by a range shifter in the treatment field. We also discuss some causes of PSQA problems. Methods: The PSQA procedure, which is designed to verify both the accuracy of the treatment planning system's (Eclipse™ v8.9) dose calculations and the dose delivery of the Hitachi PROBEAT synchrotron, consists of (1) an end-to-end test in which the beam is delivered and measured at the prescribed gantry angle, and (2) additional dose plane measurements made from gantry angle 270°. HPlusQA™ software automatically performs the gamma analysis with criteria 3% (dose tolerance), 3 mm (distance-to-agreement, DTA) and 2%, 2 mm. Passing is defined as at least 90% of the pixels having a gamma score less than 1. Results: The PSQA gamma passing rate was 96.2% for 3%, 3 mm, and 85.3% for 2%, 2 mm. The rate depended on the treatment site. For example, the 3%, 3 mm passing rate was 95% for head and neck plans, vs 100% for prostate plans. The passing rates of multi- vs. single-field optimization plans did not significantly differ. However, the rate for fields with range shifters was 94.8±0.6%, vs 99.0±0.6% for those without (p = 0.002). Longitudinal dose gradients caused most of the low scores. Overestimation of the calculated dose proximal to the spread-out Bragg peak (SOBP) caused many of the others. Conclusion: The planned and delivered doses consistently agreed within tolerance levels. Minor dose modeling deficiencies remain proximal to the SOBP. The 3% dose tolerance, 3 mm DTA, with 90% pixel passing rate is a reasonable action level for 2D gamma comparisons.

60 APPLIED LIFE SCIENCES↗