Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “IMPLEMENTATION”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART)

This report contains key findings from a project titled Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART), which was carried out through a collaborative effort of a team of researchers from Texas A&M Engineering Experiment Station, Temple University, and Quanta Technology, LLC. The in-kind support came from OSIsoft (acquired by AVEVA), which provided their PI Historian software to demonstrate the use case of streaming PMU data. The first section of the report describes the project goals and objectives related to the development of Machine Learning (ML) models capable of detecting and classifying events by processing phasor measurements captured in the field by Phasor Measurement Units (PMUs). The data for this study was contributed by the utilities/ISOs from the Western and Eastern interconnects and ERCOT, further referred to as Interconnect B (IC B), Interconnect A (IC A), and Interconnect C (IC C), respectively. The approach that the BDSMART Research Team proposed and the key research tasks defined by the team are outlined in this section. The next section describes the technical approach. We first discuss the data constraints related to the PMU measurements and data interpretation constraints imposed by the data contributors. They provided neither the topological information of the grid nor PMU placement locations and captured recorded data at very few locations in the system with the reporting rate of either 30 or 60 fps. The recordings are mostly positive sequence voltage, frequency, and ROCOF, and in some limited cases, three-phase voltages and currents. We then reflect on the bad data issues that stem from poor recording practices and vague definitions of the PMU status bits to supposedly be used for bad data identification. Finally, the data discovery points to imprecise time stamps with incomplete event start/end time, as well as inconsistent and incomplete event labeling, which combined make the implementation of the data models using supervising learning quite challenging. Following the data discovery study, we hypothesize that because the IC B data has the most complete labels, we should focus our model development on that data and then test it on data from other interconnects. We also define the common metrics used to evaluate the results from the ML algorithm tests. We concluded this section by summarizing the common ML models we used and explaining how we implemented and tested them. The issues from this section are expanded in the Training Dataset Report from this project. The final section of this report deals with the accomplishments and conclusions. As the accomplishments, we formulate the problem we are solving and what is achieved by solving the problem. We then reflect on each of the analytics tools we developed and point out the performance of each tool when applied to solving the mentioned problems. We reference this work for further details to the papers we published on each tool. In the conclusions, we give recommendations on how to improve future PMU recording practices to facilitate the ML algorithm implementation and guidance for the future standardization work aimed at clarifying the ambiguities associated with the PMU status bits. We finally list future tasks that can bring about further improvements in the proposed algorithms. The issues from this section are expanded in the Training, and Test Dataset Report filed at the project completion date.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Oak Ridge National Laboratory Pilot Demonstration of an Attestation and Anomaly Detection Framework using Distributed Ledger Technology for Power Grid Infrastructure

This report summarizes the design and pilot demonstration of a framework called Grid Guard that was created to provide increased data and device trustworthiness to electric grid devices by leveraging distributed ledger technology (DLT), specifically blockchain. Grid Guard contains a combination of core cryptographic methods such as the secure hash algorithm (SHA), and asymmetric cryptography, private permissioned blockchain, baselining configuration data, consensus algorithm (Raft) and the Hyperledger Fabric (HLF) framework. The system implements a low energy, fast, and robust enhancement to system trustworthiness within and across electric grid systems such as substations, control centers and metering infrastructures. Blockchain is a distributed database structured that provides a practically unalterable (immutable) timeline of stored transactions. By relying on hashing and the Raft consensus algorithm, if an entity tries to illegitimately alter a record at one instance of the database the other ledger nodes are not altered. They work to cross-reference each other and easily locate any incorrectly added data and remove it. The bulk raw data is stored in an off-chain storage (outside of the blockchain ledger) and a hash of this baseline data is stored in the Blockchain ledger via hashing windows of time-series and configuration data, after aggregation and filtering. The bulk off-chain data repository is then considered to be trust-anchored using the hashes stored in the blockchain. To secure the electric grid testbed devices and data, device configuration baselines were compared to those baselines that had been previously stored in the ledger. Statistical baselines for device configurations, network communication patterns, and high-speed sensor data are calculated and then stored off-chain and hashes stored in the ledger. Measurements such as three-phase voltage and current, frequency, breaker status, protection scheme settings, network configuration settings (and other device configuration artifacts) and network traffic features (packet interarrival times) are compared every minute or other selected time windows. During phase 1 of the Grid Guard DLT project different DLT technologies were studies, and an assessment was performed on DLT technology vulnerabilities, uses, and key characteristics. DLT consensus protocols were studies (e.g., RAFT, named after Reliable, Replicated, Redundant, And Fault-Tolerant). Also, cryptography, public, private and permissioned or permissionless systems were assessed. Grid Guard implements a permissioned private DLT. Consensus algorithm selection and choice of DLT implementation depended heavily on the use-case. For this use-case, parameters were selected to measure performance and existing tools for assessment. Benchmarking was performed theoretically and practically. During phase 2 hashed transactions/blocks were inserted into the ledger every second. During phase 2 of the Grid Guard DLT project, a prototype framework was developed and demonstrated for attestation of critical substation devices and data using precision timing systems that use PTP and IRIG-B protocols) on a testbed of operational devices that emulated a distribution substation, control center, and power metering infrastructure using real Operational Technology (OT). The testbed includes OT devices such as protective relays, human machine interfaces (HMI), and power meters. To determine when to collect and compare system and network baselines, an initial examination of an anomaly detection capability to identify malicious manipulation of data streams was conducted. The resulting anomaly detection was demonstrated in a set of experiments and leveraged to trigger device artifact attestation checks. Attestation checks occur against device configuration baselines when compared with the immutable blockchain-stored baselines, which provided a cryptographically supported means by which to store baselines. The electrical substation-grid testbed was created to test the Grid Guard framework. The testbed emulates the operations of a portion of a power grid and SCADA systems as closely as possible. The testbed integrates real protocols, mainly IEC 61850 standard protocols, such as the Sampled Value (SV) and the GOOSE protocols. The testbed also supports DNP3 and other layer 2 and layer 3 protocols such as Telnet, SSH, SFTP/FTP and other proprietary protocols needed to connect to industrial control system equipment. The testbed emulates real power conditions using the OpalRT hardware-in-the-loop (HIL) device which can create fault situations that cannot be easily tested on real systems. The electrical substation-grid testbed was created using real measurement, communication, and protection devices that electrical utilities commonly use.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Spiner vs EOSPAC6: capabilities, performance, and accuracy considerations

The Cross-Cutting Capabilities Project (XCAP) has significant interest in collecting and implementing a common set of libraries to be used across the various production hydro-codes at LANL. A common set of libraries will help facilitate comparisons between codes. In addition to minimizing variables for comparisons and physics validation, efforts toward optimization and porting of packages for future architectures will be more efficient. As part of this push Singularity is being considered as an inclusive materials interface library. One crucial part of this library is the importing, inverting, and interpolation of equation of state (EOS) data. Currently the Lagrangian Applications Project (LAP) and Safety Applications Project (SAP) are using EOSPAC6 for delivery EOS data. The Eulerian Applications Project (EAP) has historically used TEOS and more recently an implementation of Singularly via a package called Spiner for EOS interpolation. Singularity currently has the option of using several analytical EOS models, directly employing SESAME via Spiner, or directly using SESAME via EOSPAC6 (not optimized yet). Assuming Singularity moves forward as a common platform for the implementation of materials models, a decision will need to be made on which EOS interpolation package XCAP should move forward with, given resources and people are finite. We will attempt to address pros and cons of Spiner and EOSPAC6 and the trade-offs that should be considered during the decision making process. This report is intended as an ASC-PEM-EOS perspective on what is needed in an EOS interpolation package. Performance and Accuracy sections will mostly address data from comparison studies in LA-UR-22-22699.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Tearing parameter failure integration with the multilevel solver.

The tearing parameter criterion and material softening failure method currently used in the multilinear elastic-plastic constitutive model was added as an option to modular failure capabilities. The modular failure implementation was integrated with the multilevel solver for multi-element simulations. Currently, this implementation is only available to the J 2 plasticity model due to the formulation of the material softening approach. The implementation compared well with multilinear elastic-plastic model results for a uniaxial tension test, a simple shear test, and a representative structural problem. Necessary generalizations of the failure method to extend it as a modular option for all plasticity models are highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coastal Acoustic Buoy for Offshore Wind: Project Synthesis

The population of North Atlantic right whales is critically endangered and their habitat overlaps with offshore windfarm leases. It is therefore imperative that effective mitigation strategies be used to avoid impacts on right whales during the construction of offshore windfarms. The Department of Energy issued FOA Number DE-FOA-0001924 to encourage the development of technology that could monitor large exclusion zones for right whales in order to mitigate potential impact of construction noise on right whales. This report summarizes the past two years of the development and evaluation of the Coastal Acoustic Buoy for Offshore Wind (CABOW) project which aimed to develop technology to monitor large exclusion zones for North Atlantic right whales. Over the course of the project SMRU Consulting have implemented a rigorous design process including comparison of different approaches (e.g., single sensor vs multiple sensors), as well as consideration of placement and timing of acoustic monitoring. We have evaluated critical components of the CABOW system including, reliability, detection range, and bearing accuracy in areas adjacent to offshore windfarm leases in Maryland by conducting 3,536 playbacks of simulated right whale upcalls. The maximum call detection range was 7.5 km when noise was 99 dB re 1µPa rms (50-225 Hz), but this reduced to < 1 km when ambient noise levels were high. Our detection probability in the field was measured as a function of range as well as the source-to-noise level ratio allowing us to build a model to predict the probability of detection under various scenarios (sample size: 3,536 calls x 5 buoys = 17,680). The median bearing error was -0.25° but this is likely an underestimate of error due to experimental design. Using the published recall and precision of the two detectors we implemented in the CABOW system, we estimate that at a recall of 80%, our precision was > 80%, within the range of what we were aiming for in this project. To estimate our exclusion zone false negative and false positive rates, we built a simulation model using the empirical data from our field trial. We modelled three to nine CABOW units placed on the 10 km exclusion zone and estimated our false negative rate to be 1% or less (which was our project goal) and our false positive rate to be between seven and nine percent, slightly above our goal of 5%. However, we also modelled an equivalent PAM system that does not have bearing capabilities and found the false positive rates for that system to be six to eight times higher than the CABOW rate. This higher false positive rate of PAM systems with low spatial information could have significant cost repercussions for offshore wind developers by adding work shutdowns or delays without providing additional protection for right whales. The model we built allows us to explore the placement of PAM systems under various scenarios and will thus help facilitate planning of PAM mitigation systems to meet NOAA Incidental Harassment Authorizations for specific windfarms. We achieved an average system uptime of 98.3%, just below our goal of 99%. The issues that caused these short losses of data have been identified and fixed. Right whale detections and audio clips were typically transferred via radio from the buoys to the base station in two to four seconds. We therefore believe we have developed a highly robust real-time PAM system. Based on the above, we feel we have achieved the stated funding goal of developing a cost-effective and robust real-time PAM system that enables the monitoring of large exclusion zones for right whales during the constructions of offshore windfarms. This should lead to decreased costs and risks for the offshore wind sector while providing robust mitigation for right whales. It is important to state that PAM mitigation will need to be implemented with other mitigation strategies (e.g., visual observers) to provide a complete mitigation strategy to ensure that any effects on right whales from offshore windfarm construction is minimized.

17 WIND ENERGY↗
DOE OSTI2022

Xyce™ Parallel Electronic Simulator Users' Guide (V.7.6)

This manual describes the use of the Xyce™ Parallel Electronic Simulator. Xyce™ has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. (2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. (3) Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using modern coding practices. Xyce™ is a parallel code in the most general sense of the phrase—a message passing parallel implementation—which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel eficiency is achieved as the number of processors grows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The VTK-m Users' Guide (V.1.9)

High-performance computing relies on ever finer threading. Advances in processor technology include ever greater numbers of cores, hyperthreading, accelerators with integrated blocks of cores, and special vectorized instructions, all of which require more software parallelism to achieve peak performance. Traditional visualization solutions cannot support this extreme level of concurrency. Extreme scale systems require a new programming model and a fundamental change in how we design algorithms. To address these issues we created VTK-m: the visualization toolkit for multi-/many-core architectures. VTK-m supports a number of algorithms and the ability to design further algorithms through a top-down design with an emphasis on extreme parallelism. VTK-m also provides support for finding and building links across topologies, making it possible to perform operations that determine manifold surfaces, interpolate generated values, and find adjacencies. Although VTK-m provides a simplified high-level interface for programming, its template-based code removes the overhead of abstraction. VTK-m simplifies the development of parallel scientific visualization algorithms by providing a framework of supporting functionality that allows developers to focus on visualization operations. Consider the listings in Figure 1.1 that compares the size of the implementation for the Marching Cubes algorithm in VTK-m with the equivalent reference implementation in the CUDA software development kit. Because VTK-m internally manages the parallel distribution of work and data, the VTK-m implementation is shorter and easier to maintain. Additionally, VTK-m provides data abstractions not provided by other libraries that make code written in VTK-m more versatile.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

ExaSGD: 2022 Kernel Thrust Activities

The Kernel Thrust milestone ADSE22-407 covers the development of device-capable optimization algorithms and solvers technologies required by the ExaSGD project’s software stack in order to solve security-constrained alternating current optimal power flow (SC-ACOPF) problems on emerging exascale architectures. To this extent, in FY22 the main objective of the Kernel Thrust was (i) provide sparse optimization solver that runs efficiently on hardware accelerator devices (i.e., NVIDIA and AMD GPUs) to perform intra-node computations, (ii) strengthen the reliability and increase the performance of the mixed-dense sparse (MDS) solver of HiOp for deployment on the FY22 target architectures, Summit and Crusher, and (iii) increase performance by improving the mathematical algorithm and refining the parallel MPI-based implementation of the coarse-grain parallel solver HiOp-PriDec for capabilities deployment on the FY22 target architectures, Summit and Crusher. This document presents the developments and contributions done by the Kernels Thrust Team in FY22 toward completion of the above-mentioned objectives. These contributions progressed along four main development (sub)thrusts: (1) Design and implementation of a sparse optimization solver for use on hardware accelerators; (2) Improvement of the mathematical algorithm and of the parallel implementation of HiOp-PriDec to ensure readiness and efficient coarse-grain parallelism for FY23 target exascale machine; and (3) Support Software and Application Development Thrusts of the exaSGD project in their deployment of the project’s software stack on AMD- and NVIDIA-based architectures. The development of the sparse optimization solver (thrust 1 above) was new in FY22 and resulted in a new sparse solver in HiOp (available as of version 0.6). The second development thrust was a continuation of the efforts from FY21 and improved the mathematical algorithm and the communication strategy of the HiOp-PriDec solver. The last developement thrust is a large collaborative effort. Namely, the project’s teams from multiple labs (LLNL, PNNL, ORNL, and NREL) performed large-scale demonstration of the ExaSGD software stack, namely the optimization solvers of HiOp interfaced with the modeling front-end ExaGO and the stochastic sampler PowerScenarios. These demonstration efforts solved large-scale instances of the SC-ACOPF challenge problem of medium network sizes (10, 000-bus system) and large number of contingencies on Summit (NVIDIA accelerators) and Crusher (AMD accelerators) systems at ORNL.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Initial Engineering and Design for CO 2 Capture from Ethanol Facilities

The Energy & Environmental Research Center, in partnership with the U.S. Department of Energy (DOE) and North Dakota ethanol producer Red Trail Energy, LLC (RTE) and with technical support from Trimeric Corporation and the KLJ engineering firm, completed an initial engineering design (IED) for a hybrid capture system to estimate associated costs for retrofitting the operational RTE ethanol plant. The project goal supports the DOE objective to achieve negative emissions using commercially available technologies to attain Technology Readiness Level 6+ from a facility emitting >100,000 tonnes of carbon dioxide (CO 2 ) annually. The RTE CCS (carbon capture and storage) Project is currently operating a CO 2 capture facility, adjacent to the RTE ethanol facility in western North Dakota, and injecting the CO 2 more than a mile below RTE property for permanent storage. This novel hybrid capture system would process about 310,000 tonnes of CO 2 annually for the RTE host site and includes CO 2 capture and compression from bioprocessing as well as capturing CO 2 produced from natural gas boilers. The bioprocessing capture process, based on past RTE studies investigating ammonia liquefaction technology with CO 2 produced from natural gas boilers captured utilizing chemical absorption with amine, was shown to be technically feasible for the RTE site. Activities conducted to support the project goal included 1) designing a hybrid capture system using CO 2 emissions from both bioprocessing and heat production at the host site facility; 2) conducting a pre-front-end engineering and design analysis of the hybrid capture system to include environmental health and safety, a constructability report, identification of permits, and corporate approvals; and 3) performing a techno-economic assessment in accordance with DOE’s methodology, as demonstrated by the bituminous baseline study. In addition, a life cycle assessment (LCA) was completed for the RTE site to estimate the carbon life cycle for ethanol-CCS implementation using the designed hybrid capture system to determine potential for net-negative carbon emissions. Results of the IED showed the hybrid system to be technically viable, with a moderate estimated cost of $55/tonne CO 2 captured for the hybrid system. The cradle-to-gate LCA showed preliminary net-negative carbon emissions potential anticipated from implementing a hybrid CCS system at a commercial scale. Recommended next steps toward potential implementation include hybrid capture system demonstration testing for detailed engineering and LCA model comparisons with low-carbon fuel incentive programs for financial support.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Nanosecond Gated CMOS Camera (NSGCC) ICD (Rev. 2.1)

The Ultra-Fast X-ray Imager (UXI) program is an ongoing effort at Sandia National Laboratories to create high speed, multi-frame, time-gated Read Out Integrated Circuits (ROICs), and a corresponding suite of photodetectors to image a wide variety of High Energy Density (HED) physics experiments on both Sandia’s Z-Machine and LLNL’s National Ignition Facility (NIF). Several cameras have been designed over the length of the program; one of the most recent is the Icarus, which is an improvement on past imagers (Furi and Hippogriff). A second sensor that can be connected is the Daedalus sensor. The Icarus is a 1024 × 512-pixel array with either 25 μm or 8 µm spatial resolution containing four frames of storage per pixel and has improved timing generation and distribution components while achieved 2 ns time gating. The Daedalus sensor is also a 1024 x 512-pixel array with 25 µm special resolution containing three frames of storage per pixel and has an increased set of features for a wider variety of applications from interlacing of rows in each frame to configurability of all shutters. See Section 14 for details regarding the Icarus implementation of the firmware and Section 15 for details regarding the Daedalus implementation. Due to the unique test environments UXI sensors are targeted for, full custom hardware was required to physically mount an Icarus or Daedalus sensor, manage its various functions, and read out pixel data for transfer to a host computer. Beyond experimental functionality, the hardware also needed to accommodate sensor characterization requirements. Lawrence Livermore National Laboratory’s ‘Version 4.0 Board’ was the result of these efforts. It mounts all the components required to fully utilize the Icarus and Daedalus sensors including analog to digital converters to convert pixel data and various system voltages to digital form for readout and analysis, DAC channels for remote configuration of critical bias voltages, static random-access memories to buffer pixel data, RS422 and Gigabit Ethernet communications for remote access, and an FPGA to tie these components together. This document describes the FPGA electrical interfaces in detail to allow the reader a greater understanding of the device, and to facilitate implementation of custom software to control and manage it. The Version 4.0 Board is a continuation of the Nano-second Gated CMOS hardware design that retains much of the functionality of the Version 1.0 Board while adding features including a DAC instead of digital potentiometers, as well as sensors for pressure and radiation. The Version 4.0 board is intended for applications requiring tight form-factor enclosures. It is composed of two stacking boards; one holds the FPGA and regulators to power the various components of the board while the other contains the mating connector to the sensor, image-readoff ADCs, the DAC, and other components.

42 ENGINEERING↗
DOE OSTI2020

Covariance Testing at BNL [Slides]

This presentation discusses testing implemented in FUDGE & ADVANCE in regard to proper ENDF formatting: STAN, STANEF, CHECKR, FUDGE, NJOY, convert to BOXR format, plotting, endfcovreport.py which is an overview of ENDF library’s covariance contents, and covariance Matrix class implements that are the core of all GNDS covariance. Several checks are performed, including symmetric, real, positive definite, condition # of Eigenvalues, and too big uncertainties. The presentation states that it can fix too big and/or small uncertainties and that it can plot matrices. It states that further testing is recommended by CSEWG but has not yet been implemented and suggests methods such as documentation which needs a human to execute, completeness which is partially complete, and reasonableness which has not been started. The presentation concludes by elaborating on reasonableness and states that only cross section is documented, that similar limits on CP reactions cross sections are needed, that maximum limits are needed too which is tricky since ENDF assumes Normal PDFs, and that there may be other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Tailoring Cementitious Materials Towards Value-Added Use of Large CO 2 Volumes

Hydraulic cements with alternative chemistries were developed for large-volume and value-added use of carbon dioxide. The hydraulic cements were processed using the energy-efficient mechanochemical technique at room temperature and atmospheric pressure. Carbon dioxide was captured directly from combustion emissions during processing. For this purpose, mechanochemical processing of hydraulic cements was accomplished under a flow of combustion emissions prior to the release of emissions to the atmosphere. The process removed a significant fraction of carbon dioxide from combustion emissions. The hydraulic cements captured carbon dioxide at about 10% of their weight. Two hydraulic cement chemistries were developed, and their mechanochemical processing was successfully scaled-up. The characteristic feature of one chemistry was its relatively low (near-neutral) pH where the integration of carbon dioxide yielded clear value. The second cement chemistry was based on alkali activation of industrial wastes. This chemistry could make value-added use of carbon dioxide, but it had to be refined to control the pH drop caused by CO 2 integration. These two cements render binding effects upon hydration by forming a combination of stable carbonates and aluminosilicates or phosphates. Efforts to develop cement chemistries based solely on carbonates were not successful. The mechanochemical process was found to integrate carbon dioxide into the alternative cement chemistries in the form of disordered and metastable carbonates. During hydration reactions, the disordered/metastable carbonates are either transformed into stable carbonate phases with desired binding efforts, or carbonates get integrated into the primary inorganic binders (hydrates). When compared with mechanochemical processing in pure carbon dioxide, mechanochemical processing in combustion emissions produced hydraulic cements with improved engineering properties. The mechanochemical process was scaled-up, and its variables as well as the raw materials formulations were optimized for implementation at pilot scale. Scale-up was found to enhance the carbon capture potential and the engineering qualities of the resulting hydraulic cements. This was because scale-up raises the intensity of mechanical energy input to raw materials. Certain mechanochemical phenomena cannot be induced, irrespective of the cumulative mechanical energy input, unless the intensity of impact is raised above a minimum level that cannot be achieved in laboratory-scale implementation of the process. Due to this effect, the duration of the mechanochemical process as well as its energy demand could be reduced significantly (by an order of magnitude) upon transition from laboratory to pilot scale, with the hydraulic cements produced at pilot scale offering engineering properties that were superior to those realized in laboratory-scale mechanochemical processing. The hydraulic cements produced at pilot scale via mechanochemical processing under a flow of combustion emissions were thoroughly characterized. They were found to meet standard requirements for ‘General Use’ hydraulic cements. They were compatible with the industrial-scale concrete production and construction practices that are used with the currently prevalent hydraulic cements. The new hydraulic cements with integrated carbon dioxide were found to offer distinct advantages over the currently prevalent Portland cement in terms of net carbon footprint and energy content. The combined raw materials and energy costs of the new hydraulic cements are competitive, and major cost savings can be realized because of the simplified production process that significantly lowers the capital investment in cement production plants. Mechanochemical processing of the new hydraulic cements under a flow of combustion emissions can be implemented using some existing components of cement manufacturing plants; this facilitates adoption of the technology by the cement industry.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transient Efficiency, Flexibility, and Reliability Optimization of Coal-Fired Power Plants - Final Report

This program developed an advanced model-based monitoring and model-predictive control algorithms for a coal fired power plant (CFPP), and deployed these algorithms in a real-time platform to demonstrate performance benefits for transient flexibility and plant operation efficiency. More specifically, the objectives were successfully achieved through a combination of (i) developing a high-fidelity transient plant model in Apros, which was used as a high-fidelity plant simulation between $100-50\% TMCR$ where TMCR denotes the turbine maximum continuous rating, i.e., baseload, (ii) developing a very fast physics-based reduced-order model (ROM) of the plant, which ran more than $100\times$ faster than real-time, enabling its use as real-time embedded model for model-based estimation (MBE) and model predictive control (MPC) (iii) implementing a real-time MBE based on ROM using a robust unscented Kalman filter (UKF) to continuously tune the ROM to match the measurements from high-fidelity Apros plant model despite significant plant-model mismatch, and thus, obtain a Digital Twin of the plant (iv) designing and implementing a real-time MPC with dual objectives of transient plant load tracking with high ramp rates and minimizing coal consumption, i.e., improving plant efficiency in the baseload-partload operation range of $100-50\% TMCR$. Each key element above was developed and tested individually, and has been reported in corresponding Topical Reports. Finally, all the individual elements were integrated in an overall closed-loop system, that was successfully tested in desktop Simulink test harness simulations with ROM or high-fidelity model as the plant. Thereafter, the Simulink implementation was used to auto-generate C-code and deploy as real-time Docker microservice containers in Linux, and validate that they can run in real-time in the hardware-in-the loop (HIL) setup and produce the same results as in Simulink. The results of the integrated simulation tests in Simulink as well as the real-time HIL deployment are documented in this final report, showing good load tracking for load ramps at $3-4\%/min$ ramp rates, and achieving up to $5.5\%$ reduction in coal relative to baseline operation at $50\% TMCR$ load. The desktop and HIL simulations show successful performance of the overall model based estimation and control solution and achieve the key objectives of the program for flexible, efficient and reliable operation of subcritical coal fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

A Prototype Thick-Target Bremsstrahlung Model with Angularly-Dependent Emission in the MCNP6 ® Code

This document summarizes the current thick-target bremsstrahlung (TTB) model in MCNP and provides test results for an alternative implementation to improve the accuracy with reduced cost compared to full electron transport. It has been observed that the current TTB model produces inaccurate results in problems where the medium is thick with respect to electrons, but the photon distribution in the problem has a strong directionality. An example of such a simulation is detectors surrounding a metal target irradiated with a radiographic beam of high energy photons. The primary cause of this discrepancy is the current TTB method emits all bremsstrahlung photons in the same direction as the primary electron produced from each (γ, e ± ) interaction, leading to artificially forward peaked photon distributions for intermediate to high-energy incident photons. To improve the TTB model, we have implemented an angularly-dependent TTB model in a developer version of the MCNP6 ® code; for developers, this was done on the branch prototype/angular_ttb in the mcnp6 repo on bitbucket. The angularly-dependent TTB model accounts for the energy and scattering of electrons as they slow down in the current material, but does not sample the computationally expensive energy straggling, secondary electron events, and tracking electrons; this approach is significantly less computationally expensive than full electron transport and can be comparable to the original TTB method for problems with sufficiently complex materials and geometry. To evaluate the method, we have modeled a simple problem of a beam of 5 MeV photons incident on a sphere of plutonium surrounded by detectors at different deflection angles. For this problem, the angularly-dependent TTB produces a photon flux within 8.0% for a 90 degree deflection angle and 0.8% along the beam axis, as compared to the electron transport solution. This is an improvement compared to a 43% and 81% discrepancy with the original TTB method, respectively. The rest of this work includes the following: the first section details the current TTB treatment in MCNP, which has not been well documented elsewhere. Then, the modified TTB algorithm is detailed and the approximations compared to the condensed history algorithm are compared. Results are given comparing the two TTB methods to the condensed history transport algorithm. The appendix includes details for code developers on relevant electron transport implementation details and potential code improvements for future work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

The VTK-m Users' Guide (V.2.0)

High-performance computing relies on ever finer threading. Advances in processor technology include ever greater numbers of cores, hyperthreading, accelerators with integrated blocks of cores, and special vectorized instructions, all of which require more software parallelism to achieve peak performance. Traditional visualization solutions cannot support this extreme level of concurrency. Extreme scale systems require a new programming model and a fundamental change in how we design algorithms. To address these issues we created VTK-m: the visualization toolkit for multi-/many-core architectures. VTK-m supports a number of algorithms and the ability to design further algorithms through a top-down design with an emphasis on extreme parallelism. VTK-m also provides support for finding and building links across topologies, making it possible to perform operations that determine manifold surfaces, interpolate generated values, and find adjacencies. Although VTK-m provides a simplified high-level interface for programming, its template-based code removes the overhead of abstraction. VTK-m simplifies the development of parallel scientific visualization algorithms by providing a framework of supporting functionality that allows developers to focus on visualization operations. Consider the listings in Figure 1.1 that compares the size of the implementation for the Marching Cubes algorithm in VTK-m with the equivalent reference implementation in the CUDA software development kit. Because VTK-m internally manages the parallel distribution of work and data, the VTK-m implementation is shorter and easier to maintain. Additionally, VTK-m provides data abstractions not provided by other libraries that make code written in VTK-m more versatile.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Monitoring and Documentation of Forest Management Activities at Los Alamos National Laboratory

Forest management treatments at Los Alamos National Laboratory (LANL) are administered by the Emergency Management Division (EMD-DO) Wildland Fire Program within the Wildland Fire and Forest Management Program (WFFM Program or the Program). Open space forest inventory monitoring is conducted by Forest Health Program staff and students within the Environmental Stewardship Group (EPC-ES). The goals of the integrated Wildland Fire and Forest Management Program are to reduce impacts to Laboratory operations from climate driven events by implementing forest management and healthy forest initiatives such as forest thinning and erosion control. The LANL Wildland Fire Mitigation and Forest Health Plan (EMD PLAN-200, LA-UR-19-25122; LANL 2019) presents treatment standards for LANL property to meet the following goals: 1) Restore and maintain landscapes: LANL landscapes are resilient to disturbances; 2) Develop a fire-adapted community: LANL workforce, neighbors, and infrastructure can withstand a wildland fire without loss of life and property; 3) Ensure wildland fire mitigation implementation: All wildland fire mitigation working group organizations participate in making and implementing safe, effective, efficient risk-based wildland fire management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of FERC Order 2222 on DER Participation Rules in US Electricity Markets

Electricity markets in the bulk grid are beginning to implement market mechanisms that support the procurement of flexible capabilities from wide range of technologies, including distributed energy resources (DERs). The flexibility of these resources will help counterbalance supply uncertainties from large-scale integration of variable renewable generation. To encourage development of distributed and aggregated market participants, FERC Order 2222 was issued in September 2020 to require each Independent System Operator (ISO) in the US to implement rules that enable broader participation from aggregations of DERs in the bulk market. The following paper first describes the generic design of ISO markets before introducing the new market participation rules that ISOs have proposed for compliance with Order 2222. The paper then describes how software performance issues may continue to affect the eligibility requirements and offer structures for DER aggregations participating in ISOs, noting that continued research on computational methods may help reduce burdens for DER integration. The prospects for transmission and distribution system coordination is second major issue discussed, which will require minor changes to existing processes in the short term. In the longer term, there is more opportunity for more wide-ranging reforms, such as the development of a Distribution System Operator (DSO) framework. Newly proposed market rules may affect how Transactive Energy Systems (TES) will help facilitate efficient formation of DER aggregations and operation of the individual DERs within an aggregation. Within the TES context, the challenge is to fully understand how resource eligibility and operational and planning coordination methods will affect the design and implementation of TES.

24 POWER TRANSMISSION AND DISTRIBUTION↗