Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-C-BRF-05: Design and Geometric Validation of An Externally and Internally Deformable, Programmable Lung Motion Phantom

Purpose: Most clinically-deployed strategies for respiratory motion management in lung radiotherapy (e.g., gating, tracking) use external markers that serve as surrogates for tumor motion. However, typical lung phantoms used to validate these strategies are rigid-exterior+rigid-interior or rigid-exterior+deformable-interior. Neither class adequately represents the human anatomy, which is deformable internally as well as externally. We describe the construction and experimental validation of a more realistic, externally- and internally-deformable, programmable lung phantom. Methods: The outer shell of a commercially-available lung phantom (RS- 1500, RSD Inc.) was used. The shell consists of a chest cavity with a flexible anterior surface, and embedded vertebrae, rib-cage and sternum. A 3-axis platform was programmed with sinusoidal and six patient-recorded lung tumor trajectories. The platform was used to drive a rigid foam ‘diaphragm’ that compressed/decompressed the phantom interior. Experimental characterization comprised of mapping the superior-inferior (SI) and anterior-posterior (AP) trajectories of external and internal radioopaque markers with kV x-ray fluoroscopy and correlating these with optical surface monitoring using the in-room VisionRT system. Results: The phantom correctly reproduced the programmed motion as well as realistic effects such as hysteresis. The reproducibility of marker trajectories over multiple runs for sinusoidal as well as patient traces, as characterized by fluoroscopy, was within 0.4 mm RMS error for internal as well as external markers. The motion trajectories of internal and external markers as measured by fluoroscopy were found to be highly correlated (R=0.97). Furthermore, motion trajectories of arbitrary points on the deforming phantom surface, as recorded by the VisionRT system also showed a high correlation with respect to the fluoroscopically-measured trajectories of internal markers (R=0.92). Conclusion: We have developed a realistic externally- and internally-deformable lung phantom that will serve as a valuable tool for clinical QA and motion management research. This work was supported through funding from the NIH and VisionRT Ltd. Amit Sawant has research funding from Varian Medical Systems, VisionRT and Elekta.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

The use of RapidArc volumetric-modulated arc therapy to deliver stereotactic radiosurgery and stereotactic body radiotherapy to intracranial and extracranial targets

Twenty-three targets in 16 patients treated with stereotactic radiosurgery (SRS) or stereotactic body radiotherapy (SBRT) were analyzed in terms of dosimetric homogeneity, target conformity, organ-at-risk (OAR) sparing, monitor unit (MU) usage, and beam-on time per fraction using RapidArc volumetric-modulated arc therapy (VMAT) vs. multifield sliding-window intensity-modulated radiation therapy (IMRT). Patients underwent computed tomography simulation with site-specific immobilization. Magnetic resonance imaging fusion and optical tracking were incorporated as clinically indicated. Treatment planning was performed using Eclipse v8.6 to generate sliding-window IMRT and 1-arc and 2-arc RapidArc plans. Dosimetric parameters used for target analysis were RTOG conformity index (CI{sub RTOG}), homogeneity index (HI{sub RTOG}), inverse Paddick Conformity Index (PCI), D{sub mean} and D5-D95. OAR sparing was analyzed in terms of D{sub max} and D{sub mean}. Treatment delivery was evaluated based on measured beam-on times delivered on a Varian Trilogy linear accelerator and recorded MU values. Dosimetric conformity, homogeneity, and OAR sparing were comparable between IMRT, 1-arc RapidArc and 2-arc RapidArc plans. Mean beam-on times {+-} SD for IMRT and 1-arc and 2-arc treatments were 10.5 {+-} 7.3, 2.6 {+-} 1.6, and 3.0 {+-} 1.1 minutes, respectively. Mean MUs were 3041, 1774, and 1676 for IMRT, 1-, and 2-arc plans, respectively. Although dosimetric conformity, homogeneity, and OAR sparing were similar between these techniques, SRS and SBRT fractions treated with RapidArc were delivered with substantially less beam-on time and fewer MUs than IMRT. The rapid delivery of SRS and SBRT with RapidArc improved workflow on the linac with these otherwise time-consuming treatments and limited the potential for intrafraction organ and patient motion, which can cause significant dosimetric errors. These clinically important advantages make image-guided RapidArc useful in the delivery of SRS and SBRT to intracranial and extracranial targets.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

An externally and internally deformable, programmable lung motion phantom

Purpose: Most clinically deployed strategies for respiratory motion management in lung radiotherapy (e.g., gating and tracking) use external markers that serve as surrogates for tumor motion. However, typical lung phantoms used to validate these strategies are based on a rigid exterior and a rigid or a deformable-interior. Such designs do not adequately represent respiration because the thoracic anatomy deforms internally as well as externally. In order to create a closer approximation of respiratory motion, the authors describe the construction and experimental testing of an externally as well as internally deformable, programmable lung phantom. Methods: The outer shell of a commercially available lung phantom (RS-1500, RSD, Inc.) was used. The shell consists of a chest cavity with a flexible anterior surface, and embedded vertebrae, rib-cage and sternum. A custom-made insert was designed using a piece of natural latex foam block. A motion platform was programmed with sinusoidal and ten patient-recorded lung tumor trajectories. The platform was used to drive a rigid foam “diaphragm” that compressed/decompressed the phantom interior. Experimental characterization comprised of determining the reproducibility and the external–internal correlation of external and internal marker trajectories extracted from kV x-ray fluoroscopy. Experiments were conducted to illustrate three example applications of the phantom—(i) validating the geometric accuracy of the VisionRT surface photogrammetry system; (ii) validating an image registration tool, NiftyReg; and (iii) quantifying the geometric error due to irregular motion in four-dimensional computed tomography (4DCT). Results: The phantom correctly reproduced sinusoidal and patient-derived motion, as well as realistic respiratory motion-related effects such as hysteresis. The reproducibility of marker trajectories over multiple runs for sinusoidal as well as patient traces, as characterized by fluoroscopy, was within 0.25 mm RMS error. The motion trajectories of internal and external radio-opaque markers as measured by fluoroscopy were found to be highly correlated (R > 0.95). Using the phantom, it was demonstrated that the motion trajectories of regions-of-interest on the surface as measured by VisionRT are highly consistent with corresponding fluoroscopically acquired surface marker trajectories, with RMS errors within 0.26 mm. Furthermore, it was shown that the trajectories of external and internal marker trajectories derived from NiftyReg deformation vector fields were within 1 mm root mean square errors comparing to trajectories obtained by segmenting markers from individual fluoro frames. Finally, it was shown that while 4DCT can be used to localize internal markers for sinusoidal motion with reasonable accuracy, the localization error increases significantly (by a factor of ∼2) in the presence of cycle-to-cycle variations that are observed in patient-derived respiratory motion. Conclusions: The authors have developed a realistic externally and internally deformable, programmable lung phantom that will serve as a valuable tool for clinical and investigational motion management studies in thoracic and abdominal radiation therapies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A method of surface marker location optimization for tumor motion estimation in lung stereotactic body radiation therapy

Purpose: Accurately localizing lung tumor localization is essential for high-precision radiation therapy techniques such as stereotactic body radiation therapy (SBRT). Since direct monitoring of tumor motion is not always achievable due to the limitation of imaging modalities for treatment guidance, placement of fiducial markers on the patient’s body surface to act as a surrogate for tumor position prediction is a practical alternative for tracking lung tumor motion during SBRT treatments. In this work, the authors propose an innovative and robust model to solve the multimarker position optimization problem. The model is able to overcome the major drawbacks of the sparse optimization approach (SOA) model. Methods: The principle-component-analysis (PCA) method was employed as the framework to build the authors’ statistical prediction model. The method can be divided into two stages. The first stage is to build the surrogate tumor matrix and calculate its eigenvalues and associated eigenvectors. The second stage is to determine the “best represented” columns of the eigenvector matrix obtained from stage one and subsequently acquire the optimal marker positions as well as numbers. Using 4-dimensional CT (4DCT) and breath hold CT imaging data, the PCA method was compared to the SOA method with respect to calculation time, average prediction accuracy, prediction stability, noise resistance, marker position consistency, and marker distribution. Results: The PCA and SOA methods which were both tested were on all 11 patients for a total of 130 cases including 4DCT and breath-hold CT scenarios. The maximum calculation time for the PCA method was less than 1 s with 64 752 surface points, whereas the average calculation time for the SOA method was over 12 min with 400 surface points. Overall, the tumor center position prediction errors were comparable between the two methods, and all were less than 1.5 mm. However, for the extreme scenarios (breath hold), the prediction errors for the PCA method were not only smaller, but were also more stable than for the SOA method. Results obtained by imposing a series of random noises to the surrogates indicated that the PCA method was much more noise resistant than the SOA method. The marker position consistency tests using various combinations of 4DCT phases to construct the surrogates suggested that the marker position predictions of the PCA method were more consistent than those of the SOA method, in spite of surrogate construction. Marker distribution tests indicated that greater than 80% of the calculated marker positions fell into the high cross correlation and high motion magnitude regions for both of the algorithms. Conclusions: The PCA model is an accurate, efficient, robust, and practical model for solving the multimarker position optimization problem to predict lung tumor motion during SBRT treatments. Due to its generality, PCA model can also be applied to other imaging guidance system whichever using surface motion as the surrogates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-258: Development of a New Patient Set-Up Monitoring System Using Force Sensing Resistor (FSR) Sensor for the Radiation Therapy

Purpose: The purpose of this work is to develop a new patient set-up monitoring system using force sensing resistor (FSR) sensors that can confirm pressure of contact surface and evaluate its feasibility. Methods: In this study, we focused on develop the patient set-up monitoring system to compensate for the limitation of existing optical based monitoring system, so the developed system can inform motion in the radiation therapy. The set-up monitoring system was designed consisting of sensor units (FSR sensor), signal conditioning devices (USB cable/interface electronics), a control PC, and a developed analysis software. The sensor unit was made by attaching FSR sensor and dispersing pressure sponge to prevent error which is caused by concentrating specific point. Measured signal from the FSR sensor was sampled to arduino mega 2560 microcontroller, transferred to control PC by using serial communication. The measured data went through normalization process. The normalized data was displayed through the developed graphic user interface (GUI) software. The software was designed to display a single sensor unit intensity (maximum 16 sensors) and display 2D pressure distribution (using 16 sensors) according to the purpose. Results: Changes of pressure value according to motion was confirmed by the developed set-up monitoring system. Very small movement such as little physical change in appearance can be confirmed using a single unit and using 2D pressure distribution. Also, the set-up monitoring system can observe in real time. Conclusion: In this study, we developed the new set-up monitoring system using FSR sensor. Especially, we expect that the new set-up monitoring system is suitable for motion monitoring of blind area that is hard to confirm existing optical system and compensate existing optical based monitoring system. As a further study, an integrated system will be constructed through correlation of existing optical monitoring system. This work was supported by the Industrial R&D program of MOTIE/KEIT. [10048997, Development of the core technology for integrated therapy devices based on real-time MRI guided tumor tracking] and the Mid-career Researcher Program (2014R1A2A1A10050270) through the National Research Foundation of Korea funded by the Ministry of Science, ICT&Future Planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-BRA-02: Automated Triplet Beam Orientation Optimization for MRI-Guided Co-60 Radiotherapy

Purpose: MRI guided Co-60 provides daily and intrafractional MRI soft tissue imaging for improved target tracking and adaptive radiotherapy. To remedy the low output limitation, the system uses three Co-60 sources at 120° apart, but using all three sources in planning is considerably unintuitive. We automate the beam orientation optimization using column generation, and then solve a novel fluence map optimization (FMO) problem while regularizing the number of MLC segments. Methods: Three patients—1 prostate (PRT), 1 lung (LNG), and 1 head-and-neck boost plan (H&NBoost)—were evaluated. The beamlet dose for 180 equally spaced coplanar beams under 0.35 T magnetic field was calculated using Monte Carlo. The 60 triplets were selected utilizing the column generation algorithm. The FMO problem was formulated using an L2-norm minimization with anisotropic total variation (TV) regularization term, which allows for control over the number of MLC segments. Our Fluence Regularized and Optimized Selection of Triplets (FROST) plans were compared against the clinical treatment plans (CLN) produced by an experienced dosimetrist. Results: The mean PTV D95, D98, and D99 differ by −0.02%, +0.12%, and +0.44% of the prescription dose between planning methods, showing same PTV dose coverage. The mean PTV homogeneity (D95/D5) was at 0.9360 (FROST) and 0.9356 (CLN). R50 decreased by 0.07 with FROST. On average, FROST reduced Dmax and Dmean of OARs by 6.56% and 5.86% of the prescription dose. The manual CLN planning required iterative trial and error runs which is very time consuming, while FROST required minimal human intervention. Conclusions: MRI guided Co-60 therapy needs the output of all sources yet suffers from unintuitive and laborious manual beam selection processes. Automated triplet orientation optimization is shown essential to overcome the difficulty and improves the dosimetry. A novel FMO with regularization provides additional controls over the number of MLC segments and treatment time. Varian Medical Systems; NIH grant R01CA188300; NIH grant R43CA183390.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Matlab Based LOCO

The LOCO algorithm has been used by many accelerators around the world. Although the uses for LOCO vary, the most common use has been to find calibration errors and correct the optics functions. The light source community in particular has made extensive use of the LOCO algorithms to tightly control the beta function and coupling. Maintaining high quality beam parameters requires constant attention so a relatively large effort was put into software development for the LOCO application. The LOCO code was originally written in FORTRAN. This code worked fine but it was somewhat awkward to use. For instance, the FORTRAN code itself did not calculate the model response matrix. It required a separate modeling code such as MAD to calculate the model matrix then one manually loads the data into the LOCO code. As the number of people interested in LOCO grew, it required making it easier to use. The decision to port LOCO to Matlab was relatively easy. It's best to use a matrix programming language with good graphics capability; Matlab was also being used for high level machine control; and the accelerator modeling code AT, [5], was already developed for Matlab. Since LOCO requires collecting and processing a relative large amount of data, it is very helpful to have the LOCO code compatible with the high level machine control, [3]. A number of new features were added while porting the code from FORTRAN and new methods continue to evolve, [7][9]. Although Matlab LOCO was written with AT as the underlying tracking code, a mechanism to connect to other modeling codes has been provided.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Fluid Dynamics Effects on Microstructure Prediction in Single Laser Tracks for Additive Manufacturing

The Laser Powder Bed Fusion Additive Manufacturing (LPBFAM) is one of the most important processes for the production of lightweight, cost-effective, complex, and high-performance enduse parts. At present, the cost and time associated with LPBFAM process development is very high due to a lack of fundamental process understanding. In this project, a multi-physics model was developed on a highly parallel open-source code, Truchas, with the ultimate goal of providing experimentally validated process maps for tailoring microstructure to achieve desired performance for LPBFAM. As a critical step towards fully LPBFAM modeling, modeling of single-track laser fusion (STLF) were conducted. Multi-physics simulations were conducted using Truchas for STLF by considering heat transfer, phase-change, fluid dynamics, surface tension phenomena, and evaporation. In order to assess the effect of the fluid dynamics on the solidification and ensuing microstructure, a heat transfer-and-solidification-only (HTS) model and a fully coupled heat-transfer-solidification and fluid-dynamics (HTS-FD) model were considered. In the HTS_FD model, the fluid flow was considered to be laminar while the molten alloy surface is assumed to be flat and non-deformable. This study is one the first attempt to understand the effect of the fluid flow on microstructure in STLF and LPBFAM. The results show that the fluid flow affects the solidification and ensuing microstructure in STLF. In order to validate Truchas for STLF and LPBFAM modeling, experimental data, which was obtained at GE Global Research (GEGR) for liquid pool shape and microstructure, were compared with those from numerical simulation results. For the keyhole regime case, numerical simulation results indicate that the melt-pool shape was typical to that of the conduction case with very large deviations from the measured melt-pool depths. Two sensitivity studies were conducted by varying the evaporation flux and the surface tension coefficient at fixed laser absorptivities of 0.5 and 0.3, respectively. The minimum value for the overall combined error, between the calculated values with the HTS_FD model and measured values for both the width and height of the melt-pool, was attained for a surface tension coefficient of -0.75e-4 N/m. An analytical solidification model for the columnar-to-equiaxed transition (CET) in rapid solidification was used to assess the microstructure variation within the melt-pool. From numerical simulation results, thermal gradient, G, and solidification velocity, V, were obtained in order to predict the microstructure type (e.g., dendritic, cellular). For the fluid-flow HTS_FD model, G was found to exhibit a maximum at the extremity of the solidified pool (i.e., at the free surface). By contrast, for HTS simulations, G was found to exhibit a maximum around the entire edge of the solidified pool. For the HTS_FD simulations, the minimum values of the cooling rate, GV, were found to be approximately half than their corresponding values for HTS simulations. By contrast with the min GV values, the maximum values for GV were found to occur for the HTS_FD simulations. Thus, HTS_FD simulations were found to exhibit a wider range of cooling rates than the HTS simulations. Concerning the solidification maps, the variation of the solidification velocity, V, as function of the thermal gradient, G, was obtained. It was found that the fluid flow model results (HTS_FD) exhibited an increased spread in the V(G) variation within the melt pool with respect to the HTS model results (without the fluid flow). A preliminary correlation model for the primary dendrite arm spacing (PDAS) based on a power law dependence on thermal gradient and solidification velocity was used to estimate PDAS from HTS and HTS_FD simulations. For low laser powers and low laser speeds, the PDAS obtained 6 with the fluid dynamics model (HTS_FD) was larger by more than 30% with respect to the PDAS calculated with the simple HTS model. In the second part of the project, fifty-seven simulations were conducted in order to obtain data on microstructure variables that can be used for process map development. In order to cover the entire processing space, thirty-three cases were selected STLF process simulations at six power levels and six laser scanning speeds. The 57 STLF process simulations were conducted as follows: 33 simulations with the heat-transfer-only (HTS) model and 24 simulations with the fluid flow model (HTS_FD). The solidification map data shows that for all simulations, a columnar dendritic microstructure would be expected. It was found that the minimum, average, and maximum thermal gradient exhibit exponential variations with respect to a process variable defined as the ratio between the power and square root of scan speed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Noise-canceling quantum feedback: Non-Hermitian dynamics with applications to state preparation and magic state distillation

Time-continuous quantum measurement allows for the tracking of a quantum system in real time via sequences of short, and individually weak, measurement intervals. Such measurements are necessarily invasive, imparting backaction to the system, and allowing the observer to update their state estimate based on stochastic measurement outcomes. Feedback control then involves real-time interventions by an observer, conditioned on the time-continuous measurement signal that they receive. Here, we consider here diffusive quantum trajectories and focus on the “noise-canceling” subset of feedback protocols that aim to minimize the degree of stochasticity in the dynamics. We derive such a class of feedback operations, showing that under the idealized assumptions of pure states, unit measurement efficiency, and zero time delay in implementing feedback operations, perfectly noise-canceling feedback always exists. We consider the resulting noise-canceled dynamics generated by an effective non-Hermitian Hamiltonian; while non-Hermitian Hamiltonians from continuous monitoring generally describe rare dynamics (accessible by costly postselection), the use of noise-canceling feedback here leads to non-Hermitian dynamics that occur deterministically. We demonstrate this via examples of entangled state preparation and stabilization. We then illustrate the potential for the application of noise cancellation to boost success rates in magic state distillation protocols. We show that adding feedback based on noise cancellation into a time-continuous five-to-one distillation protocol leads to higher probabilities of successful distillation across a range of input errors and increases the threshold on input errors for which the protocol is effective. Our results highlight the efficacy of noise-canceling feedback-aided protocols for quantum state preparation and stabilization tasks.

Karmakar, Tathagata [University of California, Ber↗
DOE OSTI2026

MOD-Amp System Design Spring: Spring 2026 – Georgetown University, SYSM–5620

High-energy lasers (HELs) play an important role in both national defense and scientific research. In defense applications, HELs are used for target detection, tracking, and engagement. In research environments, they support studies of extreme physical conditions relevant to fusion energy and plasma science. These systems depend on the amplification of light through stimulated emission of radiation, allowing optical energy to be increased to the levels required for operation. This amplification occurs when light passes through an energized gain medium that receives energy from an external optical or electrical source. To achieve the desired output, laser systems often use multiple amplification stages, including high-gain preamplifiers and lower-gain power or booster amplifiers. At Lawrence Livermore National Laboratory (LLNL) and other national laboratories, many large-aperture laser amplifier systems are aging and rely on system-specific hardware, obsolete technologies, and incomplete documentation. These legacy systems create challenges for maintenance, supportability, and long-term operation. Their lack of standardization also increases the difficulty of sustaining reliable performance over time. As this infrastructure continues to age, the likelihood of unplanned downtime grows, which can negatively affect both national security missions and scientific research programs that depend on dependable HEL capabilities. The purpose of this document is to demonstrate the application of systems engineering fundamentals and design thinking through the development of a laser amplifier case study. The proposed system concept is intended as an academic exercise and not as a finalized engineering design. As a result, the development presented in this document is incomplete and may contain technical assumptions or errors that would require further investigation before any real-world implementation.

42 ENGINEERING↗
DOE OSTI2015

Sparse-view proton computed tomography using modulated proton beams

Purpose: Proton imaging that uses a modulated proton beam and an intensity detector allows a relatively fast image acquisition compared to the imaging approach based on a trajectory tracking detector. In addition, it requires a relatively simple implementation in a conventional proton therapy equipment. The model of geometric straight ray assumed in conventional computed tomography (CT) image reconstruction is however challenged by multiple-Coulomb scattering and energy straggling in the proton imaging. Radiation dose to the patient is another important issue that has to be taken care of for practical applications. In this work, the authors have investigated iterative image reconstructions after a deconvolution of the sparsely view-sampled data to address these issues in proton CT. Methods: Proton projection images were acquired using the modulated proton beams and the EBT2 film as an intensity detector. Four electron-density cylinders representing normal soft tissues and bone were used as imaged object and scanned at 40 views that are equally separated over 360°. Digitized film images were converted to water-equivalent thickness by use of an empirically derived conversion curve. For improving the image quality, a deconvolution-based image deblurring with an empirically acquired point spread function was employed. They have implemented iterative image reconstruction algorithms such as adaptive steepest descent-projection onto convex sets (ASD-POCS), superiorization method–projection onto convex sets (SM-POCS), superiorization method–expectation maximization (SM-EM), and expectation maximization-total variation minimization (EM-TV). Performance of the four image reconstruction algorithms was analyzed and compared quantitatively via contrast-to-noise ratio (CNR) and root-mean-square-error (RMSE). Results: Objects of higher electron density have been reconstructed more accurately than those of lower density objects. The bone, for example, has been reconstructed within 1% error. EM-based algorithms produced an increased image noise and RMSE as the iteration reaches about 20, while the POCS-based algorithms showed a monotonic convergence with iterations. The ASD-POCS algorithm outperformed the others in terms of CNR, RMSE, and the accuracy of the reconstructed relative stopping power in the region of lung and soft tissues. Conclusions: The four iterative algorithms, i.e., ASD-POCS, SM-POCS, SM-EM, and EM-TV, have been developed and applied for proton CT image reconstruction. Although it still seems that the images need to be improved for practical applications to the treatment planning, proton CT imaging by use of the modulated beams in sparse-view sampling has demonstrated its feasibility.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Holistic Measurement Driven Resilience: Combining Operational Fault and Failure Measurements and Fault Injection for Quantifying Fault Detection, Propagation and Impact. Final report

For HPC systems to date, application resilience to faults and failures has been accomplished by the brute- force method of checkpoint/restart, which allows an application to make forward progress in the face of system and application faults, errors, and failures independent of root cause or end result. It has remained the primary resilience mechanism because we lack a way to identify faults and anticipate consequences early enough to take meaningful mitigating action. However, checkpoint/restart implementations put a tremendous burden on system resources and on the applications themselves and is becoming less feasible at scale. Because we have not yet operated at scales at which checkpoint/restart fails to provide forward progress, despite increasing costs, vendors have had little motivation to provide the instrumentation necessary for early identification of faults and failures. However, as we move from petascale to exascale, component mean time to failure (MTTF) will render the existing techniques ineffectual and/or too expensive. Furthermore, fault recovery mechanisms such as failover and/or error correction introduce performance inconsistency. Instrumentation allowing early indication of problems and tools to enable use of such information by systems, operating systems, and applications offer an alternative, more scalable and less costly solution. In the HMDR project, we built on our experience and expertise developed and accumulated over years of research on design, monitoring, measurement, and assessment of resilient computing systems. Analysis of field data on the current and past generations of extreme-scale systems revealed several challenges that, if not addressed in increasingly larger and more complex systems, may hinder the effectiveness of future exascale computing systems. Specifically, i) file systems and interconnects in current-generation large-scale systems already operate at the margins of resiliency, including consistent performance, and may not scale to larger deployments; ii) automated, software-based failover mechanisms are frequently inadequate and can introduce wider failures, such that failures during recovery may lead to system/application failures, including system-wide outages; and iii) silent data corruption represents a critical fault mode and will require efficient detection mechanisms if next-generation applications are to take full advantage of exascale hardware. To address the above challenges, we assembled a team of world-renowned experts in resilient extreme- scale computing from the University of Illinois (Electrical and Computer Engineering, Computer Science, and NCSA), SNL, LANL, NERSC, and Cray. Our team includes representatives from centers that house many of the largest HPC resources in the world, both today and over the coming years. The team has a unique track record of research in i) system and application failure characterization based on the analysis of field data, ii) data-driven design of fault/error detection mechanisms, and iii) experimental characterization of system/application resiliency. The team includes system owners/operators who provide continuous data collection and access and ensure installation of appropriate analysis tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Muon Neutrino to Electron Neutrino Oscillation in NO$\nu$A

NOvA is a long-baseline neutrino oscillation experiment optimized for electron neu- trino ( e) appearance in the NuMI beam, a muon neutrino ( $\nu$ μ ) source at Fermilab. It consists of two functionally identical, nearly fully-active liquid-scintillator tracking calorimeters. The near detector (ND) at Fermilab is used to study the neutrino beam spectrum and composition before oscillation, and measure background rate to the $\nu$ e appearance search. The far detector, 810 km away in Northern Minnesota, observes the oscillated beam and is used to extract oscillation parameters from the data. NO$\nu$ A's long baseline, combined with the ability of the NuMI beam to operate in the anti- neutrino mode, makes NO$\nu$ A sensitive to the last unmeasured parameters in neutrino oscillations- mass hierarchy, CP violation and the octant of mixing angle 23. This thesis presents the search for e appearance in the first data collected by the NO A detectors from October 2013 till May 2015. Studies of the NuMI neutrino data collected in the NO A near detector are also presented, which show large discrepancies between the ND simulation and data. Muon- removed electron (MRE) events, constructed by replacing the muon in $\nu$ μ charged cur- rent interactions by a simulated electron, are used to correct the far detector $\nu$e appearance prediction for these discrepancies. In the analysis of the first data, a total of 6 e candidate events are observed in the far detector on a background of 1, a 3.46 σ excess, which is interpreted as strong evidence for e appearance. The results are consistent with our expectation, based on constraints from other neutrino oscillation experiments. The result presented here differs from the officially published $\nu$ e appearance result from the NO A experiment where the systematic error is assumed to cover the MRE correction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Evaluation of a Sulfate Solubility Limit Greater Than 0.65 Weight Percent in Sludge Batch 10 Glasses

Previous laboratory-scale crucible testing with batch chemicals confirmed that the sulfate (SO 4 2- ) limit for Sludge Batch 10 (SB10) was 0.65 weight percent (wt.%) in glass. This limit signifies that 0.65 wt.% SO 4 2- can be retained in the glass without the formation of a sulfate phase. The Defense Waste Processing Facility imposes this constraint in the Material Tracking Program. Based on preliminary calculations to support the Material Tracking Program, it was anticipated that transfer volumes of the monosodium titanate/sludge solids (MST/SS) stream from the Salt Waste Processing Facility (SWPF) may need to be reduced to maintain projected sulfate concentrations below 0.65 wt.% in glass. Savannah River Mission Completion requested that the Savannah River National Laboratory perform additional sulfate testing to determine whether a sulfate solubility limit greater than 0.65 wt.% is feasible for SB10, which could allow for higher transfer volumes of the MST/SS stream. This report documents the results of the testing at higher sulfate concentrations for the glass composition region defined by the most recent SB10 projection (November 2022) and Frits 473 and 625. Frit 473 was recommended for SB10 and Frit 625 was used during SB9 processing and the SB9 to SB10 transition. A total of twenty-one glass compositions were developed based on the expected compositional variables, which include sludge-only (SO) and coupled processing with the SWPF, waste loading (WL), and frit composition. The target sulfate concentrations were varied from 0.65-0.85 wt.% at 32 and 40% WL. Each glass was prepared from reagent grade chemicals and melted at 1150 °C. Visual observations were used to confirm the presence of a sulfate salt phase on the cooled glass surfaces. Representative samples of each glass were submitted for chemical composition analysis by inductively coupled plasma-optical emission spectroscopy and Cs analysis by inductively coupled plasma-mass spectrometry. Overall the majority of mean measured values are consistent with the target values for each major oxide of interest with less than 5% error. The percent errors for the measured SO 4 2- concentrations are generally less than 10%, which is comparable with previous sulfate solubility study measurements and acceptable. Only the SO glasses based on Frit 625 formed a sulfate phase at a 0.80 wt.% SO 4 2- target concentration at both 32 and 40% WL. The remainder of the glasses did not form a sulfate layer. Due to the formation of the sulfate phase, the limit is conservatively set at 0.70 wt.% based on the measured sulfate concentrations of 0.71 wt.% and 0.75 wt.% for these two glasses. None of the SO or coupled operation glasses based on Frit 473 formed a sulfate salt phase, which supports a sulfate limit of 0.80 wt.%. The following SO 4 2- concentration limits are recommended during SO and coupled SB10 processing: (1) 0.70 wt.% during processing with Frit 625, and (2) 0.80 wt.% during processing with Frit 473.

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

Fault-Tolerant Decentralized Control for Large-Scale Inverter-Based Resources for Active Power Tracking

Integration of inverter-based resources (IBRs) which lack the intrinsic characteristics such as the inertial response of the traditional synchronous-generator (SG)-based sources presents a new challenge in the form of analyzing the grid stability under their presence. While the dynamic composition of IBRs differs from that of the SGs, the control objective remains similar in terms of tracking the desired active power. This letter presents a decentralized primal-dual-based fault-tolerant control framework for the power allocation in IBRs. Overall, a hierarchical control algorithm is developed with a lower level addressing the current control and the parameter estimation for the IBRs and the higher level acting as the reference power generator to the low level based on the desired active power profile. The decentralized network-based algorithm adaptively splits the desired power between the IBRs taking into consideration the health of the IBRs transmission lines. The proposed framework is tested through a simulation on the network of IBRs and the high-level controller performance is compared against the existing framework in the literature. The proposed algorithm shows significant performance improvement in the magnitude of power deviation and settling time to the nominal value under faulty conditions as compared to the algorithm in the literature.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2018

Human performance benefits gained by dynamic instructions compared to smart PDFs

Human performance benefits gained by dynamic instructions compared to smart PDFs In recent years multiple nuclear utilities have deployed an electronic work management solution to replace the traditional paper-based processes. There are many benefits of electronic work management solutions compared to paper processes, such as a more streamlined review and approval process of work packages and a better way to track and resolve any potential bottlenecks causing delays in the review process. These benefits have been recognized where electronic work management solutions have been successfully deployed. However the utilities also hoped the new solution would help the field worker conduct his/her job more efficiently, correctly, and safely. Unfortunately, the industry is now realizing that the deployed solutions did not fully help them reach that goal. The electronic work management solutions present the work instructions as an electronic copy (i.e., a portable document format (PDF)) of the paper instruction. The PDF versions used in these electronic work management solutions have some ability for data input, notes, and mark-ups, but are unable to provide many additional benefits compared to using paper. In fact, conducting work using a PDF on a tablet can be more cumbersome than with a paper instruction. For example, flipping back and forth in a paper instruction can be time consuming, but it is not at all as frustrating as scrolling back and forth on a tablet or continuously having to navigate between tabs. The researchers at Idaho National Laboratory (INL) approached the next step beyond smart PDFs from a human factors perspective. In close collaboration with multiple utilities they investigated how to present instruction content and how to design the interaction with the instruction to improve the workers performance and efficiency. Their solution is a dynamic presentation of the instruction that guides the worker through the correct path as well as ensuring the worker is on the correct component. The solution also automatically validates input and provides guidance if unexpected conditions are identified. A well-designed solution has the opportunity to vastly reduce the risk of human errors in the field, such as preventing actions on incorrect component or unintentionally skipping steps. The human factors researchers at INL leveraged their expertise, existing human performance tools, and a great amount of involvement from field workers in multiple utilities to develop and evaluate design concepts that will provide the human performance and efficiency benefits that the electronic work management solutions fall short to deliver. This paper will describe in detail how these new design concepts improve human performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Remote Radiation Sensing Using Aerial and Ground Platforms

Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗