Search NASA⌕ Search

SEARCH · Search NASA

Results for “interactive data visualization”

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.

500 records · Page 28

Automated detection of double nuclei galaxies using GOTHIC and the discovery of a large sample of dual AGN

We present a novel algorithm to detect double nuclei galaxies (DNG) called Gothic (Graph-bOosTed iterated HIll Climbing) – that detects whether a given image of a galaxy has two or more closely separated nuclei. Our aim is to test for the presence of dual/multiple active galactic nuclei (AGN) in galaxies that visually represent a DNG. Although galaxy mergers are common, the detection of dual AGN is rare. Their detection is very important as they help us understand the formation of supermassive black hole (SMBH) binaries, SMBH growth and AGN feedback effects in multiple nuclei systems. There is thus a need for an algorithm to do a systematic survey of existing imaging data for the discovery of DNGs and dual AGNs. We have tested Gothic on an established sample of DNGs with a 100 per cent detection rate and subsequently conducted a blind search of 1 million SDSS DR16 galaxies (with spectroscopic data available) lying in the redshift range of z = 0 to 0.75. From the list of candidate DNGs found, we have detected 159 dual AGNs, of which 2 are triple AGN systems. Our results show that dual AGNs are not common, and triple AGN even rarer. The colour (u–r) magnitude plots of the DNGs indicate that star formation is quenched as the nuclei come closer and as the AGN fraction increases. The quenching is especially prominent for dual/triple AGN galaxies that lie at the extreme end of the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗

Scalable edge clustering of dynamic graphs via weighted line graphs

Timestamped relational datasets consisting of records (or connections) between pairs of entities are ubiquitous in network science. For applications like peer-to-peer communication, email, various social network interactions, and computer network security, it is useful to organize these records into groups based on how and when they are occurring. Weighted line graphs offer a natural way to model how records are related in such datasets but for large real-world graph topologies, building and utilizing the line graph is prohibitively expensive. Here, we present the framework to cluster the edges of a dynamic graph via the associated line graph that contains two major contributions. The first is a method to work with the line graph implicitly and the second is a distributed scale implementation of an agglomerative hierarchical graph clustering algorithm. We outline a novel hierarchical dynamic graph edge clustering approach that efficiently breaks massive relational datasets into small sets of edges containing events at various timescales. This is in stark contrast to traditional graph clustering algorithms that prioritize highly connected (clique-like) community structures. Our approach relies on constructing a sufficient subgraph of a weighted line graph and applying a hierarchical agglomerative clustering. This approach is related to scalable techniques from spatial clustering, nonlinear-dimension reduction, topological data analysis, and draws particular inspiration from HDBSCAN. As an edge clustering, this method yields an overlapping node clustering. Our algorithm is parallelizable and we demonstrate efficient clustering of a billion-scale, real-world dynamic graph into small edge sets that correlate in topology and time. The entire clustering process for a graph with tens of billions of edges takes just a few minutes of run time on 256 nodes of a distributed compute environment. We argue how the output of the edge clustering is useful for a multitude of data visualization and powerful machine learning tasks, both involving the original massive dynamic graph data and metadata associated with the nodes and edges. Finally, we describe how this approach can be extended to dynamic hypergraphs and dynamic graphs/hypergraphs with unstructured data living on vertices and edges.

Data Analysis↗

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗

Securing Grid-interactive Efficient Buildings (GEB) through Cyber Defense and Resilient System (CYDRES)

The DOE CYDRES project is driven by the urgent need to address critical research gaps in the domain of cyber-physical security of smart buildings, including Grid-interactive Efficient Buildings (GEBs). CYDRES, a real-time advanced building resilient platform, aims to enhance the cyber-attack-immune capabilities of buildings through multi-layered prevention, detection, and adaptation mechanisms. CYDRES consists of five key modules: a multi-layer network analyzer, an Automatic Fault Detection, Diagnosis, and Prognosis (AFDDP) framework, an intelligent mode selector, a cyber-resilient control framework, and a situation awareness platform. The Network Analyzer employs a data-driven framework that includes a protocol state learning tool and a CRF (Conditional Random Field) command validator. In Hardware-In-the-Loop (HIL) testbeds, it achieved 100% detection accuracy with a false alarm rate of 3%, validating its efficacy in identifying selected cyber-attacks. The AFDDP framework leverages pattern matching, PCA (Principal Component Analysis)-based strategies, and a DBN (Dynamic Bayesian Network)-based fault diagnosis approach to pinpoint the causes of physical system abnormalities using Building Automation System (BAS) data. In HIL experiments, the AFDDP module attained a detection accuracy of over 95% with a false alarm rate below 7%. Additionally, the fault detector utilized machine learning (Random Forest) and deep learning (Multi-Layer Perceptron) methods with acoustic sensor data to achieve a 100% fault detection accuracy in Heating, Ventilation, and Air-Conditioning (HVAC) equipment. The Mode Selector offered real-time impact analysis, allowing immediate actions to protect BASs in the face of emerging threats. The cyber-resilient control framework included an adaptive Model Predictive Control (MPC) and a measurement compensator, reducing temperature violations by up to 94% and improving the total demand flexibility by up to 70% in HIL experiments. Such HIL experiments covered a cyber-attack case and a physical fault case, showcasing CYDRES’ efficiency in maintaining operational continuity during threats. The situation awareness platform in Grafana enhanced real-time threat detection and response visualization, augmenting the operational awareness for building operators. CYDRES demonstrated high technical effectiveness in various test scenarios, particularly in HIL environments. The project's phased development approach ensured efficient use of resources, highlighting its practical feasibility and readiness for commercialization. By enhancing the security and resilience of building operations, CYDRES represents a significant advance in mitigating risks associated with cyber-physical systems, thereby enhancing public confidence in the safety of modern building infrastructure. Future directions for the project include expanding testing protocols, refining AFDDP methodologies, exploring more comprehensive resilient control strategies, and testing in real commercial buildings.

42 ENGINEERING↗

The Trip Itinerary Optimization Platform: A Framework for Personalized Travel Information

The New Concepts Incubator team at the National Renewable Energy Laboratory (NREL) developed a three-stage online platform for travel diary collection, personal travel plan optimization and travel itinerary visualization. In the first stage, users provide a travel diary for the previous day through an interactive map and calendar interface and survey for travel attitudes and behaviors. One or more days later, users are invited via email to engage in a second stage where they view a personal mobility dashboard displaying recommended travel itineraries generated from a novel framework that optimizes travel outcomes over a sequence of interrelated trips. A week or more after viewing these recommended travel itineraries on the dashboard, users are emailed again to engage in a third stage where they complete a final survey about travel attitudes and behaviors. A usability study of the platform conducted online showed that, in general, users found the system valuable for informing their travel decisions. A total of 274 individuals were recruited through Amazon Mechanical Turk, an online survey platform, to participate in a transportation study using this platform. On average, the platform distilled 65 feasible travel plans per individual into two recommended itineraries, each optimal according to one or more outcomes and dependent on the fixed times and locations from the travel diary. For 45 percent of users, the trip recommendation algorithm returned only a single, typically automobile-centric, itinerary because there were no other viable alternative transportation modes available. Platform users generally agreed that the dashboard was enjoyable and easy to use, and that it would be a helpful tool in adopting new travel behaviors. Users generally agreed most that the time, cost and user preferred recommendations 'made sense' to them, and were most willing to implement these itineraries. Platform users typically expressed low willingness to try the carbon and calories optimized itineraries. Of the platform users who viewed the dashboard, 13 percent reported changing their travel behavior, most adopting the time, calories or carbon optimized itineraries. While the algorithm incorporates a wealth of travel data obtained from online APIs pertaining to a travelers route such as historic traffic condition data, public transit time-tables, and bike path routes, open-ended responses from users expressed an interest in the integration of even more fine-grained traffic data and the ability to dynamically model the effect of changes in travel times. Users also commonly expressed concerns over the safety of walking and biking recommendations. Responses indicate that more information about the amenities available to cyclists and pedestrians (sidewalks, shade from trees, access to food) and routes that avoid areas of perceived elevated danger would reduce barriers to implementing these recommendations. More accurate representations of personal vehicle trips (based on vehicle make and model, implications of parking) and the identification of routes that optimize caloric intensity (seeking out elevation changes or longer walks to public transit) are promising avenues for future research.

24 POWER TRANSMISSION AND DISTRIBUTION↗

l3d2vtk: An MCNPTools Utility to Enable LNK3DNT File Visualization & Post-Processing

This paper describes recent work to develop and incorporate a utility program into the MCNPTools package that enables LNK3DNT file visualization & post-processing. A LNK3DNT file is a binary Cartesian, cylindrical, or spherical structured mesh file that describes material and material relative density distributions by mesh voxel. This mesh file format is principally used by the PARTISN discrete-ordinates radiation-transport code. However, the MCNP ® Monte Carlo radiation-transport code can use LNK3DNT files to specify geometry, and the MCNP code can produce LNK3DNT files representing a user-defined region of the MCNP constructive solid geometry (CSG). Because of the use of LNK3DNT files as MCNP input and output and with other computer codes, improving the ability to visualize and/or quantitatively post-process LNK3DNT files benefits end users. In particular, it is helpful to visually verify that the material and density distributions are as expected. Also, an analyst can orient viewers with respect to other in formation, e.g., by co-plotting a voxelized representation of MCNP CSG with mesh tally and/or particle track results. However, until now, the only capability to interactively interrogate LNK3DNT files is through a 2-D slice plotter provided by the Monte Carlo Application ToolKit (MCATK), which is not currently distributed outside of Los Alamos National Laboratory. The utility created as part of this work converts LNK3DNT binary files to an ASCII XML-based StructuredMesh VTK (.vts) file. The geometries supported are 1-D spherical and cylindrical ( r ), 1-D Cartesian ( x ), 2-D cylindrical ( r-z and r-θ ) and Cartesian ( x-y ), and 3-D cylindrical ( r-z-θ ) and Cartesian ( x-y-z ). This file format is chosen to minimize external dependencies (e.g., on an HDF5 library) and to provide a file that can be interrogated by the user in a text editor or standard visualization applications (e.g., ParaView or VisIt). However, other file formats may provide more utility such as XDMF, which can use HDF5 and is discussed further as something to explore as future work. This paper is organized as follows: a brief introduc tion to the MCNPTools library and utility programs is given, the LNK3DNT-to-VTK (l3d2vtk) utility functionality is de scribed, and then an example post-processing workflow using the ParaView application is given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

Galaxy Zoo: Major Galaxy Mergers Are Not a Significant Quenching Pathway

We use stellar mass functions to study the properties and the significance of quenching through major galaxy mergers. In addition to SDSS DR7 and Galaxy Zoo 1 data, we use samples of visually selected major galaxy mergers and post-merger galaxies. We determine the stellar mass functions of the stages that we would expect major-merger-quenched galaxies to pass through on their way from the blue cloud to the red sequence: (1) major merger, (2) post-merger, (3) blue early type, (4) green early type, and (5) red early type. Based on their similar mass function shapes, we conclude that major mergers are likely to form an evolutionary sequence from star formation to quiescence via quenching. Relative to all blue galaxies, the major-merger fraction increases as a function of stellar mass. Major-merger quenching is inconsistent with the mass and environment quenching model. At z∼0, major-merger-quenched galaxies are unlikely to constitute the majority of galaxies that transition through the green valley. Furthermore, between z∼0−0.5, major-merger-quenched galaxies account for 1%–5% of all quenched galaxies at a given stellar mass. Major galaxy mergers are therefore not a significant quenching pathway, neither at z∼0 nor within the last 5 Gyr. The majority of red galaxies must have been quenched through an alternative quenching mechanism that causes a slow blue to red evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Augmented Reality Tools for Nuclear Facilities - 18318

In this work we explore the possibility of improving nuclear facility operations through the use of emerging augmented reality technology. Specifically we focus on the use of augmented reality headsets. Augmented reality opens up the possibility for easy, on-the-fly access to documentation related to canisters, rooms, and carts. It also can be used to help plan complicated material moves in ways that do not violate nuclear criticality safety regulations. Modern augmented reality technology can also be interacted with using gaze, gesture and voice commands, thus making it particularly attractive for facilitating the documentation of nuclear material operations and work verification. An augmented reality-based nuclear facility criticality safety system would have a number of benefits for operations in nuclear facilities. This work could be used across the DOE complex to ensure technicians in gloveboxes work safely and securely. It could greatly reduce downtime associated with administrative criticality incidents. It would be a step forward in the process of improving the quality of documentation of the work that occurs in nuclear facilities while simultaneously increasing the ease with which this documentation can be collected. The augmented reality effectively endows the glovebox technicians with a sense for administrative nuclear criticality safety and the reactivity changes caused directly by their actions. Augmented reality also has great potential for facilitating inspections of nuclear infrastructure. Examples of Augmented Reality applications that facilitate infrastructure inspection and work verification will also be presented. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Galaxy Zoo DECaLS: Detailed visual morphology measurements from volunteers and deep learning for 314,000 galaxies

We present Galaxy Zoo DECaLS: detailed visual morphological classifications for Dark Energy Camera Legacy Survey images of galaxies within the SDSS DR8 footprint. Deeper DECaLS images (r = 23.6 versus r = 22.2 from SDSS) reveal spiral arms, weak bars, and tidal features not previously visible in SDSS imaging. To best exploit the greater depth of DECaLS images, volunteers select from a new set of answers designed to improve our sensitivity to mergers and bars. Galaxy Zoo volunteers provide 7.5 million individual classifications over 314 000 galaxies. 140 000 galaxies receive at least 30 classifications, sufficient to accurately measure detailed morphology like bars, and the remainder receive approximately 5. All classifications are used to train an ensemble of Bayesian convolutional neural networks (a state-of-the-art deep learning method) to predict posteriors for the detailed morphology of all 314 000 galaxies. We use active learning to focus our volunteer effort on the galaxies which, if labelled, would be most informative for training our ensemble. When measured against confident volunteer classifications, the trained networks are approximately 99 per cent accurate on every question. Morphology is a fundamental feature of every galaxy; our human and machine classifications are an accurate and detailed resource for understanding how galaxies evolve.

79 ASTRONOMY AND ASTROPHYSICS↗

Federal Emergency Management Agency Mobile Housing Unit Running Gear Testing (Final Report)

The Federal Emergency Management Agency (FEMA) stages manufactured housing units (MHUs) at key locations to be used during disasters when other housing accommodations cannot be arranged for individuals or families. Due to the unpredictability of natural disasters and the limitations of the manufactured housing industry to produce large numbers of MHUs in very short periods of time, it is necessary for FEMA to store a large number of units for many months or years. With these long storage periods come degradation issues with frames, tires, bearing, and braking systems. These types of issues are also of interest to the U.S. Department of Transportation (DOT) Federal Motor Carrier Safety Administration (FMCSA) because regulated carriers normally transport these units and 49 CFR 393.75 Tires specifies tire loading restrictions for manufactured homes. To compare the running gear (tires, brakes, drums, bearings, axles, springs, and spring mounting hardware) currently specified in the 2014 Rugged Based Performance Requirements (RBPR) and the newly recommended running gear specified in the upcoming revision to the RBPR, the Oak Ridge National Laboratory (ORNL) outlined a Longevity Test Plan, which focused on the assessment of component longevity (failure rates by component type). The Test Plan also studied the forces imparted to the MHU from the running gear from interaction with the roadway for both the running gear in the 2014 RBPR and the running gear specified in the upcoming RBPR revision. Since the assessment of component longevity and failures are also issues of concern for FMCSA, both agencies collaborated to support the research presented in this report. The tests and data collection were conducted at the Navistar Proving Grounds (NPG) in New Carlisle (near South Bend), Indiana, during the May-November 2018 period. Four alternatives for the suspension system of a typical three-bedroom MHU were considered in this test. Alternative 0 was the MHU that was selected from the FEMA staging area in Selma, Alabama. Empirical data had shown that, on average, this alternative experiences one tire failure every 150 miles (or 0.0067 tire failures per mile traveled). Before testing, Alternative 0 was equipped with similar running gear hardware, although in new condition and installed following the manufacturer’s specifications. This became the Baseline Unit, or Alternative 1. Another unit, also selected from the FEMA staging area in Selma, Alabama, was equipped with better axles (i.e., axles that were rated for 7,000lbs as opposed to 6,000lbs for the Baseline Unit) and better tires. This became the Upgrade Unit, or Alternative 2. At the end of the Longevity Test, the Upgrade Unit was mounted with radial tires. This became Alternative 3.Alternative 1 (Baseline Unit) had, on average, one tire failure every 251 miles (or 0.0040 tire failures per mile traveled). So, by simply following the correct manufacturer’s specifications for the assembly of the transportation system, a reduction of 41% in the per mile tire failure was achieved. In the same test, Alternative 2 (Upgrade Unit) showed, on average, one tire failure every 728 miles traveled (or 0.0014 tire failure s per mile traveled). When the Upgrade Unit was mounted with radial tires in the last part of the test, no tire failures were observed for the 2,400 miles tested. That is, in the case of the radial tires (Alternative 3) the tire failure was less than one per 2,400 miles traveled (or less than 0.0004 tire failures per mile traveled). The data collected for the Longevity Test showed that the maximum accelerations were registered in all cases at the center of the MHU in the vertical direction. The overall maximum was registered by the Baseline Unit, Alternative 1. If that value is normalized to 100%, then the Upgrade Unit (Alternative 2) showed, on average, 81% at the same location, and Alternative 3 (radial-tire test) 59%. That is, the upgraded transportation system showed a reduction of 19% and 41% in the vertical accelerations transmitted to MHU at its center when regular and radial tires, respectively, were used. As part of the Test Protocol, NPG conducted a visual inspection of the MHUs and transportation system at the end of each test day. ORNL reviewed all the inspection reports in an attempt to determine if there was any difference in the damage that the Baseline and Upgrade MHUs showed during the test. The Upgrade Unit always showed lower accelerations transmitted to the MHU than the Baseline Unit. In most cases, these differences were statistically significant. The FEMA inspection report, conducted by a FEMA certified inspector at the end of the tests, included a shorter list of damaged items for the Baseline Unit than for the Upgrade Unit, although similar type of damages was described for both units. This may be an indication that, although lower than those of the Baseline Unit, the accelerations transmitted by the Upgrade Unit were above the threshold at which MHU elements start to fail.

42 ENGINEERING↗

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗

Size-resolved ice nucleating particle (INP) concentrations from the MOSAiC campaign

Sample collection: Ground-based immersion mode ice nucleation measurements were conducted on samples collected at the U.S. DOE ARM AMF2 onboard the R/V Polarstern P-deck during MOSAiC from 23 Oct 2019 to 01 Oct 2020. Samples were collected using a Davis Rotating-drum Unit for Monitoring (DRUM; DRUMAir, LLC) cascading impactor in 4 stages (i.e., size bins). Instrument information on the model used, the DA400, can be found at the DRUMAir website (http://www.drumair.com/). Stages A, B, C, and D had particle size cuts at 2.96, 1.21, 0.34, and 0.15 &micro;m in particle diameter, respectively. Strips of sample substrate (PFA, perfluoroalkoxy; McMaster) were adhered to each disc in each stage and coated with petrolatum to enable particles to stick on impact. The discs move slowly over time, such that aerosol loading is streaked onto the PFA. Every 24 hours, a blank is created on the PFA to separate daily samples. During MOSAiC, the pump typically pulled 20 &ndash; 35 lpm (average of 32 lpm) at the inlet with all the discs in place in each stage. Daily samples were collected for each of the four stages, equaling 28,800 &ndash; 50,400 total L of air (average of 46,080 L total L of air) per sample. Discs rotated for approximately 24 days before the sampling per strip was complete and changed to a new disc. Discs were stored frozen after collection. Sample preparation for analysis: Following collection and transport from MOSAiC to the Kreidenweis/DeMott laboratory, care of PI Creamean, in Fort Collins, CO, strips were removed from the discs and cut up to separate daily samples. Each daily portion of the strip was stored in separate sterile 15-mL centrifuge tubes and stored in a standard commercial freezer until analysis (approximately 4 &ndash; 12 months after collection). Sample analysis: Drop freezing assays (DFAs) were conducted on one daily sample set every 3 &ndash; 4 days of the study for all four stages using the Colorado State University (CSU) cold plate (CP). Sample set days were chosen to cover the entire year of MOSAiC and avoiding days where the AOS purge system was activated. Immediately after removal from the freezer, 2 mL of 0.1-&micro;m filtered deionized water was added with a single-use sterile syringe directly to the tubes, then mixed using a vortex mixer for at least 30 minutes at 500 rpm to resuspend particles from the PFA into the deionized water. Following sample preparation, another single-use syringe was used to create 100 aliquots on a 3-inch diameter copper plate, then covered with a transparent plastic dome. The drops were approximately the same volume (2.5 &micro;L each), though careful inspection. The plate was cooled at approximately 3 &ndash; 10 K min-1 from room temperature until all drops froze on the plate. Temperature was measured by an Omega thermocouple probes and meter. Drop freezing was visually detected but recorded through monitoring software to provide the time frozen, channel of the meter used, temperature, and cooling rate for each drop. Not all 100 drops were always detected; the total number of rows in each data file equals the number of drops recorded, which typically was > 80 %. Each sample was tested three times (i.e., a new set of 100 drops was created for each test). Data provided are files containing these parameters. From the freezing temperatures recorded, the fraction frozen was and INP concentrations were calculated per L of air using the equation from Vali (1971): [INPs] = &minus;(ln(f)/Valiquot)*(Vsuspension/Vair) where f is the proportion of droplets not frozen, Valiquot is the volume of each aliquot (i.e., 2.5e-6 L), Vsuspension is the volume of the suspension (i.e., 0.002 L), and Vair is the volume of air per sample (i.e., 36000 L). Data usage guidance: There are three files included with this dataset: (1) dumcp_inp_concentrations.csv which includes cumulative INP concentrations binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C for all samples and size ranges (i.e., DRUM stages) processed via DFA using the CSU CP; (2) dumcp_inp_standard_deviation.csv which includes the standard deviations for all samples and temperatures in dumcp_inp_concentrations.csv; and (3) dumcp_blank_concentrations_standard_deviation.csv which includes the INP concentrations and standard deviations for the deionized water and PFA blanks binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C. The first column in each file represents the freezing temperature in &deg;C binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C (first cell in this column = header = &ldquo;temp&rdquo;). The second column in each file represents the standard deviation of each temperature bin in &deg;C (first cell in this column = header = &ldquo;temperr&rdquo;). For file (1) dumcp_inp_concentrations.csv, the remaining columns indicate the INP concentrations per L of air for each temperature bin. For each sample column, the header in the first row indicates the stage (i.e., size bin as &ldquo;A_&rdquo;, &ldquo;B_&rdquo;, &ldquo;C_&rdquo;, or &ldquo;D_&rdquo;) and date (&ldquo;_yyyymmdd) of collection. For file (2) dumcp_inp_standard_deviation.csv, the remaining columns indicate the INP concentration standard deviation per L of air for each temperature bin. For each sample column, the header in the first row indicates the stage (i.e., size bin as &ldquo;A_&rdquo;, &ldquo;B_&rdquo;, &ldquo;C_&rdquo;, or &ldquo;D_&rdquo;) and date (&ldquo;_yyyymmdd) of collection. For file (3) dumcp_blank_concentrations_standard_deviation.csv, the remaining columns indicate the INP concentrations and standard deviations for the deionized water blank (&ldquo;di_0.1&rdquo; and &ldquo;di_0.1sd&rdquo; for the concentrations and standard deviations, respectively) and PFA blank (&ldquo;pfa&rdquo; and &ldquo;pfasd&rdquo; for the concentrations and standard deviations, respectively). Please see Creamean et al. (2019a; 2018a; 2018b; 2019b) for additional details on DRUM sample collection and offline processing. For questions or help with data interpretation, contact Dr. Jessie Creamean (jessie.creamean@colostate.edu). References: Creamean, J. M., Mignani, C., Bukowiecki, N., and Conen, F.: Using spectra characteristics to identify ice nucleating particle populations during winter storms in the Alps, Atmos. Chem. Phys., 19, 8123&ndash;8140, 2019a. Creamean, J. M., Kirpes, R. M., Pratt, K. A., Spada, N. S., Maahn, M., de Boer, G., Schnell, R. C., and China, S.: Marine and terrestrial influences on ice nucleating particles during continuous springtime measurements in an Arctic oilfield location, Atmospheric Chemistry and Physics, 18, 18023-18042, 10.5194/acp-18-18023-2018, 2018a. Creamean, J. M., Primm, K. M., Tolbert, M. A., Hall, E. G., Wendell, J., Jordan, A., Sheridan, P. J., Smith, J., and Schnell, R. C.: HOVERCAT: A novel aerial system for evaluation of aerosol-cloud interactions, Atmospheric Measurement Techniques, 11, 3969-3985, 10.5194/amt-11-3969-2018, 2018b. Creamean, J. M., Cross, J. N., Pickart, R., McRaven, L., Lin, P., Pacini, A., Hanlon, R., Schmale, D. G., Ceniceros, J., Aydell, T., Colombi, N., Bolger, E., and DeMott, P. J.: Ice nucleating particles carried from below a phytoplankton bloom to the Arctic atmosphere, Geophys. Res. Lett., 46, 8572-8581, https://doi.org/10.1029/2019GL083039, 2019b. Vali, G.: Quantitative Evaluation of Experimental Results an the Heterogeneous Freezing Nucleation of Supercooled Liquids, Journal of the Atmospheric Sciences, 28, 402-409, 10.1175/1520-0469(1971)028<0402:qeoera>2.0.co;2, 1971.

4-stage DRUM aerosol collector (DA400),Ice nuclei,↗

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↗