Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software library”

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.

1,599 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A cross-platform execution engine for the quantum intermediate representation

Hybrid languages like the quantum intermediate representation (QIR) are essential for programming systems that mix quantum and conventional computing models, while execution of these programs is often deferred to a system-specific implementation. Here, we develop the QIR Execution Engine (QIR-EE) for parsing, interpreting, and executing QIR across multiple hardware platforms. QIR-EE uses LLVM to execute hybrid instructions specifying quantum programs and, by design, presents extension points that support customized runtime and hardware environments. We demonstrate an implementation that uses the XACC quantum hardware-accelerator library to dispatch prototypical quantum programs on different commercial quantum platforms and numerical simulators, and we validate execution of QIR-EE on IonQ, Quantinuum, and IBM hardware. Our results highlight the efficiency of hybrid executable architectures for handling mixed instructions, managing mixed data, and integrating with quantum computing frameworks to realize cross-platform execution.

LLVM↗
DOE OSTI2022

Precursor Analysis Report: Cyber Attack on Thyssenkrupp Blast Furnace 2014

The Cyber Attack on Thyssenkrupp Blast Furnace 2014 Precursor Analysis Report leverages publicly available information about the Thyssenkrupp Steel Mill cyber attack and catalogs anomalous observables for each technique employed in the attack. This analysis is based upon the methodology of the Cybersecurity for the Operational Technology Environment (CyOTE) program. In December 2014, the German Government’s Federal Office for Information Security (BSI) released a report detailing a cyber attack on a German steel mill that occurred earlier that year, though exact dates and details of the attack were not revealed. While the report did not specify the name of the company, multiple sources identified the victim as one of Europe’s largest steel manufacturers, Thyssenkrupp AG. Further, Thyssenkrupp announced on 16 May of that year that Europe’s largest blast furnace, “Schwelgern 2,” located at its facility in Duisburg, Germany, would be offline for several weeks for repairs and upgrades, suggesting Schwelgern 2 was likely the target of the attack. The attack began in early 2014, when adversaries infiltrated the victim steel mill’s Information Technology (IT) network via a spearphishing campaign, then worked their way into the Operational Technology (OT) environment, where they executed software that caused denial of service, denial of control, and eventually a loss of control. This led to the blast furnace shutting down without proper safety procedures, resulting in catastrophic physical damage. No lives were lost in the incident, but ThyssenKrupp suffered $4 million in damage to the blast furnace and an additional $6 million in lost revenue. The adversaries required specialized knowledge and expertise in steel production, which enabled them to compromise a variety of internal systems and components across both IT and OT networks. The attack also demonstrated detailed knowledge of the industrial control systems (ICS) and production processes being used. This combination resulted in one of the earliest known publicly reported cybersecurity incidents resulting in physical damage to ICS equipment. Researchers and analysts identified 19 unique techniques (used in a sequence of 20 steps) utilized during the attack with a total of 454 observables using MITRE ATT&CK® for Industrial Control Systems. The CyOTE program assesses observables accompanying techniques used prior to the triggering event to identify opportunities to detect malicious activity. If observables accompanying the attack techniques are perceived and investigated prior to the triggering event, earlier comprehension of malicious activity can take place. Fifteen of the identified techniques used during the Thyssenkrupp cyber attack were precursors to the triggering event. Analysis identified 369 observables associated with these precursor techniques, 316 of which were assessed to have an increased likelihood of being perceived in the 120 days preceding the triggering event. The response and comprehension time could have been reduced if the observables had been identified earlier. The information gathered in this report contributes to a library of observables tied to a repository of artifacts, data sources, and technique detection references for practitioners and developers to support the comprehension of indicators of attack. Asset owners and operators can use these products if they experience similar observables or to prepare for comparable scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

TU-AB-BRC-08: Egs-brachy, a Fast and Versatile Monte Carlo Code for Brachytherapy Applications

Purpose: To introduce egs-brachy, a new, fast, and versatile Monte Carlo code for brachytherapy applications. Methods: egs-brachy is an EGSnrc user-code based on the EGSnrc C++ class library (egs++). Complex phantom, applicator, and source model geometries are built using the egs++ geometry module. egs-brachy uses a tracklength estimator to score collision kerma in voxels. Interaction, spectrum, energy fluence, and phase space scoring are also implemented. Phase space sources and particle recycling may be used to improve simulation efficiency. HDR treatments (e.g. stepping source through dwell positions) can be simulated. Standard brachytherapy seeds, as well as electron and miniature x-ray tube sources are fully modelled. Variance reduction techniques for electron source simulations are implemented (Bremsstrahlung cross section enhancement, uniform Bremsstrahlung splitting, and Russian Roulette). TG-43 parameters of seeds are computed and compared to published values. Example simulations of various treatments are carried out on a single 2.5 GHz Intel Xeon E5-2680 v3 processor core. Results: TG-43 parameters calculated with egs-brachy show excellent agreement with published values. Using a phase space source, 2% average statistical uncertainty in the PTV ((2mm){sup 3} voxels) can be achieved in 10 s for 100 {sup 125}I or {sup 103}Pd seeds in a 36.2 cm{sup 3} prostate PTV, 31 s for 64 {sup 103}Pd seeds in a 64 cm{sup 3} breast PTV, and 56 s for a miniature x-ray tube in a 27 cm{sup 3} breast PTV. Comparable uncertainty is reached in 12 s in a (1 mm){sup 3} water voxel 5 mm away from a COMS 16mm eye plaque with 13 {sup 103}Pd seeds. Conclusion: The accuracy of egs-brachy has been demonstrated through benchmarking calculations. Calculation times are sufficiently fast to allow full MC simulations for routine treatment planning for diverse brachytherapy treatments (LDR, HDR, miniature x-ray tube). egs-brachy will be available as free and open-source software to the medical physics research community. This work is partially funded by the Canada Research Chairs program, the Natural Sciences and Engineering Research Council of Canada, and the Ontario Ministry of Research and Innovation (Ontario Early Researcher Award).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

From PeleC to PeleACC, to PeleC++

PeleC is an Exascale Computing Project application for simulating compressible combustion in complex geometries. It has been built on top of the popular AMReX library. In the beginning of the Exascale Computing Project, PeleC was focused on KNL. It uses a mixture of C++, C, and kernels written in Fortran to obtain performance by focusing on vectorization. Recently we have taken two approaches in deciding PeleC's future for obtaining performance on exascale GPU machines. In the first programming model, we decorated the Fortran kernels with OpenACC directives. This expedited our ability to run at large scales on Summit's GPUs, where we achieved a significant speedup over the CPUs on Summit. The second programming model involved rewriting the Fortran kernels in C++ and using AMReX's Kokkos-like lambda abstractions for running on the GPU. This resulted in similar speedups on Summit's GPUs over merely utilizing the CPUs. Both approaches involved AMReX's management of memory transfers between the device and host. In this work, we compare and contrast the benefits and pitfalls to both programming approaches regarding performance, performance portability, and productivity. We also discuss advantages we have found in taking the time to modernize our code and why have chosen a specific pathway to prepare our code for the future DOE exascale machines.

exascale computing↗
DOE OSTI2021

Mid Infra-Red Laser Sensor for Continuous Sulfur Trioxide Monitoring to Improve Coal-Fired Power Plant Performance during Flexible Operations

During the course of this project, we performed exhaustive research and development of SO3/H2SO4 sensing technology for coal-fired power plant applications (Figure 1). The development culminated in a successful field campaign of a prototype continuous real-time H2SO4 monitor at a coal-fired power plant (TRL 6) accomplishing the primary goal of the project. The developed sensors utilize tunable laser absorption spectroscopy (TLAS) operating in the mid-infrared (Mid-IR) wavelength region, which is the so-called “molecular fingerprint” region. Systems operating in the Mid-IR have orders of magnitude more sensitivity than systems operating at shorter wavelengths, such as near-infrared (NIR). However, NIR systems are more widespread due to more mature supporting technology (e.g., fiber optics, optical components, etc.). In this project, we not only produced a specific Mid-IR sensor, we also advanced Mid-IR sensor technology in general through the development and demonstration of such supporting technology. In this project, we also developed proprietary broad tuning lasers enabling the ability to effectively measure SO3, H2SO4, H2O, and SO2. Different molecular species have unique spectral signatures that can be probed with lasers operating at different wavelengths. Standard TLAS uses relatively narrow wavelength tuning distributed feedback (DFB) lasers, which can typically only target a single species with narrow features, and are not appropriate for species with broad features, such as SO3 or H2SO4. In contrast, by developing unique, broad-tuning laser technology, we were able to measure these species, as well as SO2 and H2O simultaneously. Furthermore, to enable real-time analysis at a power plant, we modified a commercially available heated gas cell to operate in the Mid-IR wavelength range and fiber coupled the lasers to enable remote delivery of the beams. To generate reference data (library spectra), our collaborators at the University of California Irvine (UCI) developed a catalytic SO3 generation facility. It is worth mentioning that representative H2SO4 and SO3 Mid-IR spectra are not a part of any publicly available database and the data generated under this project is a valuable resource in and of itself. In addition, based on the UCI study we determined that detection of SO3 is complicated by the very strong SO2 absorption. For that reason, we concentrated on H2SO4 detection. Since SO3 and H2SO4 exist in a flue gas in a state of equilibrium, which depends on temperature and humidity, by measuring water concentration and controlling the temperature of the gas cell, we developed an approach to determine SO3 concentration from the H2SO4 measurement. During the development phase of the project, we performed three testing campaigns at our collaborator’s FERCo flue gas facility with conditions representative of the coal-fired power plant (~ 40ppm SO3, 1700ppm to 2800 ppm SO2, 10% water) with the exception of particulate matter. After three test campaigns at FERCo we performed field testing at Harrison Power Station. The final system was mounted on a duct and measured H2SO4, SO2 and water. The tests were highly successful with a demonstrated real-time H2SO4 precision of 1 ppm with a 1 second update. Our collaborators at EPRI conducted an industry survey and determined that there is a very high interest for the SO3/H2SO4 monitoring in the power generation industry as well as in heavy industries in general. Furthermore, work performed by OptoKnowledge beyond the scope of this project under a synergistic DOE SBIR determined another approach to SO3 detection. We applied for Phase II on this SBIR for development a of versatile SO3/H2SO4 sensor but were not selected. We are currently looking for another opportunity to leverage all the technological advancements produced by this project including but not limited to the flue gas facility at UCI, the hardware and software developed, and relationships with FERCo, EPRI CEMTEK, and Harrison Station.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods like the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite like GAMESS is maintained and developed is considered. Future developments are briefly summarized.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advancements in Multiphysics Microdepletion Analysis of an eVinci TM -like Microreactor Leveraging OpenMC-CRAB Workflow

Nuclear microreactors (MRs) are a class of nuclear reactor technology, characterized by reduced dimensions, modular design, and reduced power output in contrast to conventional Light Water Reactors (LWRs). MRs are proposed for supplying electricity and eventual process heat to remote locations, such as military installations and disaster-affected areas. Current research work sponsored by the US Department of Energy Microreactor Program (MRP) is devoted to the development of novel modeling and simulation tools to better support MR vendors and regulatory bodies. Notably, the NRC is projected to utilize the CRAB multiphysics software driver for executing both design and beyond-design-basis accident analyses. Furthermore, the NRC has been utilizing the MELCOR code to calculate mechanistic source terms during accidents. Since MELCOR relies on isotopic inventory and reactor temperature/power profiles under accident conditions, which theoretically can be derived from CRAB, the goal is to establish a comprehensive CRAB-MELCOR computational framework. Past work was focused on testing and demonstrating CRAB's capability to generate results that can be used to inform mechanistic source term calculations in MELCOR. In particular, a computational workflow leveraging OpenMC-generated microscopic cross sections and CRAB was first applied to perform multiphysics microscopic depletion calculation followed by an accident scenario for a stylized microreactor problem. In fiscal year 2024, the research work has been focused on applying the OpenMC-CRAB workflow, which was first tested in fiscal year 2023, to a realistic 3D heat-pipe cooled MR problem representative of the eVinci TM design. The latter computational problem was developed with inputs from WEC to conserve selected neutronic and thermal characteristics of the eVinci TM design without releasing proprietary data. The results of this simulation, encompassing isotopic inventory, power density distribution, and kinetic parameters, will inform both MELCOR and the WEC-developed FATE code for mechanistic source terms calculations. The results from the two codes will then be compared for code verification purposes. This report contains the design characteristics of the realist heat pipe cooled microreactor developed as a use-case for the verification exercise, and the current results for the multiphysics microscopic depletion performed with the OpenMC-CRAB workflow. The results include eigenvalue as a function of time, power distribution at EOL, in addition to nuclides inventory's time evolution and spatial distribution. Finally, we report improvements to the workflow efficiency achieved through a collaboration with the NEAMS programs. Through this collaborative effort, we were able to strongly decrease the computational time for the multiphysics microdepletion calculation (i.e., from 17.4 hours to 5.7 hours on 280 processors) in addition to simplifying the interface to generate isotopics spatial distribution utilizable by FATE and MELCOR. Future work, including the improvement of the current microscopic cross-sections' library and the simulation of an accident scenario at EOL, is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A Data Processing Pipeline To Extract A Knowledge Graph From Heterogeneous Data For Socio-technical Analysis Of Critical Infrastructure Influence

The code is written in Python and consists of the following pipeline that is implemented in Apache Airflow. This pipeline intends to understand the companies that are directly or indirectly involved with a type of critical infrastructure system at some point in that system's lifecycle. The pipeline takes a configuration file that specifies a list of initial companies to consider, a geographic region of interest, and a set of SEC form types as well as other data sources (e.g. CrunchBase) from which to extract entities and relations. There are four main components to this pipeline as currently implemented: Entity Extraction, Network Construction, Analysis, and Visualization. First, Entity Extraction, is implemented as the `topear-extract_organizations` Apache Airflow workflow. Given an initial query that specifies a geographic region of interest and a time interval, the software will extract CI facilities of interest and organizations that have a direct influence relationship to those facilities (e.g. ownership). During the course of the LDRD, we focused on Electric Vehicle charging stations and this information is available via the Department of Energy (DOE) database on fueling stations maintained by NREL. Within the context of the DOE CESER project, we have focused on Battery Energy Storage Systems (BESS). Second, the Network Extraction component will iteratively construct a social network graph given the set of organizations and people extracted in the previous step. Organizations (and eventually People if desired) are then fed as a query to the `topgear-construct_social_network` Apache Airflow workflow which given a set of initial companies and data sets (e.g. SEC EDGAR form types, OpenCorporates, Crunchbase). This Airflow workflow will iteratively query such data sources to discover relationships with new organizations and people. For example, this module can iteratively query SEC EDGAR for metadata that documents the number of each type of form for the given set of companies and their location. This forms metadata represents a catalog of data sources from SEC EDGAR for the extracted social network knowledge graph. The pipeline then downloads these forms from the website and saves them in a build directory for further processing. These documents are then parsed for entities and relations. Again, we note that in additional to SEC data sources, this step can also pull in information on organizations via API services such as CrunchBase and OpenCorporates or bulk data sources. At the end of this step, the resultant social network, the Critical Infrastructure network, and the edges that encode relationships between organizations and CI facilities, form the Adversarial Socio-Technical Network (ASTN) that informs the analysis. Third, the Analysis component processes these generated ASTN. Previously, that has included the ability to compare prevalence of different vendors for a given infrastructure component type across different regions as well as identify common public and private investors across those vendors. This was demonstrated for EV Charging Stations across several different metropolitan areas within an IEEE PES GridEdge publication. More recently, we have looked at ways to identify infrastructure owners and operators of BESS with the most nameplate capacity across different states as well as other indictors of risk resulting from changes in ownership over time. Finally, the Visualization component consists of an HTML/CSS/JS framework by which users can interact geospatial, operational, and organizational relationships across a given portfolio of Critical Infrastructure facilities. The objective is to provide a library of UI/UX modules that can be repurposed for stakeholder-specific dashboards. All of the modules are related via a common event model that enables UI actions in one view to percolate across the other views.

Weaver, Gabriel [Idaho National Laboratory (INL), ↗
DOE OSTI2020

Revisiting Five Years of CASMI Contests with EPA Identification Tools

Software applications for high resolution mass spectrometry (HRMS)-based non-targeted analysis (NTA) continue to enhance chemical identification capabilities. Given the variety of available applications, determining the most fit-for-purpose tools and workflows can be difficult. The Critical Assessment of Small Molecule Identification (CASMI) contests were initiated in 2012 to provide a means to evaluate compound identification tools on a standardized set of blinded tandem mass spectrometry (MS/MS) data. Five CASMI contests have resulted in recommendations, publications, and invaluable datasets for practitioners of HRMS-based screening studies. The US Environmental Protection Agency’s (EPA) CompTox Chemicals Dashboard is now recognized as a valuable resource for compound identification in NTA studies. However, this application was too new and immature in functionality to participate in the five previous CASMI contests. In this work, we performed compound identification on all five CASMI contest datasets using Dashboard tools and data in order to critically evaluate Dashboard performance relative to that of other applications. CASMI data was accessed via the CASMI webpage and processed for use in our spectral matching and identification workflow. Relative to applications used by former contest participants, our tools, data, and workflow performed well, placing more challenge compounds in the top five of ranked candidates than did the winners of three contest years and tying in a fourth. In addition, we conducted an in-depth review of the CASMI structure sets and made these reviewed sets available via the Dashboard. Our results suggest that Dashboard data and tools would enhance chemical identification capabilities for practitioners of HRMS-based NTA.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Minimizing User Burden in Building Energy Analysis

Infrared (IR) imagery can provide valuable data to building energy auditors and the developers of Building Energy Models (BEM). IR by itself is not a panacea to the challenges facing energy auditors and modelers, but when combined with other techniques and Open Studio software, IR can play an important role in the minimization of the burden of creating building energy models and the diagnosing of building energy problems. While IR can be used to detect and quantify many aspects of a building’s thermal and energy performance, in Phase I we focused on a single capability: determining the R value of exterior walls. Two key attributes of using IR methods for R value determination is that IR methods measure the actual, installed performance in-situ – as opposed to relying on manufacturer or handbook values – and they can integrate performance over a wide area and thus can accurately account for thermal bridging and spatial anomalies such as “spotty” insulation installation. There is a trade between user burden and accuracy when using IR to determine a wall’s R value. The least burdensome approach is computing R values from a single IR snapshot. The most burdensome method we considered was to combine IR images with thermocouple and heat flux sensor (HFS) measurements made over the course of several days. The more effort users expend in making the measurements, the more accurate and robust the method becomes. Through a series of field tests at the ThermoAnalytics Inc. (TAI) HQ at Calumet, MI and a simulation of a 1980s era office building in Flagstaff, AZ we investigated methods to compute the R value of exterior walls using data from IR cameras, thermocouples and heat flux sensors. To compute R values the heat flux through the wall, along with temperatures on either side of the wall, must be known, modeled, or estimated – and the temperature difference through the wall should be at least 10 C. If only outdoor IR images are used in the method, the inside ambient air temperature must be estimated, and a simulation model created to predict heat fluxes and the indoor wall temperature. The accuracy of the resulting R value depends on the accuracy of the simulation and the assumptions it employs. Applying thermocouple sensors, either as iButtons or as part of a wireless sensor network, to the walls allows for the acquisition of temperatures over time. The longer the duration of the measurements, the more accurate and robust the method will be. Even with this level of measurement effort, modeling is still required to calculate the heat flux through the wall. The addition of a heat flux sensor alleviates the need for any modeling, and that results in improved accuracy. Driving the accuracy of the methods are several factors that can be quantified. How well are the wall surface temperatures known? Modeling the heat flux through walls requires knowing the details of the surrounding radiant environment and the heat transfer due to convection. How well are these boundary conditions known? Even when the heat flux is measured, error can arise when some of the measured heat is not conducted through the wall but rather goes into the thermal mass of the wall. Thermal mass effects occur when conditions and the heat flux values are fluctuating. A technique that is not robust will suffer in accuracy when these fluctuations occur. The R value determination method must be incorporated into the Open Studio process. In our proposed process, IR images and other measured sensor data are input into a module that would send R value estimates to a building library developed by team member Skidmore, Owings, & Merrill (SOM). This library of parameterized building models would use these R values along with information about the building’s age, type, and location to query the Open Studio BCL to identify the most likely constructions for walls, roof, and windows. Being based on the IMF feature of Energyplus, the SOM building library can quickly produce Energyplus models of buildings based on a limited number of user inputs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Status on the Development of a Modeling and Simulation Framework for the Economic Assessment of Nuclear Hybrid Energy Systems

An effort to design and build a modeling and simulation framework to assess the economic viability of Nuclear Hybrid Energy Systems (NHES) was undertaken in fiscal year 2015 (FY15). The purpose of this report is to document the various tasks associated with the development of such a framework and to provide a status on its progress. Several tasks have been accomplished. First, starting from a simulation strategy, a rigorous mathematical formulation has been achieved in which the economic optimization of a Nuclear Hybrid Energy System is presented as a constrained robust (under uncertainty) optimization problem. Some possible algorithms for the solution of the optimization problem are presented. A variation of the Simultaneous Perturbation Stochastic Approximation algorithm has been implemented in RAVEN and preliminary tests have been performed. The development of the software infrastructure to support the simulation of the whole NHES has also moved forward. The coupling between RAVEN and an implementation of the Modelica language (OpenModelica) has been implemented, migrated under several operating systems and tested using an adapted model of a desalination plant. In particular, this exercise was focused on testing the coupling of the different code systems; testing parallel, computationally expensive simulations on the INL cluster; and providing a proof of concept for the possibility of using surrogate models to represent the different NHES subsystems. Another important step was the porting of the RAVEN code under the Windows™ operating system. This accomplishment makes RAVEN compatible with the development environment that is being used for dynamic simulation of NHES components. A very simplified model of a NHES on the electric market has been built in RAVEN to confirm expectations on the analysis capability of RAVEN to provide insight into system economics and to test the capability of RAVEN to identify limit surfaces even for stochastic constraints. This capability will be needed in the future to enforce the stochastic constraints on the electric demand coverage from the NHES. The development team gained experience with many of the tools that are currently envisioned for use in the economic analysis of NHES and completed several important steps. Given the complexity of the project, preference has been given to a structural approach in which several independent efforts have been used to build the cornerstone of the simulation framework. While this is good approach in establishing such a complex framework, it may delay reaching more complete results on the performance of analyzed system configurations. The integration of the previously reported exergy analysis approach was initially proposed as part of this milestone. However, in reality, the exergy-based apportioning of cost will take place only in a second stage of the implementation since it will be used to properly allocate cost among the different NHES subsystems. Therefore, exergy does not appear at the level of the main drivers in the analysis framework; the latter development of the base framework is the focus of this report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

Advanced SOFC quality control and the role of manufacturing defects on stack reliability

Atrex Energy rolled out first commercial SOFC generator products in 2013. Following that we embarked on implementation of mass production in order to ramp up productivity and lower the costs in the meantime. To ensure the goal of mass production, an accurate, cost effective and consistent quality control methodology needs to be in place. The aim of this DOE funded project is to devise a machine and methodology for cell quality control using modern non-destructive testing methods and cell rejection criterion based on the empirical determination of the role of manufacturing flaws on stack degradation and reliability. The work scope encompasses two goals aimed at 1) developing automatable techniques for objective flaw detection and 2) developing an understanding of the ramifications of flaws on stack and cell reliability. By accomplishing these two goals it is likely that an advanced stage Quality Control (QC) device can be developed. This project scope encompasses: 1) Invention of a high speed automated instrument for spatial mapping of surface and sub-surface imperfections in production cells (ex situ) using modern imaging techniques. 2) Identification of rejection criteria upon which to screen cells from production by: a. Magnification of the effect of production defects on cell degradation through fabrication and testing of micro-cells which encompass these imperfections. b. The long term testing of full size production cells both as single cells and in small stacks, in order to correlate the (in situ) performance and degradation with the imperfection type (identified ex-situ) and hence be able to classify imperfections as actual defects. We first built a defect library with the cells containing imperfections classified and documented. The cells with various imperfections were tested to evaluate the imperfection potential impact on cell electrochemical performance. Evaluations were conducted on minicell, single cell and stack levels. On minicell level, pinhole, scratch, crack, contamination, and pop-out were investigated. We found: (1) none of the imperfections causes calamity type failure, (2) pinholes do not show measurable effect on minicell performance, (3) imperfection such as electrolyte abrasion, scratch, and pop out cause minicell potential to fluctuate by different degrees; the cells with such imperfections exhibit sensitivity of potential to high fuel utilization and “self healing” behavior. On single cell level we evaluated the effect of long crack and high population pinholes. The cell with a long crack does not show performance compromise at start. There is no potential fluctuation or sensitivity to high fuel utilization or catastrophic failure. The voltage degradation due to thermal cycle is slightly higher than normal cells. The cells with high population pinholes perform surprisingly as well as normal cells over a period of 3,500h. Even when subjected to cyclic load condition, the cells 4 show negligible degradation. Our cells appear to have “self immunization” feature which fends off the adversary effects of certain imperfections. Defects were also evaluated on stack level. Two stacks were made each with ten defect cells. The first stack was tested for 13,2000h under steady load condition while the second was tested for 9,000h under transient load conditions. We found: (1) voltage fluctuation and sensitivity to thermal cycle were not observed in contrast to minicell test, (2) of the various cracks examined, including longitudinal cracks, circumferential crack, crack near closed and open ends, none exhibited distinguishable impact on performance, (3) defect cells exhibited higher degradation rates on average under steady load condition; Individually some defect cells performed similarly or better than reference cells, (4) under steady load condition all the cells degrade monotonically with degradation rate increasing with time, (5) under cyclic (transient) load condition, both defect and reference cells exhibit complex dynamic behaviors, (6) fuel loss (leak) was measured occasionally for stack 1, the cells with crack showed the highest fuel loss increase, nevertheless, no correlation with performance degradation could be established with certainty, and (7) no catastrophic failure ever happened other than incremental degradation. Posttest inspection and microanalysis were carried out; the observations are summarized in what follows: (1) Electrolyte pinholes, some scratches and IC porosity became undetectable after testing. (2) The cathode covering the popout defects is intact without any cracking or delamination. (3) The scratch was covered by cathode in pristine state; Ni oxide or NiO-YSZ filled the space where electrolyte was missing. After >10,000h operation, the cathode that used to cover the scratch was missing but the size of the scratch remains unchanged; there is a fairly dense Ni oxide layer on top of the scratch. (4) All of the cracks remain their original dimensions. Current collector kept longitudinal cracks closed which mitigated the adversary effect of crack on cell performance/durability. (5) The crack at open end showed some deterioration manifested by a small oxidation zone in the vicinity of the crack. Wide spectrum of imaging techniques were surveyed and evaluated to detect tubular SOFC defects. Optical methods, with and without assistance of expression technique, were applied to detect all kinds of defects. Thermograph method in combination of CO2 leak, ultra-sound and electric potential was also extensively evaluated to detect surface and under-surface defects such as short circuit. Based on these works, laser camera has been selected to build the final QC machine vision inspection system. The Gocator laser profile sensor was setup to generate 2D surface intensity images of tubular SOFC surface which were analyzed with Sherlock machine vision software supplied by Teledyne Dalsa. The Gocator output was integrated with the Sherlock software and enabled live images to be imported by a Common Vision Blox driver. The Sherlock software was programmed to detect various surface defects from the surface intensity images generated by the Gocator. These measurements are subsequently used to make pass or fail decisions.

01 COAL, LIGNITE, AND PEAT↗