The credibility of an engineering model is of critical importance in large-scale projects. How concerned should an engineer be when reusing someone else's model when they may not know the author or be familiar with the tools that were used to create it? In this report, the authors advance engineers' capabilities for assessing models through examination of the underlying semantic structure of a model--the ontology. This ontology defines the objects in a model, types of objects, and relationships between them. In this study, two advances in ontology simplification and visualization are discussed and are demonstrated on two systems engineering models. These advances are critical steps toward enabling engineering models to interoperate, as well as assessing models for credibility. For example, results of this research show an 80% reduction in file size and representation size, dramatically improving the throughput of graph algorithms applied to the analysis of these models. Finally, four future problems are outlined in ontology research toward establishing credible models--ontology discovery, ontology matching, ontology alignment, and model assessment.
We present a new pre-processor tool written in Python that creates multicontinuum meshes for PFLOTRAN to simulate two-phase flow and transport in both the fracture and matrix continua. We discuss the multicontinuum modeling approach to simulate potentially mobile water and gas in the fractured volcanic tuffs at Aqueduct Mesa, at the Nevada National Security Site.
This report consolidates the requirements for an Exceedance Response Action (ERA) Level 1 and ERA Level 2 Action Plan for pH. A discharger’s baseline status for any given parameter changes to Level 1 status if sampling results indicate a Numeric Action Level (NAL) exceedance for that same parameter. NAL exceedance can be either of the following: (1) Instantaneous maximum NAL exceedance: Occurs when two or more analytical results for any single parameter within a reporting year exceed the instantaneous maximum NAL (for example for pH a value less than 6 or a value greater than 9); and (2) Annual NAL exceedance: Occurs when the average of all the analytical results for a parameter within a reporting year exceeds the annual NAL. A Discharger’s Level 1 status for any given parameter changes to Level 2 status if sampling results indicate a Numeric Action Level (NAL) exceedance for the same parameter while a Discharger is in Level 1.
We present a new method to discriminate between earthquakes and buried explosions using observed seismic data. The method is different from previous seismic discrimination algorithms in two main ways. First, we use seismic spatial gradients, as well as the wave attributes estimated from them (referred to as gradiometric attributes), rather than the conventional three-component seismograms recorded on a distributed array. The primary advantage of this is that a gradiometer is only a fraction of a wavelength in aperture com¬pared with a conventional seismic array or network. Second, we use the gradiometric attributes as input data into a machine learning algorithm. The resulting discrimination algorithm uses the norms of truncated principal components obtained from the gradio- metric data to distinguish the two classes of seismic events. Using high-fidelity synthetic data, we show that the data and gradiometric attributes recorded by a single seismic gra¬diometer performs as well as a conventional distributed array at the event type discrimi¬nation task.
The Geophysical Monitoring System (GMS) State-of-Health User Interface (SOH UI) is a web-based application that allows a user to view and acknowledge the SOH status of stations in the GMS system. The SOH UI will primarily be used by the System Controller, who monitors and controls the system and external data connections. The System Controller uses the station SOH UIs to monitor, detect, and troubleshoot problems with station data availability and quality.
This report is intended to detail the findings of our investigation of the applicability of machine learning to the task of aftershock identification. The ability to automatically identify nuisance aftershock events to reduce analyst workload when searching for events of interest is an important step in improving nuclear monitoring capabilities and while waveform cross - correlation methods have proven successful, they have limitations (e.g., difficulties with spike artifacts, multiple aftershocks in the same window) that machine learning may be able to overcome. Here we apply a Paired Neural Network (PNN) to a dataset consisting of real, high quality signals added to real seismic noises in order to work with controlled, labeled data and establish a baseline of the PNN's capability to identify aftershocks. We compare to waveform cross - correlation and find that the PNN performs well, outperforming waveform cross - correlation when classifying similar waveform pairs, i.e., aftershocks.
The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.
Computing k-cores on graphs is an important graph mining target as it provides an efficient means of identifying a graph's dense and cohesive regions. Computing k-cores on hypergraphs has seen recent interest, as many datasets naturally produce hypergraphs. Maintaining k-cores as the underlying data changes is important as graphs are large, growing, and continuously modified. In many practical applications, the graph updates are bursty, both with periods of significant activity and periods of relative calm. Existing maintenance algorithms fail to handle large bursts, and prior parallel approaches on both graphs and hypergraphs fail to scale as available cores increase.We address these problems by presenting two parallel and scalable fully-dynamic batch algorithms for maintaining k-cores on both graphs and hypergraphs. Both algorithms take advantage of the connection between k-cores and h-indices. One algorithm is well suited for large batches and the other for small. We provide the first algorithms that experimentally demonstrate scalability as the number of threads increase while sustaining high change rates in graphs and hypergraphs.
While a great deal of research has been performed to quantify and characterize the wave energy resource, there are still open questions about how a wave energy developer should use this wave resource information to design a wave energy converter device to suit a specific environment or, alternatively, to assess potential deployment locations. It is natural to focus first on the impressive magnitudes of power available from ocean waves, and to be drawn to locations where mean power levels are highest. However, a number of additional factors such as intermittency and capacity factor may be influential in determining economic viability of a wave energy converter, and should therefore be considered at the resource level, so that these factors can influence device design decisions. This study examines a set of wave resource metrics aimed towards this end of bettering accounting for variability in wave energy converter design. The results show distinct regional trends that may factor into project siting and wave energy converter design. Although a definitive solution for the optimal size of a wave energy converter is beyond the reaches of this study, the evidence presented does support the idea that smaller devices with lower power ratings may merit closer consideration.
This article examines the qualitative features of an anelasticity model associated with the bowing of dislocations in the presence of phase transition. A simple physically plausible mechanism is introduced to describe the interaction of anelasticity and the transformation. Varying the anelastic parameters results in strong differences in the deviatoric stress response. The model is applied to study the behavior of tin (Sn) and compared to data from ramp driven compression-release experiments. Tin exhibits a complex phase diagram within a relatively accessible range of temperature and pressures and the characterization of its phases is considered an open problem with significant scientific merit. The coupling between anelasticity, plasticity, and phase transformation contributes to release wave features traditionally associated with the phase transition effect alone suggesting the importance of accounting for the effects jointly. Posterior distributions of the plastic and anelastic parameters are computed using Bayesian-inference-based methods, further highlighting the importance of anelasticity in this regime.
The U.S. Strategic Petroleum Reserve is a crude oil storage system run by the U.S. Department of Energy. The reserve consists of 60 active storage caverns spread across four sites in Louisiana and Texas, near the Gulf of Mexico. Beginning in 2016, the SPR began executing U.S. congressionally mandated oil sales. The configuration of the reserve, with a total capacity of greater than 700 MMB, requires raw water to be used instead of saturated brine for oil withdrawals such as for sales. All sales will produce leaching within the caverns used for oil delivery. Twenty-five caverns had a combined total of over 39 MMB of water injected in CY 20 as part of the Exchange for Storage program; oil was withdrawn in the same manner as for congressionally mandated sales. Leaching effects were monitored in these caverns to understand how the oil withdrawals may impact the long-term integrity of the caverns. While frequent sonars are the best way to monitor changes in cavern shape, they can be resource intensive for the number of caverns involved in sales and exchanges. An intermediate option is to model the leaching effects and see if any concerning features develop. The leaching effects were modeled here using the Sandia Solution Mining Code (SANSMIC) . The results indicate that leaching induced features are not of concern in the majority of the caverns, 19 of 25. Six caverns, BH-107, BH-113, BH-114, BM-4, BM-106, and WH-114 have features that may grow with additional leaching and should be monitored as leaching continues in those caverns. Ten caverns had post sale sonars that were compared with SANSMIC results. SANSMIC was able to capture the leaching well , particularly the formation of shelves and flares. A deviation in the SANSMIC and sonar cavern shapes was observed near the cavern floor in caverns with significant floor rise, a process not captured by SANSMIC. These results suggest SANSMIC is a useful tool for monitoring changes in cavern shape due to leaching effects related to sales and exchanges.
TEMPI provides a transparent non-contiguous data-handling layer compatible with various MPIs. MPI Datatypes are a powerful abstraction for allowing an MPI implementation to operate on non-contiguous data. CUDA-aware MPI implementations must also manage transfer of such data between the host system and GPU. The non-unique and recursive nature of MPI datatypes mean that providing fast GPU handling is a challenge. The same noncontiguous pattern may be described in a variety of ways, all of which should be treated equivalently by an implementation. This work introduces a novel technique to do this for strided datatypes. Methods for transferring non-contiguous data between the CPU and GPU depends on the properties of the data layout. This work shows that a simple performance model can accurately select the fastest method. Unfortunately, the combination of MPI software and system hardware available may not provide sufficient performance. The contributions of this work are deployed on OLCF Summit through an interposer library which does not require privileged access to the system to use
In April of 2020, Sandia National Laboratories had an urgent need to identify and manage the data that could be used to create mobile applications, models, reports, and visualizations to assist management in safely bringing the workforce onsite during the COVID-19 pandemic. Multiple divisions volunteered to design and build software solutions; meanwhile, requests for new data sources, including duplicate requests, were inundating Information Technology (IT) and data owners. The Enterprise Data Governance Team was assigned to resolve obtaining and accessing new sources of data in an accelerated timeframe. Through successful collaboration with multiple stakeholders and domain owners across Sandia, the Enterprise Data Governance Team rapidly developed a centralized data strategy and solution for use in safeguarding the Sandia workforce during the COVID-19 pandemic. This foundation enabled teams to successfully develop solutions, including reports for executives and management as well as the data for modeling and scientific analysis.
Computing k-cores on graphs is an important graph mining target as it provides an efficient means of identifying a graph's dense and cohesive regions. Computing k-cores on hypergraphs has seen recent interest, as many datasets naturally produce hypergraphs. Maintaining k-cores as the underlying data changes is important as graphs are large, growing, and continuously modified. In many practical applications, the graph updates are bursty, both with periods of significant activity and periods of relative calm. Existing maintenance algorithms fail to handle large bursts, and prior parallel approaches on both graphs and hypergraphs fail to scale as available cores increase.We address these problems by presenting two parallel and scalable fully-dynamic batch algorithms for maintaining k-cores on both graphs and hypergraphs. Both algorithms take advantage of the connection between k-cores and h-indices. One algorithm is well suited for large batches and the other for small. We provide the first algorithms that experimentally demonstrate scalability as the number of threads increase while sustaining high change rates in graphs and hypergraphs.
Digital image correlation (DIC) is an established test technique in several fields including quasi-static displacement measurements. Recently there has been growing interest in using DIC to measure structural dynamic response and even extract modal parameters from that information. While high-speed cameras have become more ubiquitous, there are no commercial end-to-end packages for modal analysis based on image data, particularly when combined with traditional data acquisition systems. As such, the practitioner is left to develop several key data processing capabilities, hardware interface equipment, and testing practices themselves. This work highlights several practical aspects that have been encountered while establishing DIC as a viable modal testing capability in a laboratory environment.
The reversible computation paradigm aims to provide a new foundation for general classical digital computing that is capable of circumventing the thermodynamic limits to the energy efficiency of the conventional, non-reversible digital paradigm. However, to date, the essential rationale for, and analysis of, classical reversible computing (RC) has not yet been expressed in terms that leverage the modern formal methods of non-equilibrium quantum thermodynamics (NEQT). In this paper, we begin developing an NEQT-based foundation for the physics of reversible computing. We use the framework of Gorini-Kossakowski-Sudarshan-Lindblad dynamics (a.k.a. Lindbladians) with multiple asymptotic states, incorporating recent results from resource theory, full counting statistics and stochastic thermodynamics. Important conclusions include that, as expected: (1) Landauer’s Principle indeed sets a strict lower bound on entropy generation in traditional non-reversible architectures for deterministic computing machines when we account for the loss of correlations; and (2) implementations of the alternative reversible computation paradigm can potentially avoid such losses, and thereby circumvent the Landauer limit, potentially allowing the efficiency of future digital computing technologies to continue improving indefinitely. We also outline a research plan for identifying the fundamental minimum energy dissipation of reversible computing machines as a function of speed.
Radar is by its basic nature a ranging instrument. If radar range and range-rate measurements from multiple directions can be made and assembled, then multilateration allows locating a feature common to the set of Synthetic Aperture Radar (SAR) images to an accurate 3-D coordinate. The ability to employ effective multilateration algorithms is highly dependent on the geometry of the data collections, and the accuracy with which relative range measurements can be made. The problem can be cast as a least-squares exercise, and the concept of Dilution of Precision can describe the accuracy and precision with which a 3-D location can be made.
The electric grid is becoming increasingly cyber-physical with the addition of smart technologies, new communication interfaces, and automated grid-support functions. Because of this, it is no longer sufficient to only study the physical system dynamics, but the cyber system must also be monitored as well to examine cyber-physical interactions and effects on the overall system. To address this gap for both operational and security needs, cyber-physical situational awareness is needed to monitor the system to detect any faults or malicious activity. Techniques and models to understand the physical system (the power system operation) exist, but methods to study the cyber system are needed, which can assist in understanding how the network traffic and changes to network conditions affect applications such as data analysis, intrusion detection systems (IDS), and anomaly detection. In this paper, we examine and develop models of data flows in communication networks of cyber-physical systems (CPSs) and explore how network calculus can be utilized to develop those models for CPSs, with a focus on anomaly and intrusion detection. This provides a foundation for methods to examine how changes to behavior in the CPS can be modeled and for investigating cyber effects in CPSs in anomaly detection applications.
Performing terrain classification with data from heterogeneous imaging modalities is a very challenging problem. The challenge is further compounded by very high spatial resolution. (In this paper we consider very high spatial resolution to be much less than a meter.) At very high resolution many additional complications arise, such as geometric differences in imaging modalities and heightened pixel-by-pixel variability due to inhomogeneity within terrain classes. In this paper we consider the fusion of very high resolution hyperspectral imaging (HSI) and polarimetric synthetic aperture radar (PolSAR) data. We introduce a framework that utilizes the probabilistic feature fusion (PFF) one-class classifier for data fusion and demonstrate the effect of making pixelwise, superpixel, and pixelwise voting (within a superpixel) terrain classification decisions. We show that fusing imaging modality data sets, combined with pixelwise voting within the spatial extent of superpixels, gives a robust terrain classification framework that gives a good balance between quantitative and qualitative results.
The supercritical carbon dioxide (s-CO2) Brayton cycle is currently being explored as a replacement for the steam Rankine cycle due to its potential for higher efficiency and lower cycle cost. 316 stainless steel is a candidate alloy for use in s-CO2 up to roughly 600 °C, but the mechanical effects of prolonged exposure of base and welded material in s-CO2 have not been analyzed. The potential for carburization makes this an important concern for the implementation of 316 and similar austenitic stainless steels in the s-CO2 environment. In this study, welded and base material of two types of 316–316L and 316H–were exposed in either s-CO2 or argon at 550 °C or 750 °C for 1000 h. 550 °C s-CO2 exposure yielded a thin (< 1 µm) Cr oxide with occasional nodules of duplex Fe oxide and Fe–Cr spinel that were up to 5 microns thick. However, tensile results from s-CO−2 exposure matched those of 550 °C thermal aging in Ar, indicating that no mechanically detrimental carburization occurred in either 316 variant after 1000 h exposure. Conversely, 750 °C s-CO2 exposure produced roughly 10 × the oxide thickness, with a more substantial Fe oxide (3–5 µm) on the majority of the surface and nodules of up to 40 µm thick. In comparison to aged samples, tensile testing of 750 °C CO2-exposed samples revealed ductility loss attributed to carburization. Projections of 316L performance in s-CO2 indicate that mechanically detrimental carburization—equal to that shown here for 750 °C, 1000 h—will likely be present after 7–14 years of service at 550 °C.
Electron emission from thick polished samples of polycrystalline molybdenum (Mo) and single crystalline 〈111〉 silver (Ag) was measured with hard x-ray photoemission spectroscopy. Six different excitation x-ray energies were used, nominally 8.0, 11.0, 13.0, 15.0, 18.0, and 21.5 keV. Survey spectra were recorded with each excitation to a kinetic energy of at most 15 keV, often capturing the entire emission range. The Mo 1s core peak was measured. Detailed LMM Auger spectra of Mo show marked increases in intensity and altered shape when x-ray energy exceeds the Mo 1s binding energy. The Mo and Ag L-shell photoelectron peaks are measured at four x-ray energies up to 18 keV showing the transition from 2p3/2 to 2s photoionization dominance.
We present the results of large scale molecular dynamics simulations aimed at understanding the origins of high friction coefficients in pure metals, and their concomitant reduction in alloys and composites. We utilize a series of targeted simulations to demonstrate that different slip mechanisms are active in the two systems, leading to differing frictional behavior. Specifically, we show that in pure metals, sliding occurs along the crystallographic slip planes, whereas in alloys shear is accommodated by grain boundaries. In pure metals, there is significant grain growth induced by the applied shear stress and the slip planes are commensurate contacts with high friction. However, the presence of dissimilar atoms in alloys suppresses grain growth and stabilizes grain boundaries, leading to low friction via grain boundary sliding. Graphic Abstract: [Figure not available: see fulltext.]
For machine vision, one of the most important operations is fast and effective object cueing or segmentation. Sandia National Labs has a long history of development and implementation of very fast and effective cueing/segmentation algorithms. This report covers the history, motivation and implementation of evolving frameworks (Sandia FOA Frameworks) upon which this long legacy of successful algorithms are built. The report describes the innovative microprocessor implementation, enabling extremely fast morphological processing, combined with a novel adaptive quantization front - end and a feature - based backend that resulted in Sandia developing fast and effective cueing in a wide variety of applications, from defect detection to SAR ATR. The report covers evolution from Sandia FOA 1.0 Framework (1995) to current Sandia FOA 4.0 Framework (2021). Requirements for the cueing algorithm for SIS - AOP (FOA_SIS - AOP) that drove the Sandia FOA 4.0 Framework development are discussed, along with information on how to use the Sandia FOA Frameworks.
Titanium hydride of varying TiH stoichiometry is used in pyrotechnic compositions. In order to yield consistent performance, manufacturing processes must be developed to ensure precise and reproducible material properties, including composition and morphology. Legacy synthesis protocols are not comprehensive nor are the required apparatuses still available. To guide the development of novel production procedures, this report reviews literature on relevant chemical reactions and diffusion events occurring at elevated temperature in the TiH(2-x)/TiOy system. Titanium hydride exposed to air spontaneously forms a passivating oxide layer. Upon heating, significant hydrogen release, which is accompanied by changes to the surface oxide layer, is noted by 375–400°C. At higher temperatures (above about 500°C) the oxide layer is reported to be essentially nonexistent as a result of oxide-layer dissolution processes and, potentially, oxide-layer reduction due to water formation. Based on the reviewed literature, we hypothesize that, by 500°C, the surface layer consists of an oxyhydride phase, which is a solid solution of oxygen in titanium hydride. We believe that hydrogen release from titanium hydride is controlled by the kinetics of molecular hydrogen desorption on the oxyhydride surface. No literature data is available for corresponding activation energies of the dynamic desorption process, and the equilibrium phase diagram of this three-component system remains largely unexplored as well. These gaps in knowledge might be addressed through coordinated computational modeling and experimental efforts.
Large volumes of data are being collected by Sandia National Laboratories as part of an active commercial-off-the-shelf (COTS) part testing and surveillance program. This user manual documents Python-based COTS Data Analytics software that has been developed for standardizing, displaying, visualizing, and analyzing the resulting COTS part testing and surveillance data. It is the objective of these software tools to streamline the analysis of COTS testing and surveillance data and improve the efficiency with which test engineers and data analytics experts can pinpoint possible performance and reliability problems in COTS parts.
The critical pitting temperature (CPT) of selective laser melted (SLM) 316 L stainless steel in 1.0 M NaCl was measured and compared with a commercial wrought alloy. Potentiostatic measurements determined a mean CPT value of 16 ± 0.7 °C, 27.5 ± 0.8 °C and 31 ± 1 °C for the wrought alloy, the SLM alloy normal to the build direction and parallel to the build direction, respectively. The lead-in pencil electrode technique was used to study the pit chemistry of the two alloys and to explain the higher CPT values observed for the SLM alloy. A lower critical current density required for passivation in a simulated pit solution was measured for the SLM alloy. Moreover, the ratio of the critical concentration to saturated concentration of dissolving metal cations was found to be higher for the SLM alloy, which was related to its different salt film properties, possibly as a result of the SLM's distinct microstructure.
Streaks in the buffer layer of wall-bounded turbulence are tracked in time to study their life cycle. Spatially and temporally resolved direct numerical simulation data are used to analyze the strong wall-parallel movements conditioned to low-speed streamwise flow. The analysis of the streaks shows that there is a clear distinction between wall-attached and detached streaks, and that the wall-attached streaks can be further categorized into streaks that are contained in the buffer layer and the ones that reach the outer region. The results reveal that streaks are born in the buffer layer, coalescing with each other to create larger streaks that are still attached to the wall. Once the streak becomes large enough, it starts to meander due to the large streamwise-to-wall-normal aspect ratio, and consequently the elongation in the streamwise direction, which makes it more difficult for the streak to be oriented strictly in the streamwise direction. While the continuous interaction of the streaks allows the superstructure to span extremely long temporal and length scales, individual streak components are relatively small and short-lived. Tall-attached streaks eventually split into wall-attached and wall-detached components. These wall-detached streaks have a strong wall-normal velocity away from the wall, similar to ejections or bursts observed in the literature. Conditionally averaging the flow fields to these split events show that the detached streak has not only a larger wall-normal velocity compared to the wall-attached counterpart, it also has a larger (less negative) streamwise velocity, similar to the velocity field at the tip of a vortex cluster.
Distributed controllers play a prominent role in electric power grid operation. The coordinated failure or malfunction of these controllers is a serious threat, where the resulting mechanisms and consequences are not yet well-known and planned against. If certain controllers are maliciously compromised by an adversary, they can be manipulated to drive the system to an unsafe state. The authors present a strategy for distributed controller defence (SDCD) for improved grid tolerance under conditions of distributed controller compromise. The work of the authors’ first formalises the roles that distributed controllers play and their control support groups using controllability analysis techniques. With these formally defined roles and groups, the authors then present defence strategies for maintaining or regaining system control during such an attack. A general control response framework is presented here for the compromise or failure of distributed controllers using the remaining, operational set. The SDCD approach is successfully demonstrated with a 7-bus system and the IEEE 118-bus system for single and coordinated distributed controller compromise; the results indicate that SDCD is able to significantly reduce system stress and mitigate compromise consequences.