A collection of MATLAB functions and class definitions called System Workflow Tools (SWFT) are available to semi-automate steps in the simulation process. Some of these steps are often simple and routine for smaller finite element models, but if done directly by an analyst can quickly become labor intensive, cumbersome, and error prone for larger, system level models. Some of SWFT’s capabilities demonstrated in this report includes writing Sierra input decks and processing Quantities of Interest (QOI) from results files. SWFT also writes scripts in order to utilize other software programs such as Cubit (separating system level CAD into subassemblies and components, creating nodesets and sidesets), DAKOTA (ensemble management), and ParaView (contour plots and animations). Detailed commands and workflows from mesh generation to report generation are provided as examples for analysts to utilize SWFT capabilities.
This report presents our work to model the workloads of a linear electromagnetic application based on the method of moments in the frequency domain to effectively load balance the matrix assembly. This application is particularly challenging to load balance due to its lack of persistent iterative behavior, its operation under tight memory constraint (where the matrix may fill 80% of memory on each node), and the algorithmic complexity of the computational method. This report describes the first step in our work to apply an inspector-executor approach for load balancing workloads where key parameters are exposed during the inspector phase and a pre-trained model is applied to predict relative task weights for the load balancer.
Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. We found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.
Plasma distribution in 3D space is heavily influenced by complex surfaces and the coupling interactions between plasma properties and interfacing material properties. For example, guided streamers that transition to surface ionization waves (SIWs) and propagate over structured dielectrics experience field enhancements that can lead to localized increases in ionization rates and complex 3D configurations that are difficult to analyze. Investigating these configurations requires techniques than can provide a more complete 3D picture. To help address this capability gap, a tomographic optical emission spectroscopy (tomo-OES) diagnostic system has been developed at Sandia National Laboratories that can resolve SIWs. The system includes four intensified cameras that measure the angular projections of the plasma light emission through bandpass filters. A dot calibration target co-registers each angular projection to the same voxel grid and an algebraic reconstruction technique (ART) recovers the light intensity at each voxel. An atmospheric pressure plasma jet (APPJ), provided by Peter Bruggeman, has been investigated and representative results are shown in Figure 1. Here, a bandpass filter was used to isolate emission from the N2 second positive system (SPS) at 337.1 nm to capture the transition of the streamer to SIW on a planar dielectric surface (relative permittivity 3.3) located 3 mm below the APPJ [3]. The surface wave velocity was 3.5x104 (m/s), consistent with measurements made by Steven Shannon. Characterization of this APPJ will support the group effort of standing up a reproducible APPJ across institutions for applications such as liquid treatment, catalysis, and plasma aided combustion. Future work will investigate non-planar surfaces and eventually develop tomographic laser-induced fluorescence (tomo-LIF) approaches.
Accurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.
In the dynamic landscape of Operational Technology (OT), and specifically the emerging landscape for Advanced Reactors, the establishment of trust between digital assets emerges as a challenge for cybersecurity modernization. This report reviews existing approaches to authentication in Enterprise environments, and proposed methods for authentication in OT, and analyzes each for its applicability to future Advanced Reactor digital networks. Principles of authentication ranging from underlying cryptographic mechanisms to trust authorities are evaluated through the lens of OT. These facets emphasize the importance of mutual authentication in real-time environments, enabling a paradigm shift from the current approach of strong boundaries to a more malleable network that allows for flexible operation. This work finds that there is a need for evaluation and decision making by industry stakeholders, but current technologies and approaches can be adapted to fit needs and risk tolerances.
Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. The effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.
This project has identified opportunities to bring further reductions in the mass and cost of modern wind turbine blades through the use of alternative material systems and manufacturing processes. The fiber reinforced polymer material systems currently used by the wind industry have stagnated as the technology continues to mature and as a means to reduce risk while introducing new products with continually increasing blade lengths. However, as blade lengths continue to increase, the challenge of controlling blade mass becomes even more critical to enabling the associated levelized cost of energy reductions. Stiffer and stronger reinforcement fibers can help to resolve the challenges of meeting the loading demands while limiting the increase in weight, but these materials are substantially more expensive than the traditional E-glass fiber systems. One goal of this project and associated work is to identify pathways that improve the cost-effectiveness of carbon fiber such that it is the reinforcement of choice in the primary structural elements of wind blades. The use of heavy-tow textile carbon fiber material systems has been shown to reduce the blade mass by 30-31% when used in the spar cap and by up to 7% when used in edgewise reinforcement. A pultrusion cost model was developed to enable a material cost comparison that includes an accurate estimate of the intermediate manufacturing step of pultrusion for the carbon fiber composite. Material cost reductions were revealed in most cases for the heavy-tow textile carbon fiber compared to infused fiberglass. The use of carbon fiber in the edgewise reinforcement produced the most notable material cost reduction of 33% for the heavy-tow textile carbon fiber. The mass and cost savings observed when using carbon fiber in edgewise reinforcement demonstrate a clear opportunity of this design approach. A carbon fiber conversion cost model was expanded to include a characterization of manufacturing costs when using advanced conversion processes with atmospheric plasma oxidation. This manufacturing approach was estimated to reduce the cost of carbon fiber material systems by greater than 10% and can be used with textile carbon systems or traditional carbon fiber precursors. The pultrusion cost model was also used to assess the opportunity for using pultruded fiberglass in wind blades, studying conventional E-glass fiber reinforcement. When using pultruded fiberglass as the spar cap material for two design classifications, the blade weight was reduced by 6% and 9% compared to infused fiberglass. However, due to the relatively large share of the pultrusion manufacturing cost compared to fiber cost, the spar cap material cost increased by 12% and 7%. When considering the system benefits of reduced blade mass and potentially lower blade manufacturing costs for pultruded composites, there may be opportunity for pultruded E-glass in wind blade spar caps, but further studies are needed. There is a clearer outcome for using pultruded fiberglass in the edgewise reinforcement where it resulted in a blade mass reduction of 2% and associated reinforcement material cost reduction of 1% compared to infused E-glass. The use of higher performing glass fibers, such as S-glass and H-glass systems, will produce greater mass savings but a study is needed to assess the cost implications for these more expensive systems. The most likely opportunity for these high-performance glass fibers is in the edgewise reinforcement, where the increased strength will reduce the damage accumulation of this fatigue-driven component. The blade design assessments in this project characterize the controlling material properties for the primary structural components in the flapwise and edgewise directions for modern wind blades. The observed trends with low and high wind speed turbine classifications for carbon and glass fiber reinforced polymer systems help to identify where cost reductions are needed, and where improvements in mechanical properties would help to reduce the material demands.
The ability to track the concentrations of specific molecules in the body in real time would significantly improve our ability to study, monitor, and respond to diseases. To achieve this, we require sensors that can withstand the complex environment inside the body. Electrochemical aptamer-based sensors are particularly promising for in vivo sensing, as they are among the only generalizable sensing technologies that can achieve real-time molecular monitoring directly in blood and the living body. In this project, we first focused on extending the application space of aptamer sensors to support minimally-invasive wearable measurements. To achieve this, we developed individually-addressable sensors with commercial off-the-shelf microneedles. We demonstrated sensor function in buffer, blood, and porcine skin (a common proxy for human skin). In addition to the applied sensing project, we also worked to improve fundamental understanding of the aptamer sensing platform and how it responds to biomolecular interferents. Specifically, we explored the interfacial dynamics of biofouling – a process impacting sensors placed in complex fluids, such as blood.
As a follow-up to our more comprehensive report on Adversarial Machine Learning (AML), here we provide demonstrations of AML attacks against the Limbo image database of UF6 cylinders in a variety of orientations and amongst a variety of distractor images. We demonstrate the Carlini & Wagner AML attack against a subset of Limbo images, with 100% attack success rate; meaning all attacked images were misclassified by a highly accurate trained model, yet the image changes were imperceptible to the human eye. We also demonstrate successful attacks against segmented images (images with more than one targeted object). Finally, we demonstrated the Fast Fourier Transform countermeasure that can be used to detect AML attacks on images. The intent of this and our previous report is to inform the IAEA and stakeholders of both the promise of machine learning, which could greatly improve the efficiency of surveillance monitoring, but also of the real threat of AML and potential defenses.
In this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N 2 + , and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N 2 + and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N 2 + that spatially converged to form the ‘plasma bullet’ or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.
Accurate targeting of radioisotope classifiers and estimators requires an understanding of the target problem space. In order to facilitate clear communication on expected model behavior and performance between practitioners and stakeholders on their problems, this questionnaire was created. Stakeholder responses form the basis of a trained model as well as the start of usage requirements for the model as it is integrated with analysis processes or detection systems. This questionnaire may also be useful to machine learning practitioners and gamma spectroscopists developing new algorithms as a starting point for characterizing their problem space, especially if they are using PyRIID.
Tritium exhibits unique environmental behavior because of its potential interactions with water and organic substances. Modeling the environmental consequences of tritium releases can be relatively complex and thus an evaluation of MACCS is needed to understand what updates, if any, are needed in MACCS to account for the behavior of tritium. We examine documented tritium releases and previous benchmarking assessments to perform a model intercomparison between MACCS and state-of-practice tritium-specific codes UFOTRI and ETMOD to quantify the difference between MACCS and state of practice models for assessing tritium consequences. Additionally, information to assist an analyst in judging whether a postulated tritium release is likely to lead to significant doses is provided.
There is an increasing aspiration to utilize machine learning (ML) for various tasks of relevance to national security. ML models have thus far been mostly applied to tasks and domains that, while impactful, have sufficient volume of data. For predictive tasks of national security relevance, ML models of great capacity (ability to approximate nonlinear trends in input-output maps) are often needed to capture the complex underlying physics. However, scientific problems of relevance to national security are often accompanied by various sources of sparse and/or incomplete data, including experiments and simulations, across different regimes of operation, of varying degrees of fidelity, and include noise with different characteristics and/or intensity. State-of-the-art ML models, despite exhibiting superior performance on the task and domain they were trained on, may suffer detrimental loss in performance in such sparse data environments. This report summarizes the results of the Laboratory Directed Research and Development project entitled Trust-Enhancing Probabilistic Transfer Learning for Sparse and Noisy Data Environments. The objective of the project was to develop a new transfer learning (TL) framework that aims to adaptively blend the data across different sources in tackling one task of interest, resulting in enhanced trustworthiness of ML models for mission- and safety-critical systems. The proposed framework determines when it is worth applying TL and how much knowledge is to be transferred, despite uncontrollable uncertainties. The framework accomplishes this by leveraging concepts and techniques from the fields of Bayesian inverse modeling and uncertainty quantification, relying on strong mathematical foundations of probability and measure theories to devise new uncertainty-aware TL workflows.
