Publications

Abstract not provided.

Using DTW for HPC / LDMS data clustering

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

The Canada-US Blended Cyber-Physical Exercise was a successful, first of its kind, multiorganization and multi-laboratory exercise that culminated years of complex system development and planning. The project aimed to answer three driving research questions, (1) How do cyberattacks support malicious acts leading to theft or sabotage [at a nuclear site]? (2) What are aspects of an effective combined cyber-physical response? (3) How to evaluate effectiveness of that response? Which derived the following primary objectives, 1. The May 2023 Cyber-Physical Exercise shall present a cyber-attack scenario that supports malicious acts leading to theft or sabotage. 2. The May 2023 Cyber-Physical Exercise shall define aspects of an effective combined cyber-physical response. 3. Analysis of the May 2023 Cyber-Physical Exercise shall evaluate the effectiveness of the incident response against pre-established exercise evaluation criteria. 4. Analysis of the May 2023 Cyber-Physical Exercise shall assess the effectiveness of the evaluation criteria itself. 5. Exercises shall be performed in a real-life environment. The team believes these objectives were met, and the evidence will be presented in this report. Due to the novelty of the exercise, there were several lessons learned that will be presented in this report.

Generative AI models garnered a large amount of public attention and speculation with the release of OpenAI’s chatbot, ChatGPT in November of 2022. At least two opinion camps exist – one that is excited about the possibilities these models offer for fundamental changes to human tasks, and another that is highly concerned about the power these models seem to have – especially since the release of GPT-4, which was trained on multimodal data and has ~1.7 trillion (T) parameters. We evaluated some concerns regarding these models’ power by assessing GPT 3.5 using standard, normed, and validated cognitive and personality measures. These measures come from the tradition of psychometrics in experimental psychology and have a long history of providing valuable insights and predictive distinctions in humans. For this seedling project, we developed a battery of tests that allowed us to estimate the boundaries of some of these models’ capabilities, how stable those capabilities are over a short period of time, and how they compare to humans.

In this report, we document the process related to developing a regional geologic model of a 605 x 1334 km area centered around Utah and encompassing surrounding states. This model is developed to test the effect that composition of a model has on the generation of synthetic data with the intent of using this information to improve upon full waveform moment tensor inversions. We compare observed data from three seismic events and five stations to the synthetic data generated by a preliminary model derived from a geologic framework model (GFM) developed by the USGS. The synthetic data and observed data comparisons indicate that our preliminary model performs well at smaller offset distances in the northern and central sections of the model. However, the southern stations consistently display synthetic data P- and S-wave arrival times that do not match the observed data arrival times, indicating that the velocity structure of the southern part of the model especially is inaccurate.

Abstract not provided.

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Reτ=505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an N-way rank-R parametric tensor M ∈ (0, ∞)I×.....×I generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately IR2 log22(I) non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material's mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. To construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer's glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Traditionally, electric grid planning seeks to maintain safe, reliable, efficient, and affordable service for current and future customers. As policies, expectations of the energy system, and the threat landscape evolve, additional objectives for power system planners are emerging, including decarbonization, resilience, and equity. Renewable and clean energy goals, especially in the context of deep decarbonization strategies, are changing the mix of resources on the electric grid and prompting new considerations for grid architecture. The increased frequency and severity of extreme weather events over the last two decades, coupled with cybersecurity concerns, have elevated resilience as a key system need. More recently, there has been greater focus on equity and energy justice in grid planning to ensure that disadvantaged communities are not adversely affected by grid modernization and have equal access to its benefits. In response, new thinking around multi-objective decision planning is exploring improvements in grid planning processes to better integrate approaches to meet decarbonization, resilience, and equity objectives. To provide a foundation for this work, a series of white papers was produced to summarize these emerging objectives.

This SAND report collects the results from the LDRD project “SHAZAM”, which aimed to push the limits of performance for self-healing, self-assembling power systems whose sectionalizing and load-control agents rely on local measurements only (i.e., only what they can measure at their own terminals, with no data sharing between agents). This work includes self-networking microgrids. The key objectives of this work were a) to demonstrate how high the performance of local-measurement-only self-assembling power systems can be; and b) to solve certain technical problems associated with such systems, such as their inability to prevent the accidental formation of closed loops and their tendency to thermally overload some conductors. “SHAZAM” investigators a) demonstrated that the performance of such systems can be surprisingly high, b) demonstrated that such systems are quite robust to all kinds of variations, and c) developed and demonstrated solutions to several key challenges associated with this type of system.

The modern global economy relies heavily on carbon-based products that are derived from petroleum, which presents sustainability, resource management, and greenhouse gas exacerbated climate change challenges. Due to these challenges, there is the need for a global industrial transition towards green and sustainable production. Microbial production of valuable chemicals from renewable biomass represents one promising route. However, high-volume low-value products such as commodity chemicals are still difficult to make profitable. One fundamental bottleneck is a waste of more than 1/3 of the feedstock carbon as CO₂ in the fermentation process. Here the project focuses on fundamentally reconfiguring the metabolism to reduce CO₂ loss in central metabolic pathways thereby also improving bioproduct yields. Here we present technologies to prevent CO₂ loss and balance reducing equivalents within the cell to enable complete conversion of glucose from renewable feedstocks into bioproducts.

A discussion of systems like ChatGPT, what their legal issues may be, how they are affecting society, and what the ethical considerations of their existence and use are.

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

We combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S(3P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S(3P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O(3P) reaction.

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

This paper explores the utility of organizational system modeling frameworks to provide valuable insight into information flows within organizations and subsequently the opportunities for increasing resilience against disinformation campaigns targeting the system's ability to utilize information within its decision making. Disinformation is a growing challenge for many organizations and in recent years has created delay in decision making. Here the paper has utilized the viable systems model (VSM) to characterize organizational systems and used this approach to outline potential subsystem requirements to promote resilience of the system. The results of this paper can support the development of simulations and models considering the human elements within the system as well as support the development of quantitative measures of resilience.

The US Department of Energy (DOE) is investigating the use of different materials that could be used to fill the void space inside a dual-purpose canister (DPC) loaded with spent nuclear fuel (SNF) just before it is emplaced in a deep geologic repository. The purpose of adding filler material is to maintain subcritical conditions in the repository during the postclosure period, which can span up to 1,000,000 years. Several types of materials have been proposed, including metals, cements, particulates, and glass. Part of this investigation addresses how the presence of filler material inside a DPC will affect the performance of the repository with respect to the repository features; the consequences of events that may occur; and the multiple thermal, hydrologic, chemical, and mechanical processes that may occur in a deep geologic repository over long timescales. This report describes some of the filler materials that have been proposed and studied; identifies 11 features, 6 events, and 25 processes that may be affected by the presence of filler materials; and discusses the effects that may require consideration for each feature, event, or process. The results of this study can be used to direct appropriate research and to develop suitable models if the DOE decides to use fillers to maintain subcritical conditions in DPCs used to dispose of SNF.

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.