Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Programming and software”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Automated Cache Performance Analysis And Optimization

While there is no lack of performance counter tools for coarse-grained measurement of cache activity, there is a critical lack of tools for relating data layout to cache behavior to application performance. Generally, any nontrivial optimizations are either not done at all, or are done ”by hand” requiring significant time and expertise. To the best of our knowledge no tool available to users measures the latency of memory reference instructions for particular addresses and makes this information available to users in an easy-to-use and intuitive way. In this project, we worked to enable the Open|SpeedShop performance analysis tool to gather memory reference latency information for specific instructions and memory addresses, and to gather and display this information in an easy-to-use and intuitive way to aid performance analysts in identifying problematic data structures in their codes. This tool was primarily designed for use in the supercomputer domain as well as grid, cluster, cloud-based parallel e-commerce, and engineering systems and middleware. Ultimately, we envision a tool to automate optimization of application cache layout and utilization in the Open|SpeedShop performance analysis tool. To commercialize this software, we worked to develop core capabilities for gathering enhanced memory usage performance data from applications and create and apply novel methods for automatic data structure layout optimizations, tailoring the overall approach to support existing supercomputer and cluster programming models and constraints. In this Phase I project, we focused on infrastructure necessary to gather performance data and present it in an intuitive way to users. With the advent of enhanced Precise Event-Based Sampling (PEBS) counters on recent Intel processor architectures and equivalent technology on AMD processors, we are now in a position to access memory reference information for particular addresses. Prior to the introduction of PEBS counters, cache behavior could only be measured reliably in the aggregate across tens or hundreds of thousands of instructions. With the newest iteration of PEBS technology, cache events can be tied to a tuple of instruction pointer, target address (for both loads and stores), memory hierarchy, and observed latency. With this information we can now begin asking questions regarding the efficiency of not only regions of code, but how these regions interact with particular data structures and how these interactions evolve over time. In the short term, this information will be vital for performance analysts understanding and optimizing the behavior of their codes for the memory hierarchy. In the future, we can begin to ask how data layouts might be changed to improve performance and, for a particular application, what the theoretical optimal performance might be. The overall benefit to be produced by this effort was a commercial quality easy-to- use and scalable performance tool that will allow both beginner and experienced parallel programmers to automatically tune their applications for optimal cache usage. Effective use of such a tool can literally save weeks of performance tuning effort. Easy to use. With the proposed innovations, finding and fixing memory performance issues would be more automated and hide most to all of the performance engineer expertise ”under the hood” of the Open|SpeedShop performance tool. One of the biggest public benefits from the proposed innovations is that it makes performance analysis more usable to a larger group of application developers. Intuitive reporting of results. The Open|SpeedShop performance analysis tool has a rich set of intuitive, yet detailed reports for presenting performance results to application developers. Our goal was to leverage this existing technology to present the results from our memory performance addition to Open|SpeedShop. Suitable for experts as well as novices. Application performance is getting more difficult to measure as the hardware platforms they run on become more complicated. This makes life difficult for the application developer, in that they need to know more about the hardware platform, including the memory system hierarchy, in order to understand the performance of their application. Some application developers are comfortable in that sceario, while others want to do their scientific research and not have to understand all the nuances in the hardware platform they are running their application on. Our proposed innovations were aimed to support both experts and novice performance analysts. Useful in many markets. The enhancement to Open|SpeedShop would appeal to a broader market space, as it will be useful in scientific, commercial, and cloud computing environments. Our goal was to use technology developed initially at the and Lawrence Livermore National Laboratory combined with the development and commercial software experience of the Argo Navis Technologies, LLC (ANT) to form a powerful combination to delivery these objectives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

TU-H-BRC-01: A New Fabrication Method for Secondary Skin Collimation Using 3D Scanner and 3D Printer

Purpose: Using secondary skin collimation (SSC) to protect the critical organ surrounding the tumors is always desirable for electron and/or ortho-voltage treatments. However, the time-consuming and labor-intensive fabrication processes of current method (manually fabricate the lead mask) restrict the general use of SSC. The aim of this study is to develop a new SSC fabrication method using 3D scanning and 3D printing technologies in order to largely decrease the human labor involvement and fabrication time, also improve the fabrication accuracy. Methods: First, the patient surface was scanned with a 3D scanner (Structure Sensor, Occipital, CO). The raw scan data was subsequently transferred to a 3D modeling software (Rhinoceros, Ver. 5.0, McNeel North America, Seattle, WA). The tumor contour was then digitized and shielding region was determined by clinicians in the same software. The corresponding SCC conformed to the skin surface was then automatically generated by the software with the proper shielding thickness. The shell of the SCC (with hollow inside) was consequently printed by a 3D printer (Lulzbot TAZ, Aleph Objects, CO) using plastic material. Finally, the hollow mold of SCC was filled up with a melted cerrobend alloy. Once the cerrobend alloy cooled down, the fabrication of SCC was accomplished. Results: The results indicated the proposed method can achieve a much shorter time on making a SCC compared with tradition fabrication method. The processes of making a skin contour model for patients have been eliminated with the new method. SCC created by the new method possessed better accuracy and better conformality to patient’s skin contours. Conclusion: In this study, we have demonstrated a new method for the SCC fabrication. It is anticipated that our method can be an alternative way to replace conventional manual-based methods for electron and/or ortho-voltage SCC fabrication. This research was supported by the Global Ph.D. Fellowship Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (Grant No. 2015H1A2A1034071), Republic of Korea.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

WE-D-BRA-04: Online 3D EPID-Based Dose Verification for Optimum Patient Safety

Purpose: To develop an online 3D dose verification tool based on EPID transit dosimetry to ensure optimum patient safety in radiotherapy treatments. Methods: A new software package was developed which processes EPID portal images online using a back-projection algorithm for the 3D dose reconstruction. The package processes portal images faster than the acquisition rate of the portal imager (∼ 2.5 fps). After a portal image is acquired, the software seeks for “hot spots” in the reconstructed 3D dose distribution. A hot spot is in this study defined as a 4 cm{sup 3} cube where the average cumulative reconstructed dose exceeds the average total planned dose by at least 20% and 50 cGy. If a hot spot is detected, an alert is generated resulting in a linac halt. The software has been tested by irradiating an Alderson phantom after introducing various types of serious delivery errors. Results: In our first experiment the Alderson phantom was irradiated with two arcs from a 6 MV VMAT H&N treatment having a large leaf position error or a large monitor unit error. For both arcs and both errors the linac was halted before dose delivery was completed. When no error was introduced, the linac was not halted. The complete processing of a single portal frame, including hot spot detection, takes about 220 ms on a dual hexacore Intel Xeon 25 X5650 CPU at 2.66 GHz. Conclusion: A prototype online 3D dose verification tool using portal imaging has been developed and successfully tested for various kinds of gross delivery errors. The detection of hot spots was proven to be effective for the timely detection of these errors. Current work is focused on hot spot detection criteria for various treatment sites and the introduction of a clinical pilot program with online verification of hypo-fractionated (lung) treatments.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

SciPy 1.0: fundamental algorithms for scientific computing in Python

SciPy is an open-source scientific computing library for the Python programming language. Since its initial release in 2001, SciPy has become a de facto standard for leveraging scientific algorithms in Python, with over 600 unique code contributors, thousands of dependent packages, over 100,000 dependent repositories and millions of downloads per year. In this work, we provide an overview of the capabilities and development practices of SciPy 1.0 and highlight some recent technical developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Detection and quantitation of right ventricular reversible perfusion defects by stress SPECT myocardial perfusion imaging: A proof-of-principle study

Background: In patients with right dominant coronary circulation, the right ventricular (RV) myocardium and the inferior region of the left ventricular (LV) myocardium share a common source of blood flow. We hypothesized that stress/rest SPECT myocardial perfusion imaging (MPI) could detect reversible perfusion defects in the RV in some patients with LV inferior wall perfusion abnormalities. Material and Methods: We identified 2 groups of patients with LV inferior wall perfusion defects (with or without defects in other regions of LV myocardium) from our database. Patients in group 1 (n = 17) had reversible perfusion defects in the RV free wall by visual analysis, while patients in group 2 (n = 17) did not. The images were processed with filtered back projection and, separately, with iterative reconstruction. The images were then re-processed using an automated quantitative software that is specifically designed to include the RV in the region of interest. Results: There were 76% men in group 1 and 94% in group 2 (P <0.05). The mean age was 65±20 in group 1 vs. 63±18 years in group 2 (P < 0.05). The stress type was exercise in 30% in group 1 and 35% in group 2, with the remaining patients studied with pharmacological stress testing (P = NS). The presence of RV reversible perfusion defects using filtered back projection was more evident in 13 patients (75%), while it was better seen with iterative reconstruction in 4 patients (25%). By automated analysis, the RV reversible perfusion defect size was 19 ± 14% of RV myocardium. Conclusion: This proof-of-principle study demonstrates that reversible RV perfusion defects suggestive of ischemia can be detected by SPECT myocardial perfusion imaging in some patients with LV inferior ischemia by visual analysis and can be quantitated by automated programs. Further studies on the diagnostic and prognostic relevance of assessing RV ischemia on SPECT MPI are needed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

MCNP6 Monte Carlo Code Optimization [Slides]

The MCNP6.1 Monte Carlo code, released in 2013, offers many features not available in previous versions, but runs 20-30% slower for most problems and 2-5x slower for some problems. The MCNP 2020 initiative was established in mid-2013 to address code performance and other issues, with 3 basic components: (1) Improve MCNP performance, with the specific goal of achieving a 2x performance gain within 2 years; (2) Upgrade the core MCNP software and algorithms in an evolutionary fashion, with the goal of achieving more agile and sustainable code; (3) Prepare for future exascale systems that will include massive parallelism and limited memory per core, with the goal of providing flexible and adaptable coding for possible exascale systems. This talk reviews the initial performance improvements to MCNP6.1 that have been incorporated into the 2014 update, MCNP6.1.1. The performance improvements to date have included both classic code optimizations and algorithm improvements. Testing on a variety of problems has demonstrated that the performance improvements were effective, yielding speedups by factors of 1.2x - 4x compared to MCNP6.1, depending on the type of problem. For criticality problems, MCNP6.1.1 runs 1.5x - 1.7x faster than MCNP6.1. Much more work is planned for improving MCNP6 performance, structure, and algorithms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

MO-D-213-07: RadShield: Semi- Automated Calculation of Air Kerma Rate and Barrier Thickness

Purpose: To develop the first Java-based semi-automated calculation program intended to aid professional radiation shielding design. Air-kerma rate and barrier thickness calculations are performed by implementing NCRP Report 147 formalism into a Graphical User Interface (GUI). The ultimate aim of this newly created software package is to reduce errors and improve radiographic and fluoroscopic room designs over manual approaches. Methods: Floor plans are first imported as images into the RadShield software program. These plans serve as templates for drawing barriers, occupied regions and x-ray tube locations. We have implemented sub-GUIs that allow the specification in regions and equipment for occupancy factors, design goals, number of patients, primary beam directions, source-to-patient distances and workload distributions. Once the user enters the above parameters, the program automatically calculates air-kerma rate at sampled points beyond all barriers. For each sample point, a corresponding minimum barrier thickness is calculated to meet the design goal. RadShield allows control over preshielding, sample point location and material types. Results: A functional GUI package was developed and tested. Examination of sample walls and source distributions yields a maximum percent difference of less than 0.1% between hand-calculated air-kerma rates and RadShield. Conclusion: The initial results demonstrated that RadShield calculates air-kerma rates and required barrier thicknesses with reliable accuracy and can be used to make radiation shielding design more efficient and accurate. This newly developed approach differs from conventional calculation methods in that it finds air-kerma rates and thickness requirements for many points outside the barriers, stores the information and selects the largest value needed to comply with NCRP Report 147 design goals. Floor plans, parameters, designs and reports can be saved and accessed later for modification and recalculation. We have confirmed that this software accurately calculates air-kerma rates and required barrier thicknesses for diagnostic radiography and fluoroscopic rooms.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

SU-F-T-262: Commissioning Varian Portal Dosimetry for EPID-Based Patient Specific QA in a Non-Aria Environment

Purpose: Development of an in-house program facilitates a workflow that allows Electronic Portal Imaging Device (EPID) patient specific quality assurance (QA) measurements to be acquired and analyzed in the Portal Dosimetry Application (Varian Medical Systems, Palo Alto, CA) using a non-Aria Record and Verify (R&V) system (MOSAIQ, Elekta, Crawley, UK) to deliver beams in standard clinical treatment mode. Methods: Initial calibration of an in-house software tool includes characterization of EPID dosimetry parameters by importing DICOM images of varying delivered MUs to determine linear mapping factors in order to convert image pixel values to Varian-defined Calibrated Units (CU). Using this information, the Portal Dose Image Prediction (PDIP) algorithm was commissioned by converting images of various field sizes to output factors using the Eclipse Scripting Application Programming Interface (ESAPI) and converting a delivered configuration fluence to absolute dose units. To verify the algorithm configuration, an integrated image was acquired, exported directly from the R&V client, automatically converted to a compatible, calibrated dosimetric image, and compared to a PDIP calculated image using Varian’s Portal Dosimetry Application. Results: For two C-Series and one TrueBeam Varian linear accelerators, gamma comparisons (global 3% / 3mm) of PDIP algorithm predicted dosimetric images and images converted via the inhouse system demonstrated agreement for ≥99% of all pixels, exceeding vendor-recommended commissioning guidelines. Conclusion: Combinations of a programmatic image conversion tool and ESAPI allow for an efficient and accurate method of patient IMRT QA incorporating a 3rd party R&V system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Convolution Correction Factor Adjustments on Static PRA Models for Event Assessment

Current probabilistic risk assessment (PRA) models will contain cut sets that have time-constrained basic events such as switched components, recovery of failed components, and failures in time. For these cut sets that have multiple time-constrained events that interact (e.g., offsite power recover at the same time as multiple diesel generators fail to run), there is the possibility that the minimal cut set quantified result is too large. Since current PRA models typically assume that the accident scenario starts at “time zero,” the fails-to-run basic events use this starting time as the basis for their respective reliability models. Over time, the standard practice was to recognize this issue and accept it as part of the conservatism built into the event tree/fault tree PRA models. However, now that these same PRA models are being used in a regulatory forum, this conservatism needs to be addressed. Further, this issue exists not only for the baseline PRA results, but for deficiencies that have been observed or those that have the potential of occurrence, whether the deficiency is complete failure or just a degraded condition. We have found that this conservatism becomes very important when evaluating scenarios such as station blackout sequences and other loss of offsite power sequences. This paper will address the theory and implemented process used by the Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE) PRA software. Starting with SAPHIRE version 8, the process that was developed is to calculate a “convolution correction factor” based on three distinct calculations. The first calculation is the cut set based frequency calculated via the standard event tree/fault tree PRA model (assuming all run failures start at time equal to zero). The second calculation is to perform the exact time-constrained calculation by integrating the component failure models over the time period of interest (e.g., time to recover, a mission time). For example, this exact calculation for the case of the diesel generators is to convolve their failure models with the offsite power recovery model to obtain the exact frequency. The third calculation is to produce the convolution correction factor by dividing the exact frequency calculation by the standard cut set frequency. This convolution correction factor is then applied to the PRA cut sets in order to remove the time-constrained conservatism. The process used to obtain the convolution correction factor for the baseline PRA is important; however, in this paper we also describe what happens when there is degradation of a component or a potential deficiency. These issues become very important in order to ensure that the convolution correction factor is re-calculated correctly. If the correction factor does not get re-calculated and the same baseline convolution correction factor is used, the final analysis result could be over conservative or under conservative. This paper will discuss how these issues are addressed and provide examples on how SAPHIRE automatically makes the necessary adjustments. MathCad verification results can be computed using the provided equations and input data.

99 GENERAL AND MISCELLANEOUS↗
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 OSTI2020

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

SU-E-T-437: Four-Dimensional Treatment Planning for Lung VMAT-SBRT

Purpose: To assess optimal treatment planning approach of Volumetric Modulated Arc Therapy for lung Stereotactic Body Radiation Therapy (VMAT-SBRT). Methods: Subjects were 10 patients with lung cancer who had undergone 4DCT. The internal target volume (ITV) volume ranged from 2.6 to 16.5cm{sup 3} and the tumor motion ranged from 0 to 2cm. From 4DCT, which was binned into 10 respiratory phases, 4 image data sets were created; maximum intensity projection (MIP), average intensity projection (AIP), AIP with the ITV replaced by 0HU (RITV-AIP) and RITV-AIP with the planning target volume (PTV) minus the internal target volume was set to −200 HU (HR-AIP). VMAT-SBRT plans were generated on each image set for a patient. 48Gy was prescribed to 95% of PTV. The plans were recalculated on all phase images of 4DCT and the dose distributions were accumulated using a deformable image registration software MIM Maestro™ as the 4D calculated dose to the gross tumor volume (GTV). The planned dose to the ITV and 4D calculated dose to the GTV were compared. Results: In AIP plan, 10 patients average of all dose parameters (D1%, D-mean, and D99%) discrepancy were 1Gy or smaller. MIP and RITV-AIP plans resulted in having common tendency and larger discrepancy than AIP plan. The 4D dose was lower than the planned dose, and 10 patients average of all dose parameters discrepancy were in range 1.3 to 2.6Gy. HR-AIP plan had the largest discrepancy in our trials. 4D calculated D1%, D-mean, and D99% were resulted in 3.0, 4.1, and 6.1Gy lower than the expected in plan, respectively. Conclusion: For all patients, the dose parameters expected in AIP plan approximated to 4D calculated. Using AIP image set seems optimal treatment planning approach of VMAT-SBRT for a mobile tumor. Funding Support: This work was supported by the Japan Society for the Promotion of Science Core-to-Core program (No. 23003)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-F-BRA-12: End-User Oriented Tools and Procedures for Testing Brachytherapy TPSs Employing MBDCAs

Purpose: To develop user-oriented tools for commissioning and dosimetry testing of {sup 192}Ir brachytherapy treatment planning systems (TPSs) employing model based dose calculation algorithms (MBDCAs). Methods: A software tool (BrachyGuide) has been developed for the automatic generation of MCNP6 input files from any CT based plan exported in DICOM RT format from Elekta and Varian TPSs. BrachyGuide also facilitates the evaluation of imported Monte Carlo (MC) and TPS dose distributions in terms of % dose differences and gamma index (CT overlaid colormaps or relative frequency plots) as well as DVHs and related indices. For users not equipped to perform MC, a set of computational models was prepared in DICOM format, accompanied by treatment plans and corresponding MCNP6 generated reference data. BrachyGuide can then be used to compare institutional and reference data as per TG186. The model set includes a water sphere with the MBDCA WG {sup 192}Ir source placed centrically and in two eccentric positions, a water sphere with cubic bone and lung inhomogeneities and a five source dwells plan, and a patient equivalent model with an Accelerated Partial Breast Irradiation (APBI) plan. Results: The tools developed were used for the dosimetry testing of the Acuros and ACE MBDCAs implemented in BrachyVision v.13 and Oncentra Brachy v.4.5, respectively. Findings were consistent with previous results in the literature. Besides points close to the source dwells, Acuros was found to agree within type A uncertainties with the reference MC results. Differences greater than MC type A uncertainty were observed for ACE at distances >5cm from the source dwells and in bone. Conclusion: The tools developed are efficient for brachytherapy MBDCA planning commissioning and testing. Since they are appropriate for distribution over the web, they will be put at the AAPM WG MBDCA’s disposal. Research co-financed by the ESF and Greek funds. NSRF operational Program: Education and Lifelong Learning Investing in Knowledge Society-Aristeia. Varian Medical Systems and Nucletron, an Elekta company provided access to TPSs for research purposes. Miss Peppa was supported by IKY-fellowships of excellence for postgraduate studies in Greece,Siemens Program.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Custom surface reflectance, shade mask, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study (2025)

This dataset contains land surface reflectance estimates and additional derived products generated from NEON Imaging Spectrometer (NIS) data collected in the Upper Gunnison river basin during June and July of 2025. Data was collected over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). These products were derived from radiance and LiDAR data collected by the NEON Airborne Observation Platform (AOP) campaign funded by the Colorado Headwaters Ecological Spectroscopy Study (CHESS) (doi:10.15485/3017965). Products include per-pixel surface reflectance (rfl) and reflectance uncertainty (rfl_unc), observational data (obs), canopy equivalent water thickness (ewt), and shade masks. Atmospheric correction was performed per flightline using the ISOFIT (Imaging Spectrometer Optimal FITting) optimal estimation framework to estimate surface reflectance and the associated per-band reflectance uncertainty. Reflectance retrievals achieved a mean absolute error of 1.5% across diverse validation surfaces (see validation report.pdf). Equivalent water thickness was calculated from surface reflectance using the Beer–Lambert absorption of liquid water. Shade masks were generated based on the geometry between the sun angle, ground surface, and sensor at the time of flight. Data products are provided per-flightline and as mosaics for each domain. Flightline data products are provided as ENVI-formatted binary files (rfl, rfl_unc, ewt) and GeoTIFFs (shade). Reflectance and uncertainty mosaics are provided as tiled NetCDFs, while all other mosaicked products are provided as cloud-optimized GeoTIFFs. These formats are supported by common geospatial software (e.g., QGIS, ArcGIS, ENVI) and programmatic libraries in Python (e.g., rasterio, xarray, spectral, netCDF4) and R (e.g., terra, ncdf4). Processing workflows were designed to be equivalent to those used to generate the 2018 CHESS campaign airborne imaging spectroscopy data products (doi:10.15485/3013527). All outputs were co-registered to a common spatial grid to support time series analyses. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgment: Data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). Computational research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004) and was funded by EMIT Extended Mission Phase E Science.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2016

RAVEN User Manual, Revision 4

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The initial development was aimed to provide dynamic risk analysis capabilities to the Thermo-Hydraulic code RELAP-7, currently under development at the Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose probabilistic and uncertainty quantification platform, capable to agnostically communicate with any system code. This agnosticism includes providing Application Programming Interfaces (APIs). These APIs are used to allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by inputs files or via python interfaces. RAVEN is capable of investigating the system response, and investigating the input space using Monte Carlo, Grid, or Latin Hyper Cube sampling schemes, but its strength is focused to- ward system feature discovery, such as limit surfaces, separating regions of the input space leading to system failure, using dynamic supervised learning techniques. The development of RAVEN has started in 2012, when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework became stronger. RAVEN principal assignment is to provide the necessary software and algorithms in order to employ the concept developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just the individuation of the frequency of an event potentially leading to a system failure, but the closeness (or not) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. for an important process such as peak pressure in a pipe) is exceeded under certain conditions. The initial development of RAVEN has been focused on providing dynamic risk assess- ment capability to RELAP-7, currently under development at the INL and, likely, future replacement of the RELAP5-3D code. Most the capabilities that have been implemented having RELAP-7 as principal focus are easily deployable for other system codes. For this reason, several side activates are currently ongoing for coupling RAVEN with soft- ware such as RELAP5-3D, etc. The aim of this document is the explanation of the input requirements, focalizing on the input structure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Insert Modeling in UNF ST&DARDS

The Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS) is a software tool that integrates a used nuclear fuel (UNF) or spent nuclear fuel (SNF) relational database and key analysis capabilities to simplify and automate numerous UNF management and fuel cycle–related activities. UNF-ST&DARDS is being developed for the US Department of Energy’s Office of Nuclear Energy Spent Fuel and Waste Disposition program. UNF-ST&DARDS provides an integrated framework that uses advanced modeling and simulation to predict the behavior of SNF over the timescales associated with permanent disposal in a geologic repository. After leaving the spent fuel pool, SNF is transferred to dry storage in a dual-purpose canister (DPC). DPCs are considered “dual purpose” because they are designed for both storage and transportation, removing the need to transfer the fuel to a separate transportation cask. However, much research has been conducted investigating the feasibility of directly disposing of DPCs in geologic repositories. Direct disposal of DPCs could reduce worker exposure during repackaging, reducing the amount of low-level waste from the discarded DPCs and potentially saving billions of dollars. Therefore, direct disposal of as-loaded DPCs is desirable if it can be done safely.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Augmenting Traditional Static Analysis With Commonly Available Metadata

Developers and security analysts have been using static analysis for a long time to analyze programs for defects and vulnerabilities with some success. Generally a static analysis tool is run on the source code for a given program, flagging areas of code that need to be further inspected by a human analyst. These areas may be obvious bugs like potential bu er over flows, information leakage flaws, or the use of uninitialized variables. These tools tend to work fairly well - every year they find many important bugs. These tools are more impressive considering the fact that they only examine the source code, which may be very complex. Now consider the amount of data available that these tools do not analyze. There are many pieces of information that would prove invaluable for finding bugs in code, things such as a history of bug reports, a history of all changes to the code, information about committers, etc. By leveraging all this additional data, it is possible to nd more bugs with less user interaction, as well as track useful metrics such as number and type of defects injected by committer. This dissertation provides a method for leveraging development metadata to find bugs that would otherwise be difficult to find using standard static analysis tools. We showcase two case studies that demonstrate the ability to find 0day vulnerabilities in large and small software projects by finding new vulnerabilities in the cpython and Roundup open source projects.

97 MATHEMATICS AND COMPUTING↗