Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static Code Analysis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Methodological approach to the sensitivity analysis of failure effects in modern digital I and C systems - 136

Modern I and C systems of nuclear power plants increasingly use digital equipment, software-based applications for I and C functions and different types and topologies of communication networks as well. The structure, operation and communication of digital safety-related I and C systems are strongly influenced by the implemented automatic fault detection and fault treatment procedures as barriers against failure propagation. Thus, these I and C systems can change their operating states automatically during the operation time, e.g. due to the detection and correction of certain failure types. Reasonably because of these dynamic properties of digital I and C systems the use of classical fault tree and event tree methods solely for the reliability analysis of modern I and C systems is less suitable. These methods basically apply static models and do not consider the temporal changes of the operating state of the I and C system or equipment appropriately. In this paper, a recently developed approach for the sensitivity analysis of digital I and C systems will be introduced that allows to investigate the dynamic behavior of such systems in cases of internal faults and to quantify the effects on the system reliability. The development and testing of the analysis methodology is carried out using models of generic I and C systems with different typical architectures of stepwise increased complexity. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Develop and verify soil/structure interaction for pile/foundation interaction

Phase II of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was used to verify the implementation of a new soil-structure interaction (SSI) model for use within offshore wind turbine modeling software. The REDWIN Macro-element model implemented and verified in this study enables a computationally efficient way to model the linear and nonlinear SSI problem, including hysteretic damping, of a monopile structure. The modeling approach was integrated into several modeling tools and a series of increasingly complex simulations was conducted using the IEA 10MW reference turbine mounted on a monopile support structure to verify the coupling between the tools and the REDWIN Macro-element SSI model. This campaign includes only numerical verification between various software and modeling approaches so no experimental measurements are available. The load cases (LC) considered include: LC1 – static response of the tower and substructure LC2 – frequency and mode-shape analysis of the tower and substructure LC3 – response of the tower and substructure due to wind-only loading LC4 – response of the tower and substructure due to wave-only loading LC5 – response of the tower and substructure due to wind and wave loading. Detailed properties of the modeled system are found in the following reference, “Bergua, Roger, Amy Robertson, Jason Jonkman, and Andy Platt. 2021. "Specification Document for OC6 Phase II: Verification of an Advanced Soil-Structure Interaction Model for Offshore Wind Turbines.” Golden, CO: National Renewable Energy Laboratory. NREL/TP-5000-79938. https://www.nlr.gov/docs/fy21osti/79938.pdf. Details on the results from the OC6 Phase II project can be found in the following reference, “Bergua R, Robertson A, Jonkman J, et al. OC6 Phase II: Integration and verification of a new soil–structure interaction model for offshore wind design.” Wind Energy. 2022;25(5):793-810. doi:10.1002/we.2698

17 WIND ENERGY↗
DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

CMB and BAO constraints for an induced gravity dark energy model with a quartic potential

We study the predictions for structure formation in an induced gravity dark energy model with a quartic potential. By developing a dedicated Einstein-Boltzmann code, we study self-consistently the dynamics of homogeneous cosmology and of linear perturbations without using any parametrization. By evolving linear perturbations with initial conditions in the radiation era, we accurately recover the quasi-static analytic approximation in the matter dominated era. We use PLANCK 2013 data and a compilation of baryonic acoustic oscillation (BAO) data to constrain the coupling γ to the Ricci curvature and the other cosmological parameters. By connecting the gravitational constant in the Einstein equation to the one measured in a Cavendish-like experiment, we find γ < 0.0012 at 95% CL with PLANCK 2013 and BAO data. This is the tightest cosmological constraint on γ and on the corresponding derived post-Newtonian parameters. Because of a degeneracy between γ and the Hubble constant H{sub 0}, we show how larger values for γ are allowed, but not preferred at a significant statistical level, when local measurements of H{sub 0} are combined in the analysis with PLANCK 2013 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

RISK MONITORING CAPABILITES FROM DYNAMIC PRA DATA

Risk Monitors (RMs) are software tools that are employed in nuclear power plants to monitor the safety properties of the plant itself during operation. These safety properties are generated using classical PRA methods which are typically based on static Boolean structures such as Fault-Trees (FTs) and Event-Trees (ETs). The outcomes of these methods are Core Damage Frequency (CDF) and risk importance values of plant components (e.g., valves and pumps). Examples of Risk Importance Measures (RIMs) commonly used are: Risk Achievement Worth (RAW), Risk Reduction Worth (RRW), Birnbaum (B) and Fussell-Vesely (FV) [1,2]. In a RM these values are updated in real time once plant components status is updated (e.g., under maintenance, off-line, testing). In addition, RMs can be used to optimize maintenance schedules so those safety features of the plant are below plant standard values. This paper presents an approach to create temporal profiles of RIMs (i.e., RM capabilities) given a temporal change of the plant configuration. Instead of employing data generated by ETs and FTs, the proposed approach will employ data generated by Dynamic PRA methods. In contrast to classical PRA methods, Dynamic PRA methods couple stochastic methods with safety analysis codes to determine risk associate to complex systems such as nuclear plants. Compared to classical PRA methods, they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. We will describe how classical RIMs can be generated by Dynamic PRA data and then we will show how the temporal profiles of the RIMs can be generated given a temporal profile of the system configuration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Investigation of static and dynamic behavior of functionally graded piezoelectric actuated Poly-Si micro cantilever probe

In the present paper a novel functionally graded piezoelectric (FGP) actuated Poly-Si micro cantilever probe is proposed for atomic force microscope. The shear piezoelectric coefficient d{sub 15} has much higher value than coupling coefficients d{sub 31} and d{sub 33}, hence in the present work the micro cantilever beam actuated by d{sub 15} effect is utilized. The material properties are graded in the thickness direction of actuator by a simple power law. A three dimensional finite element analysis has been performed using COMSOL Multiphysics® (version 4.2) software. Tip deflection and free vibration analysis for the micro cantilever probe has been done. The results presented in the paper shall be useful in the design of micro cantilever probe and their subsequent utilization in atomic force microscopes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Ab Initio Molecular Dynamics Analysis of Xenon Incorporation in Uranium Dioxide

The fission process in nuclear fuel produces volatile fission product gases, such as xenon, that migrate out of the fuel and into surrounding structures. Consequently, these fission product gases may affect thermal performance and integrity of the fuel. Fission gas release models, used in modern fuel performance codes such as BISON, are currently under investigation to improve their accuracy and ability to simulate the fission gas release phenomenon and its consequences. For uranium dioxide (UO{sub 2}), the most commonly used fuel in commercial power reactors, fission gas release models are known to have significant uncertainties. In the interest of improving such models, several complex diffusional mechanisms have been examined using a variety of techniques such as Nudge Elastic Band (NEB), utilizing both Density Functional Theory (DFT) and semi-empirical inter-atomic potential functions. One proposed diffusion mechanism in UO{sub 2} is through tetra-vacancy (i.e., 2U{sub V} + 2O{sub V}) structures, which is expected to provide a suitable means for the incorporation and diffusion of xenon in the fuel crystal. In this work, a model was developed to study the incorporation of xenon in UO{sub 2} via a tetra-vacancy using ab initio molecular dynamics (AIMD). This is a first principles technique that abandons the use empirical force fields as usually is the case of classical MD analysis. Furthermore, in contrast with static transition analysis tools such as NEB, AIMD is able to capture the dynamics of a given process. The translation of xenon within a tetra-vacancy has been investigated using ab-initio molecular dynamics. The transition frequency, barrier energy, and incorporation energy were extracted from the model, demonstrating its ability to capture a variety of parameters and its potential for use in studying fission gas diffusive behavior. Future refinements to the model will seek to develop additional insight into Xe and Kr mobility and will likely involve executing more elaborate AIMD simulations to investigate this phenomenon. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

SU-F-T-526: A Comparative Study On Gating Efficiency of Varian RPM Device and Calypso System

Purpose: In general, the linear accelerator is gated using respiratory signal obtained by way of external sensors to account for the breathing motion during radiotherapy. One of the commonly used gating devices is the Varian RPM device. Calypso system that uses electromagnetic tracking of implanted or surface transponders could also be used for gating. The aim of this study is to compare the gating efficiency of RPM device and the calypso system by phantom studies. Methods: An ArcCheck insert was used as the phantom with a Gafchromic film placed in its holder. The ArcCheck insert was placed on a Motion Sim platform and moved in the longitudinal direction simulating a respiratory motion with a period of 5 seconds and amplitude of ±6mm. The Gafchromic film was exposed to a 2 × 2cm{sup 2} field, i) with the phantom static, ii) phantom moving but ungated iii) gated with gating window of 2mm and 3mm. This was repeated with Calypso system using surface transponders with the same gating window. The Gafchromic films were read with an EPSON 11000 flatbed scanner and analysed with ‘Medphysto’ software. Results: The full width at half maximum (FWHM) as measured with film at the level of the film holder was 1.65cm when the phantom was static. FWHM measured with phantom moving and without gating was 1.16 cm and penumbra was 7 mm (80–20%) on both sides. When the beam was gated with 2 mm gating window the FWHM was 1.8 cm with RPM device and 1.9 cm with Calypso. Similarly, when the beam was gated with 3 mm window, the FWHM was 1.9cm with RPM device and 2cm with Calypso. Conclusion: This work suggests that the gating efficiency of RPM device is better than that of the Calypso with surface transponder, with reference to the latency in gating.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Study of spatial spin-modulated structures by Mössbauer spectroscopy using SpectrRelax

SpectrRelax is an application for analysis and fitting of absorption and emission Mössbauer spectra. It includes a large selection of static and relaxation spectrum models, and allows fitting and searching for optimal model parameters. Recently, we have added new models for Mössbauer spectra of nuclides in spatial spin modulated structures. In these structures, spin density or direction changes in a periodic way along a single direction, and this wave is incommensurate with the underlying lattice. The models include Spin/Charge density wave, where the shape of this wave is represented as a sum of odd harmonics, Anharmonic spin modulation where the spin direction has a cycloidal type modulation, and a Spiral-like spin structure, in which magnetic moments rotate in a plane perpendicular to the wave propagation vector, forming a spiral.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

SU-E-T-83: A Study On Evaluating the Directional Dependency of 2D Seven 29 Ion Chamber Array Clinically with Different IMRT Plans

Purpose: To evaluate the directional dependency of 2D seven 29 ion chamber array clinically with different IMRT plans. Methods: 25 patients already treated with IMRT plans were selected for the study. Verification plans were created for each treatment plan in eclipse 10 treatment planning system using the AAA algorithm with the 2D array and the Octavius CT phantom. Verification plans were done 2 times for a single patient. First plan with real IMRT (plan-related approach) and second plan with zero degree gantry angle (field-related approach). Measurements were performed on a Varian Clinac-iX, linear accelerator equipped with a millennium 120 multileaf collimator. Fluence was measured for all the delivered plans and analyzed using the verisoft software. Comparison was done by selecting the fluence delivered in static gantry (zero degree gantry) versus IMRT with real gantry angles. Results: The gamma pass percentage is greater than 97 % for all IMRT delivered with zero gantry angle and between 95%–98% for real gantry angles. Dose difference between the TPS calculated and measured for IMRT delivered with zero gantry angle was found to be between (0.03 to 0.06Gy) and with real gantry angles between (0.02 to 0.05Gy). There is a significant difference between the gamma analysis between the zero degree and true angle with a significance of 0.002. Standard deviation of gamma pass percentage between the IMRT plans with zero gantry angle was 0.68 and for IMRT with true gantry angle was found to be 0.74. Conclusion: The gamma analysis for IMRT with zero degree gantry angles shows higher pass percentage than IMRT delivered with true gantry angles. Verification plans delivered with true gantry angles lower the verification accuracy when 2D array is used for measurement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors

Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

SU-C-BRD-06: Sensitivity Study of An Automated System to Acquire and Analyze EPID Exit Dose Images

Purpose: The dosimetric consequences of errors in patient setup or beam delivery and anatomical changes are not readily known. A new product, PerFRACTION (Sun Nuclear Corporation), is designed to detect these errors by comparing the EPID exit dose image from each field of each fraction to those from baseline fraction images. This work investigates the sensitivity of PerFRACTION to detect the deviation of induced errors in a variety of realistic scenarios. Methods: Eight plans were created mimicking potential delivery or setup errors. The plans consisted of a nominal field and the field with an induced error. These were used to irradiate the EPID simulating multiple fractions with and without the error. Integrated EPID images were acquired in clinical mode and saved in ARIA. PerFRACTION automatically pulls the images into its database and performs the user defined comparison. In some cases, images were manually pushed to PerFRACTION. We varied the distance-to-agreement or dose tolerance until PerFRACTION showed failing pixels in the affected region and recorded the values. We induced errors of 1mm and greater in jaw, MLC, and couch position, 2 degree collimation rotation (patient yaw), and 0.5% to 1.5% in machine output. Both static and arc fields with the rails in or out were also acquired and compared. Results: PerFRACTION detected position errors of the jaws, MLC, and couch with an accuracy of better than 0.5 mm, and 0.2 degrees for collimator and gantry error. PerFRACTION detected a machine output error within 0.2% and detected the change in rail position. Conclusion: A new automated system for monitoring daily treatments for machine or patient variations from the first fraction using integrated EPID images was found to be sensitive enough to detect small positional, angular, and dosimetric errors within 0.5mm, 0.2 degrees, and 0.2%, respectively. Sun Nuclear Corporation has provided a software license for the product described.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Nevada National Security Site 2018 Data Report: Groundwater Monitoring Program, Area 5 Radioactive Waste Management Site

This report is a compilation of the groundwater and waste cell leachate sampling results from the Area 5 Radioactive Waste Management Site (RWMS) at the Nevada National Security Site, Nye County, Nevada. Groundwater samples from the aquifer immediately below the Area 5 RWMS have been collected and analyzed and static water levels have been measured in this aquifer since 1993. These groundwater data are evaluated for evidence of effects on the aquifer related to the Area 5 RWMS. Leachate from one lined mixed waste cell at the Area 5 RWMS (Cell 18) has been sampled since the cell opened in 2011. These leachate data are analyzed for hazardous contaminants to determine appropriate leachate handling and disposal. A second lined mixed waste cell at the Area 5 RWMS (Cell 25) began receiving waste in August 2018, but no leachate samples were collected or analyzed during 2018. This report includes data from 2014 through 2018 but focuses on the 2018 data. During 2018, groundwater samples were collected and static water levels were measured at three wells surrounding the Area 5 RWMS. Groundwater samples were collected at wells UE5PW-1, UE5PW-2, and UE5PW-3 on March 6 and August 14, 2018. Static water levels were measured at each of these wells on March 5, June 4, August 13, and October 25, 2018. Groundwater monitoring at the Area 5 RWMS complies with Title 40 Code of Federal Regulations (CFR) Part 264.97, General ground-water monitoring requirements , and 40 CFR 264.98, Detection monitoring program. Groundwater samples were analyzed for indicators of contamination (pH, specific conductance [SC], total organic carbon, total organic halides, and tritium) and, beginning in 2017, toxicity characteristic metals (arsenic, barium, cadmium, chromium, lead, mercury, selenium, and silver). General water chemistry parameters (calcium, iron, magnesium, manganese, potassium, sodium, sulfate, chloride, fluoride, bicarbonate, and silicate) were also measured. Results from samples collected in 2018 are within the limits for each parameter established by the Nevada Division of Environmental Protection in Resource Conservation and Recovery Act (RCRA) Permit NEV HW0101, Revision 6. Data from the shallow aquifer indicate there has been no measurable impact to the uppermost aquifer from the Area 5 RWMS, and there were no significant changes in measured groundwater parameters compared to previous years. Leachate from lined mixed waste Cell 18 drains into a sump and is pumped into an above ground tank. Samples are collected from the tank when the leachate volume approaches the 3,000-gallon tank capacity. During 2018, leachate samples were collected on March 1, 2018 and August 1, 2018. All leachate analysis results are below the regulatory levels for toxicity characteristic contaminants (40 CFR 261.24) and below the investigation levels for polychlorinated biphenyls (PCBs), tritium, pH, and SC specified in RCRA Permit NEV HW0101, Revision 6. No quantifiable PCB levels were detected in any leachate sample. Based on the leachate analysis results, RCRA Permit NEV HW0101 allowed all leachate collected from Cell 18 to be pumped from the collection tank and used for dust control at Cell 18.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Measuring Impact: Evaluating Thermal Zoning Simplification on Energy Efficiency Measures Analysis

Building Energy Modeling (BEM) is a versatile tool for designing, retrofitting, ensuring code compliance, obtaining certifications, qualifying for incentives, and enabling real-time building control. However, capturing all the details of building geometry for thermal zoning can be time-consuming, costly, and sometimes computationally challenging. As a result, modelers have been applying zoning simplification based on factors such as space functions and internal loads, as well as relying on their experience and judgment while adhering to zoning rules outlined in industry standards. Despite the prevalence of this common practice, a notable gap exists in the literature regarding studies quantifying the influence of simplified thermal zoning on the evaluation of Energy Efficiency Measures (EEMs). Recognizing this gap, this paper seeks to contribute to the field by enhancing the understanding of how the simplification of thermal zoning influences the evaluation of EEMs against a baseline design. The study utilized a medium office prototype model with a detailed floor plan featuring over 20 zones per floor covering diverse functional spaces with varying internal loads and occupancy schedules. A standard thermal zoning strategy outlined in ASHRAE Standard 90.1 Appendix G was employed as the simplified zoning method. This strategy condenses the zoning into a core zone and four perimeter zones per floor. It was compared with the detailed zoning approach, which involves one zone per space. Common Energy EEMs, such as enhanced envelope, high-efficiency appliances and equipment, and HVAC controls, were individually implemented and evaluated. The results indicate that the performance comparison between the two zoning methods varies depending on the type of measures considered. Basic measures, such as adding wall insulation, demonstrate similar energy impacts, while advanced HVAC control measures, such as static pressure reset, exhibit a more substantial difference that cannot be overlooked.

Xie, Jiarong↗
DOE OSTI2014

SU-E-T-03: 3D GPU-Accelerated Secondary Checks of Radiation Therapy Treatment Plans

Purpose: Redundant treatment verifications in conformal and intensity-modulated radiation therapy techniques are traditionally performed with single point calculations. New solutions can replace these checks with 3D treatment plan verifications. This work describes a software tool (Mobius3D, Mobius Medical Systems) that uses a GPU-accelerated collapsed cone algorithm to perform 3D independent verifications of TPS calculations. Methods: Mobius3D comes with reference beam models for common linear accelerators. The system uses an independently developed collapsed cone algorithm updated with recent enhancements. 144 isotropically-spaced cones are used for each voxel for calculations. These complex calculations can be sped up by using GPUs. Mobius3D calculate dose using DICOM information coming from TPS (CT, RT Struct, RT Plan RT Dose). DVH-metrics and 3D gamma tests can be used to compare both TPS and secondary calculations. 170 patients treated with all common techniques as 3DCFRT (including wedged), static and dynamic IMRT and VMAT have been successfully verified with this solution. Results: Calculation times are between 3–5 minutes for 3DCFRT treatments and 15–20 for most complex dMLC and VMAT plans. For all PTVs mean dose and 90% coverage differences are (1.12±0.97)% and (0.68±1.19)%, respectively. Mean dose discrepancies for all OARs is (0.64±1.00)%. 3D gamma (global, 3%/3 mm) analysis shows a mean passing rate of (97.8 ± 3.0)% for PTVs and (99.0±3.0)% for OARs. 3D gamma pasing rate for all voxels in CT has a mean value of (98.5±1.6)%. Conclusion: Mobius3D is a powerful tool to verify all modalities of radiation therapy treatments. Dose discrepancies calculated by this system are in good agreement with TPS. The use of reference beam data results in time savings and can be used to avoid the propagation of errors in original beam data into our QA system. GPU calculations permit enhanced collapsed cone calculations with reasonable calculation times.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Dosimetric quality, accuracy, and deliverability of modulated radiotherapy treatments for spinal metastases

Cancer often metastasizes to the vertebra, and such metastases can be treated successfully using simple, static posterior or opposed-pair radiation fields. However, in some cases, including when re-irradiation is required, spinal cord avoidance becomes necessary and more complex treatment plans must be used. This study evaluated 16 sample intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT) treatment plans designed to treat 6 typical vertebral and paraspinal volumes using a standard prescription, with the aim of investigating the advantages and limitations of these treatment techniques and providing recommendations for their optimal use in vertebral treatments. Treatment plan quality and beam complexity metrics were evaluated using the Treatment And Dose Assessor (TADA) code. A portal-imaging–based quality assurance (QA) system was used to evaluate treatment delivery accuracy, and radiochromic film measurements were used to provide high-resolution verification of treatment plan dose accuracy, especially in the steep dose gradient regions between each vertebral target and spinal cord. All treatment modalities delivered approximately the same doses and the same levels of dose heterogeneity to each planning target volume (PTV), although the minimum PTV doses in the vertebral plans were substantially lower than the prescription, because of the requirement that the plans meet a strict constraint on the dose to the spinal cord and cord planning risk volume (PRV). All plans met required dose constraints on all organs at risk, and all measured PTV-cord dose gradients were steeper than planned. Beam complexity analysis suggested that the IMRT treatment plans were more deliverable (less complex, leading to greater QA success) than the VMAT treatment plans, although the IMRT plans also took more time to deliver. The accuracy and deliverability of VMAT treatment plans were found to be substantially increased by limiting the number of monitor units (MU) per beam at the optimization stage, and thereby limiting beam modulation complexity. The VMAT arcs that were optimized with MU limitation had higher QA pass rates as well as higher modulation complexity scores (less complexity), lower modulation indices (less modulation), lower MU per beam, larger beam segments, and fewer small apertures than the VMAT arcs that were optimized without MU limitation. It is recommended that VMAT treatments for vertebral volumes, where the PTV abuts or surrounds the spinal cord, should be optimized with MU limitation. IMRT treatments may be preferable to the VMAT treatments, for dosimetry and deliverability reasons, but may be inappropriate for some patients because of their increased treatment delivery time.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Plasma Surface Interactions: Bridging from the Surface to the Micron Frontier through Leadership Class Computing (Final Technical Progress Report)

The performance demands on plasma-facing components (PFCs) in future fusion power plants are beyond the capability of current materials. In defining the plasma-surface interaction (PSI), it is evident that three coupled spatial regions influence PFC materials evolution and performance. These regions consist of (1) the edge and scrape-off layer region of the plasma, (2) the near-surface material response to extreme thermal and particle fluxes under the influence of, and feedback to, the plasma sheath, and (3) the structural material’s response to an intense, 14 MeV peaked neutron spectrum, which produces very high concentrations of transmuted elements through (n,p) and (n,α) reactions and structural material property degradation. The objective of this project is to develop and deploy validated, high-performance simulation tools capable of predicting the performance of tungsten-based plasma-facing components (PFCs) in a burning fusion plasma environment, which includes modeling surface morphology evolution in either erosion or re-deposition regimes, and the recycling of hydrogen species. This requires the development of a leadership-scale computational code to predict PFC behavior that is well integrated to a suite of multiscale modeling techniques to bridge the scales needed to address complex physical and computational issues at the plasma surface interface and the transition region below the surface where neutron damage processes in the bulk dominate material behavior. Successful completion of this project will provide PFC simulation tools to evaluate steady-state performance of tungsten-based PFC and divertor components in burning plasma environments. This will enable the identification of critical experiments to confirm whether practical PFC solutions exist for magnetic fusion energy beyond ITER. More specifically, the research activities within this project focus on two broadly defined research thrusts: 1) Developing a new simulation code, Xolotl-PSI to predict PFC operating lifetime and performance, that is specifically designed to take advantage of leadership class computing facilities; and 2) Integrating and applying discrete particle-based, as well as continuum-based, multiscale modeling techniques to provide scientific discovery of the mechanisms controlling PFC and bulk materials evolution under fusion plasma and 14-MeV neutron exposure; in this thrust, the challenge of scale-bridging between atomistic/microstructural modeling and the continuum-scale of PFC device performance is achieved through multiscale modeling techniques. Throughout the grant period, research at UMass Amherst focused on the atomic-scale analysis of the dynamics and kinetics (transport and reactions) of small mobile helium clusters near surfaces and grain boundaries in plasma-exposed tungsten, as well as the proper parameterization and incorporation of these processes into our SciDAC team’s continuum-level cluster dynamics code (Xolotl). Our work also included characterization of large-scale molecular-dynamics (MD) simulation results of helium implantation into tungsten based on the findings of our analysis and comparisons of Xolotl simulation predictions with the large-scale MD simulation results for Xolotl validation purposes; emphasis was placed on the impact of helium cluster dynamics/kinetics on tungsten surface morphology, near-surface tungsten structure, and helium retention. Furthermore, our work included initiation of a research program on evaluating the impact of structural defects in tungsten due to plasma exposure (such as helium nanobubbles) on tungsten’s thermophysical properties, focusing on lattice thermal conductivity. Our studies contributed to a fundamental understanding of the effects of plasma-surface interactions on the dynamical response of plasma-facing materials in nuclear fusion devices and the development of predictive hierarchical multiscale computational tools for the quantitative description of the material’s dynamical response to plasma exposure.

36 MATERIALS SCIENCE↗