The Mt. Pinatubo eruption on 15 June 1991 is often associated with surface warming in the subsequent Northern Hemisphere winter. Employing E3SMv2 with prognostic aerosol modifications, we generated an ensemble of simulations initialized on 1 June 1991 to limit the intra-ensemble variability at the time of the eruption and a more traditional ensemble representing the full range of intra-ensemble variability. For each ensemble member we generated a paired counterfactual simulation with the Pinatub forcing removed allowing for isolation of the Pinatubo impact. In general, the limited variability ensemble has greater coherence in the Pinatubo impact across ensemble members which leads to more statistically robust signals compared to the full variability ensemble. Stratospheric warming patterns from Pinatubo were approximately zonally symmetric and confined between 30°S and 50°N. Isolating localized surface temperature impacts was more difficult, but the limited variability simulation did identify a preferential region of cooling between 20°S to 50°N.
For multi-scale multi-physics applications e.g., the turbulent combustion code Pele, robust and accurate dimensionality reduction is crucial to solving problems at exascale and beyond. A recently developed technique, Co-Kurtosis based Principal Component Analysis (CoK-PCA) which leverages principal vectors of co-kurtosis, is a promising alternative to traditional PCA for complex chemical systems. To improve the effectiveness of this approach, we employ Artificial Neural Networks for reconstructing thermo-chemical scalars, species production rates, and overall heat release rates corresponding to the full state space. Our focus is on bolstering confidence in this deep learning based non-linear reconstruction through Uncertainty Quantification (UQ) and Sensitivity Analysis (SA). UQ involves quantifying uncertainties in inputs and outputs, while SA identifies influential inputs. One of the noteworthy challenges is the computational expense inherent in both endeavors. To address this, we employ the Monte Carlo methods to effectively quantify and propagate uncertainties in our reduced spaces while managing computational demands. Our research carries profound implications not only for the realm of combustion modeling but also for a broader audience in UQ. By showcasing the reliability and robustness of CoK-PCA in dimensionality reduction and deep learning predictions, we empower researchers and decision-makers to navigate complex combustion systems with greater confidence.
The motivation behind this research is the growing challenge of handling the massive amounts of data generated by modern imaging systems. Conventional digital image processing techniques are struggling to keep pace with the demands of high-resolution and high-speed imaging systems for remote sensing due to their high-power consumption and data storage requirements. We present a novel approach based on analog photonics to address this challenge. The proposed system utilizes a silicon-photonics-based image encoder positioned after image formation and initial optical-to-electrical conversion. The photonic encoder compresses image data using a passive disordered photonic structure to perform kernel-type random projections of the raw data. The compressed data is then processed by a back-end neural network, which reconstructs the original image with high fidelity (structural similarity exceeding 90%). Our proposed approach has the potential to compress images with ~ 1000X lower power consumption compared to digital approaches with data rates exceeding 1 terapixel/second.
This report summarizes the findings of a four months FY24 Advanced Science & Technology (AS&T) LDRD Quick Targeted Investigation (QTI) project focused on the exploration of supervised dimension reduction approaches based on autoencoders. Autoencoders have been extensively employed in literature for unsupervised learning tasks, however, their use for supervised regression tasks, which are common within scientific applications, has been limited. Motivated by linear dimension reduction strategies like Active Subspaces and Adaptive Basis, we explored the possibility of employing autoencoders to discover a non-linear manifold able to represent the original function in fewer dimensions. In this report, we discuss a neural network architecture and we perform a numerical campaign on several problems ranging from simple two-dimensional functions to a model problem for magnetohydrodynamics in five dimensions. In our preliminary results, we show that the proposed approach is found to be superior to linear dimension reduction strategies in representing the target function even with a single latent variable.
Numerical linear algebra (NLA) underpins huge swaths of computational science and engineering. For scientists and engineers to make the most of the DOE’s computing resources, it is essential that they have access to high-performance implementations of algorithms with best-in-class scalability and reliability. Despite this, prevailing NLA libraries have little to no support for breakthrough algorithms from the field of randomized numerical linear algebra (RandNLA) that have been developed over the past twenty years. The goal of this LDRD was to break a log-jam that had prevented broad adoption of RandNLA. Our work had two thrusts. The first was to develop RandBLAS: a trustworthy and high-performance C++ library for randomized dimension reduction (an operation widely known as sketching). The second was the development of a novel randomized algorithm for computing a challenging type of matrix decomposition known as Householder QR with column pivoting (Householder QRCP). In this one-year late-start LDRD we successfully delivered RandBLAS 1.0 and new CPU and GPU codes for Householder QRCP. RandBLAS has extensive documentation at https://randblas.readthedocs.io/en/stable/. Papers on RandBLAS and and our high-performance QRCP codes are forthcoming.
Modern computing systems are capable of exascale calculations, which are revolutionizing the development and application of high-fidelity numerical models in computational science and engineering. While these systems continue to grow in processing power, the available system memory has not increased commensurately, and electrical power consumption continues to grow. A predominant approach to limit the memory usage in large-scale applications is to exploit the abundant processing power and continually recompute many low-level simulation quantities, rather than storing them. However, this approach can adversely impact the throughput of the simulation and diminish the benefits of modern computing architectures. We present three novel contributions to reduce the memory burden while maintaining, and sometimes improving, performance in simulations based on finite element discretizations. The first contribution develops dictionary-based data compression schemes that detect and exploit the structure of the discretization, due to redundancies across the finite element mesh. While these schemes are shown to reduce memory requirements by more than 99% on meshes with large numbers of identical mesh cells, there are applications where this structure does not exist. The second contribution leverages a recently developed augmented Lagrangian optimization algorithm to enable r-adaptivity for meshes with the goal of enhancing the redundancies in the mesh. The third contribution extends these methods to patch-based linear solvers and preconditioners by compressing local matrices. Numerical results demonstrate the effectiveness of the proposed methods to detect, enhance and exploit mesh structure on a suite of examples inspired by large-scale applications.
Bugs in digital logic have led to some significant security vulnerabilities. Hardware bugs are particularly troublesome since they cannot be easily patched. Additionally, if the bug is in the root of trust, all trust built upon it can be vulnerable. Traditional testing either require a deep knowledge of the system, creative attack vectors and lots of human interaction. This is not scalable as there are very few engineers that can wear the hat of a designer, a verification engineer, and a cybersecurity expert. Hardware fuzzing is a relatively new research area in dynamic hardware testing. It has proven to be an effective method for discovering bugs, unexpected behaviors, and security vulnerabilities in software. While hardware fuzzing is new to the hardware domain, it has a strong track record in software testing. Fuzzing is a testing technique that randomly mutates the input data to uncover bugs or vulnerabilities in the design. It is especially good at finding corner cases that test engineers can not envision. Another advantage over other dynamic testing techniques is that, if done well, deep knowledge of the design is not required. Additionally, fuzzing scales well. If the system is set up correctly, it can run unsupervised for weeks if necessary. In this work, we propose using hardware fuzzing to improve the input vector generation for an information flow tracking tool. To get reasonable throughput of test vectors, an emulator is targeted as the execution platform. Efficient emulator execution has some specific requirements.
Poster to be presented at the Applied Superconductivity Conference, 2024 in Salt Lake City, Utah in September. This poster details our calculations on the propagation of ballistic fluxon solitons in long Josephson junctions.
MACCS is used by the Nuclear Regulatory Commission (NRC) and various national and international organizations for probabilistic consequence analysis of nuclear power accidents. This user guide is intended to assist analysts in understanding the MACCS/MACCS-UI User Interface (UI) model and to provide information regarding the code. This user guide version describes MACCS Version 5.0, model history, explains how to set up and execute a problem, and informs the user of the definition of various input parameters and any constraints placed on those parameters. This report is part of a series of reports documenting MACCS. Other reports include the MACCS Theory Manual, MACCS Verification Report, Technical Bases for Consequence Analyses Using MACCS, as well as documentation for preprocessor codes including SecPop, MelMACCS, and COMIDA2.
In this project we uncovered structure-function relationships of donor-acceptor co-crystals used to develop next-generation optoelectronic devices. Unraveling the photodynamics of molecular crystalline materials poses many challenges for spectroscopy due to broad, overlapping features representing numerous underlying dynamical processes. This leads researchers to make many assumptions about the dynamics of a system in choosing an appropriate kinetic fitting model. Computationally, electronic structure methods are either prohibitively expensive or underdeveloped for computing the excited state structure of molecular materials, especially states that exhibit charge transfer. Researchers must therefore perform calculations of excited electronic states using truncated models of molecular materials. Here we present a joint experimental-theoretical approach to bridging the gap between the photodynamics of a molecular material and its constituent molecules. We focus our efforts on quantifying the timescales and mechanisms of photoexcitation in donor-acceptor co-crystals and donor-acceptor dimers where the lowest-lying excited state is characterized by charge transfer from the donor to the acceptor. We employ ultrafast UV pump, UV-Vis probe transient absorption spectroscopy to unravel the time-resolved spectroscopic signatures of the photodynamics in both the crystalline material and donor-acceptor dimers in solution. We perform electronic structure and excited state dynamics calculations of the dimers to inform kinetic fitting models and assign the spectral features. The photodynamics of the crystal vs. dimer systems have many similarities, enabling unprecedented insights into the formation and evolution of charge transfer excitons in the crystalline systems.
This project develops a Bayesian optimization approach to extracting insights from Z Machine experimental data to determine if and how these insights can be used to extrapolate to a larger facility. The primary goal is to address the scientific challenge of informing how confidently experimental conditions can be predicted on a next generation facility, the design of which requires the reliable extrapolation of current high energy density technologies to regimes yet unobserved, except by costly high-fidelity computational models. Maximizing the use of presently available data and understanding how it informs future endeavors is critically important to enable transformative pulsed power and the science of extreme conditions. We explore a Bayesian optimization approach to experimental design which combines information theory, experimental data, and computational modeling to explore how information gain can be maximized.
This project focused on developing topology optimization software to design advanced metal thermal insulators. Initially, solid designs were created that matched the thermal performance of current baseline designs but were significantly heavier. To address this, cellular materials were incorporated, specifically the octet structure which is known for its high strength-to-weight ratio and thermal properties. By leveraging these cellular designs at various densities, superior thermal and mechanical performance was achieved without added weight. This novel approach enhances thermal management and structural integrity under extreme conditions, offering promising advancements for thermal protection systems.
Nitrogen-containing organic molecules play essential roles in nearly every aspect of chemistry. The importance of these species within gas-phase chemistry has increased in recent years, due to their emissions from wildfires, crude biofuels, ammonia combustion, and CO2 capture facilities. However, there is a lack of detailed relevant mechanistic and kinetic studies of N-compounds, resulting in poor representation of these pathways within chemical models. In this report, we detail our exploration of the gas-phase oxidation mechanisms of pyrrole, imidazole, pyrrolidine, methylamine, and dimethylamine, done by employing multiplexed VUV photoionization mass spectrometry coupled with KinBot, a computational tool that automatically explores multi-well potential energy surfaces. Overall, we vastly expanded the known gas-phase oxidation mechanisms and chemical kinetics of these species. This work will improve the accuracy and completeness of atmospheric and combustion chemistry models, ultimately leading to better-informed decision-makers in the energy and environmental sectors.
FLEXO (Flux-Limited Extended-MHD Ohm's Law) is a production-line multiphysics code developed at Sandia to enable more predictive modeling of target physics on pulsed-power devices. FLEXO uses an extended magnetohydrodynamics (XMHD) model which includes a generalized Ohm's law (GOL), an electron inertia term, and Hall physics. This report describes the code's numerical methods, its computational performance, and test problems of interest.
