In order to impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.
The domain wall-magnetic tunnel junction (DW-MTJ) is a versatile device that can simultaneously store data and perform computations. These three-terminal devices are promising for digital logic due to their nonvolatility, low-energy operation, and radiation hardness. Here, we augment the DW-MTJ logic gate with voltage-controlled magnetic anisotropy (VCMA) to improve the reliability of logical concatenation in the presence of realistic process variations. VCMA creates potential wells that allow for reliable and repeatable localization of domain walls (DWs). The DW-MTJ logic gate supports different fanouts, allowing for multiple inputs and outputs for a single device without affecting the area. We simulate a systolic array of DW-MTJ multiply-accumulate (MAC) units with 4-bit and 8-bit precision, which uses the nonvolatility of DW-MTJ logic gates to enable fine-grained pipelining and high parallelism. The DW-MTJ systolic array provides comparable throughput and efficiency to state-of-the-art CMOS systolic arrays while being radiation-hard. These results improve the feasibility of using DW-based processors, especially for extreme-environment applications such as space.
Earthquake location algorithms typically require travel time calculation. Doing this calculation in 3D, despite advances in algorithm efficiency and computational power, can still be prohibitively expensive in terms of resources and storage. Implementation of high-resolution 3D models in routine earthquake location would be a significant step forward in most of the world. Machine learning algorithms have potential to act as substitutes for travel time calculation algorithms or stored travel time tables. We investigate EikoNet - a physics informed neural network machine learning model that estimates travel times very quickly and comes with negligible memory-overhead. Specifically, we apply EikoNet to the Wasatch Fault Community Velocity Model (WFCVM), a highly detailed and complex 3D velocity model of the Salt Lake City, UT region. While routine locations in the area and studies of the 2020 Magna, UT earthquake sequence used a 1D velocity model, a 3D model may help better our understanding the structure of the major fault in the region. Our primary goal was to test the speed, memory requirements, and accuracy of EikoNet compared to a reference eikonal solver. We find that while the EikoNet is exceedingly fast and requires little memory overhead, achieving acceptable accuracy in estimated travel times is difficult and requires extensive computational resources.
Hobbs, Michael L.; Britt, Phillip F.; Hobbs, David T.; Kaneshige, Michael J.; Minette, Michael; Mintz, Jessica; Pennebaker, Frank M.; Parker, Gary R.; Rosenberg, David; Schwantes, Jon; Williams, Audrey; Pierce, Robert
Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, we systematically investigate variations in dynamic electron-ion collision frequencies ν ( ω ) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν ( ω ) . These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.
New technology for electric vehicles (EVs) must meet the requirements of higher energy usage, lower costs, and more sustainable source materials. One promising material for EV power system components is iron nitride (IN) soft magnetic composites (SMCs) because of their competitive magnetic properties and high abundance of the source materials. As part of an ongoing program at Sandia National Laboratories, this project focused on using computer modeling to optimize the prototyping process for an iron nitride SMC toroidal inductor to reach a target inductance of 600 μH. Four inductors with different combinations of wiring (26 AWG and 20 AWG) and vol% loading of iron nitride (65 vol% and 50 vol%) were fabricated at Cal Poly and characterized using an LCR meter. These inductors were also modeled using COMSOL Multiphysics™ with the Magnetic Fields module. The inductance data from the experiment and the model show that the 65 vol% IN prototypes and models agree with about 8% difference, while the 50 vol% IN samples show about a 9% difference between the prototype and the model. These results suggest that the model can predict inductance with both accuracy and precision with low confidence for the given sample size of four. An additional parameter of AC resistance is studied but the AC resistance results from the inductors and from the model generally do not agree closely, suggesting that the current model used in the project does not fully capture the mechanisms behind AC resistance of the inductor. With the focus of the project on inductance, the percent difference results of less than 9% across the four inductors that were tested increases confidence in the model’s predictive capabilities for inductance only. Using the inductance results from both the model and experiment, the final suggested inductor design is a 65 vol% core with 150 windings of 20 AWG wire that is 8 cm across and 1.5 cm tall to reach the inductance goal of 600 μH based on analysis using the optimized COMSOLTM model.
