Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

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)/TiO_y 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.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60-100 at a Mach number of 0.3. Conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This paper describes results from an optical-engine investigation of oxygenated fuel effects on ducted fuel injection (DFI) relative to conventional diesel combustion (CDC). Three fuels were tested: a baseline, non-oxygenated No. 2 emissions certification diesel (denoted CFB), and two blends containing potential renewable oxygenates. The first oxygenated blend contained 25 vol% methyl decanoate in CFB (denoted MD25), and the second contained 25 vol% tri-propylene glycol mono-methyl ether in CFB (denoted T25). Whereas DFI and fuel oxygenation primarily curtail soot emissions, intake-oxygen mole fractions of 21% and 16% were employed to explore the potential additional beneficial impact of dilution on engine-out emissions of nitrogen oxides (NOx). It was found that DFI with an oxygenated fuel can attenuate soot incandescence by ~100X (~10X from DFI and an additional ~10X from fuel oxygenation) relative to CDC with conventional diesel fuel, regardless of dilution level and without large effects on other emissions or efficiency. This breaks the soot/NOx trade-off with dilution, enabling simultaneous reductions in both soot and NOx emissions, even with conventional diesel fuel. Significant cyclic variability in soot incandescence for both CDC and DFI suggests that additional improvements in engine-out soot emissions may be possible via improved control of in-cylinder mixture formation and evolution.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Kinetic simulations of Sandia National Laboratories' Z machine are conducted to understand particle transport in the highly magnetized environment of a multi-MA accelerator. Joule heating leads to the rapid formation of electrode surface plasmas. These plasmas are implicated in reducing accelerator efficiency by diverting current away from the load [M.R. Gomez et al., Phys. Rev. Accel. Beams 20, 010401 (2017)PRABCJ2469-988810.1103/PhysRevAccelBeams.20.010401, N. Bennett et al., Phys. Rev. Accel. Beams 22, 120401 (2019)PRABCJ2469-988810.1103/PhysRevAccelBeams.22.120401]. The fully-relativistic, electromagnetic simulations presented in this paper show that particles emitted in a space-charge-limited manner, in the absence of plasma, are magnetically insulated. However, in the presence of plasma, particles are transported across the magnetic field in spite of being only weakly collisional. The simulated cross-gap currents are well-approximated by the Hall current in the generalized Ohm's law. The Hall conductivities are calculated using the simulated particle densities and energies, and the parameters that increase the Hall current are related to transmission line inductance. Analogous to the generalized Ohm's law, we extend the derivation of the magnetized diffusion coefficients to include the coupling of perpendicular components. These yield a Hall diffusion rate, which is equivalent to the empirical Bohm diffusion.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

White Dwarf (WD) stars are the most common stellar remnant in the universe. WDs usually have a hydrogen or helium atmosphere, and helium WD (called DB) spectra can be used to solve outstanding problems in stellar and galactic evolution. DB origins, which are still a mystery, must be known to solve these problems. DB masses are crucial for discriminating between different proposed DB evolutionary hypotheses. Current DB mass determination methods deliver conflicting results. The spectroscopic mass determination method relies on line broadening models that have not been validated at DB atmosphere conditions. We performed helium benchmark experiments using the White Dwarf Photosphere Experiment (WDPE) platform at Sandia National Laboratories' Z-machine that aims to study He line broadening at DB conditions. Using hydrogen/helium mixture plasmas allows investigating the importance of He Stark and van der Waals broadening simultaneously. Accurate experimental data reduction methods are essential to test these line-broadening theories. In this paper, we present data calibration methods for these benchmark He line shape experiments. We give a detailed account of data processing, spectral power calibrations, and instrument broadening measurements. Uncertainties for each data calibration step are also derived. We demonstrate that our experiments meet all benchmark experiment accuracy requirements: WDPE wavelength uncertainties are <1 Å, spectral powers can be determined to within 15%, densities are accurate at the 20% level, and instrumental broadening can be measured with 20% accuracy. Fulfilling these stringent requirements enables WDPE experimental data to provide physically meaningful conclusions about line broadening at DB conditions.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The On-Line Waste Library is a website that contains information regarding United States Department of Energy-managed high-level waste, spent nuclear fuel, and other wastes that are likely candidates for deep geologic disposal, with links to supporting documents for the data. This report provides supporting information for the data for which an already published source was not available.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Neural networks (NNs) are known as universal function approximators and can interpolate nonlinear functions between observed data points. However, when the target domain for deployment shifts from the training domain and NNs must extrapolate, the results are notoriously poor. Prior work Martinez et al. (2019) has shown that NN uncertainty estimates can be used to correct binary predictions in shifted domains without retraining the model. We hypothesize that this approach can be extended to correct real-valued time series predictions. As an exemplar, we consider two mechanical systems with nonlinear dynamics. The first system consists of a spring-mass system where the stiffness changes abruptly, and the second is a real experimental system with a frictional joint that is an open challenge for structural dynamicists to model efficiently. Our experiments will test whether 1) NN uncertainty estimates can identify when the input domain has shifted from the training domain and 2) whether the information used to calculate uncertainty estimates can be used to correct the NN’s time series predictions. While the method as proposed did not significantly improve predictions, our results did show potential for modifications that could improve models’ predictions and play a role in structural health monitoring systems that directly impact public safety.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Protecting against multi-step attacks of uncertain duration and timing forces defenders into an indefinite, always ongoing, resource-intensive response. To effectively allocate resources, a defender must be able to analyze multi-step attacks under assumption of constantly allocating resources against an uncertain stream of potentially undetected attacks. To achieve this goal, we present a novel methodology that applies a game-theoretic approach to the attack, attacker, and defender data derived from MITRE´s ATT&CK^® Framework. Time to complete attack steps is drawn from a probability distribution determined by attacker and defender strategies and capabilities. This constraints attack success parameters and enables comparing different defender resource allocation strategies. By approximating attacker-defender games as Markov processes, we represent the attacker-defender interaction, estimate the attack success parameters, determine the effects of attacker and defender strategies, and maximize opportunities for defender strategy improvements against an uncertain stream of attacks. This novel representation and analysis of multi-step attacks enables defender policy optimization and resource allocation, which we illustrate using the data from MITRE´ s APT3 ATT&CK^® Framework.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The study focuses on source region electromagnetic pulse (EMP), and provides considerations regarding the impact from an EMP due to a low-altitude nuclear denotation.

Abstract not provided.

We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poisson’s and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poisson’s problem and the linear peridynamic solid model.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report documents the progress made under the Terry Turbine Expanded Operating Band (TTEXOB) program's modeling and simulation (MODSIM) initiative at Sandia National Laboratories (SNL ). It describes the US Federal Fiscal Year 2021 (FY21) extended period-of-performance MODSIM work completed since the closure of FY20 with due reference to the Texas A&M University (TAMU) hybrid milestone 5/6 experimental program. This work, which falls under Milestone 7 of the program, provides a counterpart to the various experiments. The overall TTEXOB program and its milestone-based approach are described in the program's Summary Plan. Details of the individual milestone test plans can be found in the corresponding detailed test plan, e.g. the Milestone 3 and 4 Detailed Test Plan. SNL MODISM is conducted alongside experiments performed at TAMU, and SNL technical staff regularly consults with TAMU on the experimental program. In FY21, MELCOR code models and capabilities were exercised in two different contexts: experimental comparisons to the TAMU ZS-1 and GS-2, and stand-alone analyses of a station black-out (SBO) scenario in a generic boiling water reactor (BWR). Code to experiment comparisons met with fair success when turbine losses were well characterized as for the ZS-1 turbine. Both deterministic and Bayesian calibration processes were used to find a recommended turbine torque multiplier for ZS-1 type turbines. This process could be repeated for GS-2 type turbines if GS-2 losses were better understood. Stand-alone generic BWR SBO calculations revealed that three different modes of self-regulating turbopump behavior may be observed depending on certain modeling parameters and choices having to do with turbine nozzles. Aspects of this predicted behavior may have been observed in TAMU GS-2 experiments.

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.

We present scalable stochastic optimization approaches for improving power systems' resilience to extreme weather events. We consider both proactive redispatch and transmission line hardening as alternatives for mitigating expected load shed due to extreme weather, resulting in large-scale stochastic linear programs (LPs) and mixed-integer linear programs (MILPs). We solve these stochastic optimization problems with progressive hedging (PH), a parallel, scenario-based decomposition algorithm. Our computational experiments indicate that our proposed method for enhancing power system resilience can provide high-quality solutions efficiently. With up to 128 scenarios on a 2,000-bus network, the operations (redispatch) and investment (hardening) resilience problems can be solved in approximately 6 min and 2 h of wall-clock time, respectively. Additionally, we solve the investment problems with up to 512 scenarios, demonstrating that the approach scales very well with the number of scenarios. Moreover, the method produces high quality solutions that result in statistically significant reductions in expected load shed. Our proposed approach can be augmented to incorporate a variety of other operational and investment resilience strategies, or a combination of such strategies.

The IEA PVPS Task 13 group, experts who focus on photovoltaic performance, operation, and reliability from several leading R&D centers, universities, and industrial companies, is developing a framework for the calculation of performance loss rates of a large number of commercial and research photovoltaic (PV) power plants and their related weather data coming across various climatic zones. The general steps to calculate the performance loss rate are (i) input data cleaning and grading; (ii) data filtering; (iii) performance metric selection, corrections, and aggregation; and finally, (iv) application of a statistical modeling method to determine the performance loss rate value. In this study, several high-quality power and irradiance datasets have been shared, and the participants of the study were asked to calculate the performance loss rate of each individual system using their preferred methodologies. The data are used for benchmarking activities and to define capabilities and uncertainties of all the various methods. The combination of data filtering, metrics (performance ratio or power based), and statistical modeling methods are benchmarked in terms of (i) their deviation from the average value and (ii) their uncertainty, standard error, and confidence intervals. It was observed that careful data filtering is an essential foundation for reliable performance loss rate calculations. Furthermore, the selection of the calculation steps filter/metric/statistical method is highly dependent on one another, and the steps should not be assessed individually.

Self-assembly of iron oxide nanoparticles (IONPs) into 1D chains is appealing, because of their biocompatibility and higher mobility compared to 2D/3D assemblies while traversing the circulatory passages and blood vessels for in vivo biomedical applications. In this work, parameters such as size, concentration, composition, and magnetic field, responsible for chain formation of IONPs in a dispersion as opposed to spatially confining substrates, are examined. In particular, the monodisperse 27 nm IONPs synthesized by an extended LaMer mechanism are shown to form chains at 4 mT, which are lengthened with applied field reaching 270 nm at 2.2 T. The chain lengths are completely reversible in field. Using a combination of scattering methods and reverse Monte Carlo simulations the formation of chains is directly visualized. The visualization of real-space IONPs assemblies formed in dispersions presents a novel tool for biomedical researchers. This allows for rapid exploration of the behavior of IONPs in solution in a broad parameter space and unambiguous extraction of ​the parameters of the equilibrium structures. Additionally, it can be extended to study novel assemblies formed by more complex geometries of IONPs.

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.

In the present work, premixed combustion in a turbulent boundary layer under auto-ignitive conditions is investigated using direct numerical simulation (DNS). The turbulent inflow of the reactive DNS is obtained by temporal sampling of a corresponding inert DNS of a turbulent boundary layer at a location with Reτ= 360, where Reτ is the friction Reynolds number. The reactants of the DNS are determined by mixing the products of lean natural gas combustion and a H2/N2 fuel jet, resulting in a lean mixture of high temperature with a short ignition delay time. In the free stream the reaction front is stabilized at a streamwise location which can be predicted using the free stream velocity U∞ and the ignition delay time τig. Inside the boundary layer, combustion modifies the near-wall coherent turbulent structures considerably and turbulence results in reaction front wrinkling. The combustion modes in various regions were examined based on the results of displacement velocity, species budget and chemical explosive mode analysis (CEMA). It was indicated that flame propagation prevails in the near-wall region and auto-ignition becomes increasingly important as the wall-normal distance increases. The interactions of turbulence and combustion were studied through statistics of reaction front normal vector and strain rate tensor. It was found that the reaction front normal preferentially aligns with the most compressive strain rate in regions where the effects of heat release on the strain rate are minor and with the most extensive strain rate where its effects are significant. Negative correlations between the wall heat flux and flame quenching distance were observed. A new quenching mode, back-on quenching, was identified. It was found that the heat release rate at the wall is the highest when head-on quenching occurs and lowest when back-on quenching occurs.

Optimally-shaped electromagnetic fields have the capacity to coherently control the dynamics of quantum systems and thus offer a promising means for controlling molecular transformations relevant to chemical, biological, and materials applications. Currently, advances in this area are hindered by the prohibitive cost of the quantum dynamics simulations needed to explore the principles and possibilities of molecular control. However, the emergence of nascent quantum-computing devices suggests that efficient simulations of quantum dynamics may be on the horizon. In this article, we study how quantum computers could be employed to design optimally-shaped fields to control molecular systems. We introduce a hybrid algorithm that utilizes a quantum computer for simulating the field-induced quantum dynamics of a molecular system in polynomial time, in combination with a classical optimization approach for updating the field. Qubit encoding methods relevant for molecular control problems are described, and procedures for simulating the quantum dynamics and obtaining the simulation results are discussed. Numerical illustrations are then presented that explicitly treat paradigmatic vibrational and rotational control problems, and also consider how optimally-shaped fields could be used to elucidate the mechanisms of energy transfer in light-harvesting complexes. Resource estimates, as well as a numerical assessment of the impact of hardware noise and the prospects of near-term hardware implementations, are provided for the latter task.

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.

A discrete direct (DD) model calibration and uncertainty propagation approach is explained and demonstrated on a 4-parameter Johnson-Cook (J-C) strain-rate dependent material strength model for an aluminum alloy. The methodology’s performance is characterized in many trials involving four random realizations of strain-rate dependent material-test data curves per trial, drawn from a large synthetic population. The J-C model is calibrated to particular combinations of the data curves to obtain calibration parameter sets which are then propagated to “Can Crush” structural model predictions to produce samples of predicted response variability. These are processed with appropriate sparse-sample uncertainty quantification (UQ) methods to estimate various statistics of response with an appropriate level of conservatism. This is tested on 16 output quantities (von Mises stresses and equivalent plastic strains) and it is shown that important statistics of the true variabilities of the 16 quantities are bounded with a high success rate that is reasonably predictable and controllable. The DD approach has several advantages over other calibration-UQ approaches like Bayesian inference for capturing and utilizing the information obtained from typically small numbers of replicate experiments in model calibration situations—especially when sparse replicate functional data are involved like force–displacement curves from material tests. The DD methodology is straightforward and efficient for calibration and propagation problems involving aleatory and epistemic uncertainties in calibration experiments, models, and procedures.

In a typical optical test, a stereo camera pair is required to measure the three-dimensional motion of a test article; one camera typically only measures motions in the image plane of the camera, and measurements in the out-of-plane direction are missing. Finite element expansion techniques provide a path to estimate responses from a test at unmeasured degrees of freedom. Treating the case of a single camera as a measurement with unmeasured degrees of freedom, a finite element model is used to expand to the missing third dimension of the image data, allowing a full-field, three-dimensional measurement to be obtained from a set of images from a single camera. The key to this technique relies on the mapping of finite element deformations to image deformations, creating a set of mode shape images that are used to filter the response in the image into modal responses. These modal responses are then applied to the finite element model to estimate physical responses at all finite element model degrees of freedom. The mapping from finite element model to image is achieved using synthetic images produced by a rendering software. The technique is applied first to a synthetic deformation image, and then is validated using an experimental set of images.

Phase-based motion processing and the associated Motion Magnification that it enables has become popular not only for the striking videos that it can produce of traditionally stiff structures visualized with very large deflections, but also for its ability to pull information out of the noise floor of images so that they can be processed with more traditional optical techniques such as digital image correlation or feature tracking. While the majority of papers in the literature have utilized the Phase-based Image Processing approach as a pre-processor for more quantitative analyses, the technique itself can be used directly to extract modal parameters from an image, noting that the extracted phases are proportional to displacements in the image. Once phases are extracted, they can be fit using traditional experimental modal analysis techniques. This produces a mode “shape” where the degrees of freedom are phases instead of physical motions. These phases can be scaled to produce on-image visualizations of the mode shapes, rather than operational shapes produced by bandpass filtering. Modal filtering techniques can also be used to visualize motions from an environment on an image using the modal phases as a basis for the expansion.

Magnetic loading was used to shocklessly compress four different metals to extreme pressures. Velocimetry monitored the behavior of the material as it was loaded to a desired peak state and then decompressed back down to lower pressures. Two distinct analysis methods, including a wave profile analysis and a novel Bayesian calibration approach, were employed to estimate quantitative strength metrics associated with the loading reversal. Specifically, we report for the first time on strength estimates for tantalum, gold, platinum, and iridium under shockless compression at strain rates of ∼ 5 × 10 5/s in the pressure range of ∼ 100–400 GPa. The magnitude of the shear stresses supported by the different metals under these extreme conditions are surprisingly similar, representing a dramatic departure from ambient conditions.

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.

Abstract not provided.

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.]

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

Tritium diffusion in α-Zr containing point defects such as vacancies or self-interstitial atoms (SIAs) is simulated using molecular dynamics. Point defects rapidly aggregate to form extended defects, such as 3D nanoclusters and Frank loops. The geometry of extended defects is affected by the presence of tritium. At low temperature and in the absence of tritium, vacancies aggregate to form stacking fault pyramids. Addition of tritium at these temperatures promotes aggregation of vacancies to form 3D nanoclusters, within which the tritium concentration can be sufficiently high to suggest that these defects may serve as nucleation sites for hydride precipitation. Trapping of tritium in vacancy nanocluster reduces the calculated bulk diffusivity by an amount proportional to the vacancy concentration. At high temperature, vacancy clusters change shape to form planar basal dislocation loops, which bind tritium less strongly, leading to a sharp reduction in the fraction of trapped tritium and a corresponding increase in tritium diffusivity at high temperature. In contrast, SIAs increase tritium diffusion through α-Zr. Analysis of atomic trajectories shows that tritium does not interact directly with SIAs. In conclusion, diffusion enhancement is instead related to expansion of the lattice.

The Human Readiness Level scale complements and supplements the existing technology readiness level scale to support comprehensive and systematic evaluation of human system aspects throughout a system’s life cycle. The objective is to ensure humans can use a fielded technology or system as intended to support mission operations safely and effectively. This article defines the nine human readiness levels in the scale, explains their meaning, and illustrates their application using a helmet-mounted display example.

We present a combined experimental and theoretical investigation of the autoignition chemistry of a prototypical cyclic hydrocarbon, cyclopentane. Experiments using a high-pressure photolysis reactor coupled to time-resolved synchrotron VUV photoionization mass spectrometry directly probe the short-lived radical intermediates and products in cyclopentane oxidation reactions. We detect key peroxy radical intermediates ROO and OOQOOH, as well as several hydroperoxides, formed by second O2 addition. Automated quantum chemical calculations map out the R + O2 + O2 reaction channels and demonstrate that the detected intermediates belong to the dominant radical chain-branching pathway: ROO (+ O2) → γ-QOOH + O2 → γ-OOQOOH → products. ROO, OOQOOH, and hydroperoxide products of second-O2 addition undergo extensive dissociative ionization, making their experimental assignment challenging. We use photoionization dynamics calculations to aid in their characterization and report the absolute photoionization spectra of isomerically pure ROO and γ-OOQOOH. A global statistical fit of the observed kinetics enables reliable quantification of the time-resolved concentrations of these elusive, yet critical species, paving the way for detailed comparisons with theoretical predictions from master-equation-based models.

This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSM) from Orano (formerly Transnuclear Inc.) and provides a description of the surface characterization activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. The goal of the testing is to collect highly defensible and detailed surface deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion crack (CISCC) research. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the thermal-hydraulic-environment. Instrumentation throughout the canister, storage module, and environment will provide an extensive amount of information for the use of model validation. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer a high-fidelity quantification of the conditions experienced in a harsh yet realistic environment.

Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this study we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.

On May 26, 2021, Sandia National Laboratories (SNL) convened a diverse group of experts spanning private industry, academia, the United States military and federal government, and the national laboratories, and hosted a series of panels to gain their insight on critical emergent research and capability development needs to support national cyber strategy objectives. Two panelists of experts presented their prepared remarks, followed by open discussion from over 250 audience members.

A demonstration of 2D hole gases in GeSn/Ge heterostructures with a mobility as high as 20 000 cm² V^–1 s^–1 is given. Both the Shubnikov–de Haas oscillations and integer quantum Hall effect are observed, indicating high sample quality. The Rashba spin-orbit coupling (SOC) is investigated via magneto-transport. Further, a transition from weak localization to weak anti-localization is observed, which shows the tunability of the SOC strength by gating. The magneto-transport data are fitted to the Hikami–Larkin–Nagaoka formula. The phase-coherence and spin-relaxation times, as well as spin-splitting energy and Rashba coefficient of the k-cubic term, are extracted. Furthermore, the analysis reveals that the effects of strain and confinement potential at a high fraction of Sn suppress the Rashba SOC caused by the GeSn/Ge heterostructures.

We report a detailed study of the tunnel barriers within a single-hole GaAs/AlGaAs double quantum dot device (DQD). For quantum information applications as well as fundamental studies, careful tuning and reliable measurements of the barriers are important requirements. In order to tune a DQD device adequately into the single-hole electric dipole spin resonance regime, one has to employ a variety of techniques to cover the extended range of tunnel couplings. In this work, we demonstrate four separate techniques, based upon charge sensing, quantum transport, time-resolved pulsing, and electron dipole spin resonance spectroscopy to determine the couplings as a function of relevant gate voltages and magnetic field. Measurements were performed under conditions of both symmetric and asymmetric tunnel couplings to the leads. Good agreement was observed between different techniques when measured under the same conditions. The results indicate that even in this relatively simple circuit, the requirement to tune multiple gates and the consequences of real potential profiles result in non-intuitive dependencies of the couplings as a function of the plunger gate voltage and the magnetic field.

With evolving landscape of DC power transmission and distribution, a reliable and fast protection against faults is critical, especially for medium- and high-voltage applications. Thus, solid-state circuit breakers (SSCB), consisting of cascaded silicon carbide (SiC) junction field-effect transistors (JFET), utilize the intrinsic normally-ON characteristic along with their low ON-resistance. This approach provides an efficient and robust protection solution from detrimental short-circuit events. However, for applications that require high-voltage blocking capability, a proper number of JFETs need be connected in series to achieve the desired blocking voltage rating. Ensuring equal voltage balancing across the JFETs during the switching transitions as well as the blocking stage is critical and hence, this paper presents a novel passive balancing network for series connected JFETs for DC SSCB applications. The dynamic voltage balancing network to synchronize both the turn ON and OFF intervals is described analytically. Moreover, the static voltage balancing network is implemented to establish equal sharing of the total blocking voltage across the series connection of JFETs. The proposed dynamic and steady-state balancing networks are validated by SPICE simulation and experimental results.