Lehoucq, Rich; Mckinley, Scott A.; Miles, Christopher E.; Ding, Fangyuan
Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. These snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.
This slide is intended to serve as an informational display located on the machine high bay for visitors on Z tours. This display reflects some of our updated research capabilities that Z performs to support our nuclear deterrence mission.
This manuscript presents a complete framework for the development and verification of physics-informed neural networks with application to the alternating-current power flow (ACPF) equations. Physics-informed neural networks (PINN)s have received considerable interest within power systems communities for their ability to harness underlying physical equations to produce simple neural network architectures that achieve high accuracy using limited training data. The methodology developed in this work builds on existing methods and explores new important aspects around the implementation of PINNs including: (i) obtaining operationally relevant training data, (ii) efficiently training PINNs and using pruning techniques to reduce their complexity, and (iii) globally verifying the worst-case predictions given known physical constraints. The methodology is applied to the IEEE-14 and 118 bus systems where PINNs show substantially improved accuracy in a data-limited setting and attain better guarantees with respect to worst-case predictions.
This data documentation report describes geologic and hydrologic laboratory analysis and data collected in support of site characterization of the Physical Experiment 1 (PE1) testbed, Aqueduct Mesa, Nevada. The documentation includes a summary of laboratory tests performed, discussion of sample selection for assessing heterogeneity of various testbed properties, methods, and results per data type.
This is a presentation for a short course on accelerators at the 2024 International Power Modulator and High Voltage Conference (Indianapolis, IN, May 28-Jun 1, 2024), and largely consists of already-published material.
Crystal plasticity finite element method (CPFEM) has been an integrated computational materials engineering (ICME) workhorse to study materials behaviors and structure-property relationships for the last few decades. These relations are mappings from the microstructure space to the materials properties space. Due to the stochastic and random nature of microstructures, there is always some uncertainty associated with materials properties, for example, in homogenized stress-strain curves. For critical applications with strong reliability needs, it is often desirable to quantify the microstructure-induced uncertainty in the context of structure-property relationships. However, this uncertainty quantification (UQ) problem often incurs a large computational cost because many statistically equivalent representative volume elements (SERVEs) are needed. In this article, we apply a multi-level Monte Carlo (MLMC) method to CPFEM to study the uncertainty in stress-strain curves, given an ensemble of SERVEs at multiple mesh resolutions. By using the information at coarse meshes, we show that it is possible to approximate the response at fine meshes with a much reduced computational cost. We focus on problems where the model output is multi-dimensional, which requires us to track multiple quantities of interest (QoIs) at the same time. Our numerical results show that MLMC can accelerate UQ tasks around 2.23×, compared to the classical Monte Carlo (MC) method, which is widely known as ensemble average in the CPFEM literature.
Photonic Doppler Velocimetry (PDV) is a fiber-based measurement amenable to a wide range of experimental conditions. Interference between two optical signals—one Doppler shifted and the other not—is the essential principle in these measurements. A confluence of commercial technologies, largely driven by the telecommunication industry, makes PDV particularly convenient at near-infrared wavelengths. This discussion considers how measurement time scales of interest relate to the design, operation, and analysis of a PDV measurement, starting from the steady state through nanosecond resolution. Benefits and outstanding challenges of PDV are summarized, with comparisons to related diagnostics.
A thermally driven, micrometer-scale switch technology has been created that utilizes the ErH3/Er2O3 materials system. The technology is comprised of novel thin film switches, interconnects, on-board micro-scale heaters for passive thermal environment sensing, and on-board micro-scale heaters for individualized switch actuation. Switches undergo a thermodynamically stable reduction/oxidation reaction leading to a multi-decade (>11 orders) change in resistance. The resistance contrast remains after cooling to room temperature, making them suitable as thermal fuses. An activation energy of 290 kJ/mol was calculated for the switch reaction, and a thermos-kinetic model was employed to determine switch times of 120 ms at 560 °C with the potential to scale to 1 ms at 680 °C.
The apparent velocity measured by an interferometric surface velocimeter is a function of both the surface velocity and the time derivative of the refractive index along the measurement path. We employed this dual sensitivity to simultaneously measure km/s surface velocities and 1018 cm−3 average plasma densities with combined VISAR (velocity interferometer system for any reflector) and PDV (photonic Doppler velocimetry) measurements in experiments performed on the Z Pulsed Power Facility. We detail the governing equations, associated assumptions, and analysis specifics and show that the surface velocity can be extracted without knowledge of the specific plasma density profile.
Hydrogen energy storage can be used to achieve goals of national energy security, renewable energy integration, and grid resilience. Adapting underground natural gas storage facility (UNGSF) infrastructure for underground hydrogen storage (UHS) is one method of storing large quantities of hydrogen that has already largely been proven to work for natural gas. There are currently some underground salt caverns in the United States that are being used for hydrogen storage by commercial entities, but it is still a fairly new concept in that it has not been widely deployed nor has it been done with other geologic formations like depleted hydrocarbon reservoirs. Assessments of UHS systems can help identify and evaluate risks to people both working within the facility and residing nearby. This report provides example risk assessment methodologies and analyses for generic wellhead and processing facility configurations, specifically in the context of the risks of unintentional hydrogen releases into the air. Assessment of the hydrogen containment in the subsurface is also critically important for a safety assessment for a UHS facility, but those geomechanical assessments are not included in this report.
Sandia National Laboratories partnered with a multi-disciplinary group of subject matter experts to evaluate a stratigraphic geothermal resource in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The stratigraphic reservoir in northern Steptoe Valley was previously discovered during oil and gas exploration. Subsequent studies, such as the Nevada Play Fairway Analysis, included data which further highlighted potential resource targets in the basin. Geophysical surveys, complemented with refined geologic mapping and geochemical sampling, were deployed to further characterize the resource. The resulting 3D geologic interpretation, conceptual model refinements, and reservoir simulations suggest that a >100MWe power-capable reservoir is likely economically accessible using conventional well placement and stimulation techniques in the Paleozoic carbonates of the deep/central basin of northern Steptoe Valley. Additional geophysical characterization and exploration drilling efforts are recommended to calibrate interpretation and determine where/how to potentially develop the northern Steptoe resource. The geophysical tools, interpretations, lessons learned, and public data generated by this study establish an exploration methodology to inform decisions for characterization and development of northern Steptoe Valley and other stratigraphic geothermal reservoirs in the western U.S.
The Single Volume Scatter Camera (SVSC) Collaboration aims to develop portable neutron imaging systems for a variety of applications in nuclear non-proliferation. Conventional double-scatter neutron imagers are composed of several separate detector volumes organized in at least two planes. A neutron must scatter in two of these detector volumes for its initial trajectory to be reconstructed. As such, these systems typically have a large footprint and poor geometric efficiency. We report on the design and characterization of a prototype monolithic neutron scatter camera that is intended to significantly improve upon the geometrical shortcomings of conventional neutron cameras. The detector consists of a 50 mm×56 mm× 60 mm monolithic block of EJ-204 plastic scintillator instrumented on two faces with arrays of 64 Hamamatsu S13360-6075PE silicon photomultipliers (SiPMs). The electronic crosstalk is limited to < 5% between adjacent channels and < 0.1% between all other channel pairs. SiPMs introduce a significantly elevated dark count rate over PMTs, as well as correlated noise from after-pulsing and optical crosstalk. In this article, we characterize the dark count rate and optical crosstalk and present a modified event reconstruction likelihood function that accounts for them. We find that the average dark count rate per SiPM is 4.3 MHz with a standard deviation of 1.5 MHz among devices. The analysis method we employ to measure internal optical crosstalk also naturally yields the mean and width of the single-electron pulse height. We calculate separate contributions to the width of the single-electron pulse-height from electronic noise and avalanche fluctuations. We demonstrate a timing resolution for a single-photon pulse to be (128 ± 4) ps. Finally, coincidence analysis is employed to measure external (pixel-to-pixel) optical crosstalk. We present a map of the average external crosstalk probability between 2×4 groups of SiPMs, as well as the in-situ timing characteristics extracted from the coincidence analysis. Further work is needed to characterize the performance of the camera at reconstructing single- and double-site interactions, as well as image reconstruction.
Artificial intelligence (AI) and machine learning (ML) are near-ubiquitous in day-to-day life; from cars with automated driver-assistance, recommender systems, generative content platforms, and large language chatbots. Implementing AI as a tool for international safeguards could significantly decrease the burden on safeguards inspectors and nuclear facility operators. The use of AI would allow inspectors to complete their in-field activities quicker, while identifying patterns and anomalies and freeing inspectors to focus on the uniquely human component of inspections. Sandia National Laboratories has spent the past two and a half years developing on-device machine learning to develop both a digital and robotic assistant. This combined platform, which we term INSPECTA, has numerous on-device machine learning capabilities that have been demonstrated at the laboratory scale. This work describes early successes implementing AI/ML capabilities to reduce the burden of tedious inspector tasks such as seal examination, information recall, note taking, and more.
Somoye, Idris O.; Plusquellic, Jim; Mannos, Tom J.; Dziki, Brian
Recent evaluations of counter-based periodic testing strategies for fault detection in Microprocessor(μP) have shown that only a small set of counters is needed to provide complete coverage of severe faults. Severe faults are defined as faults that leak sensitive information, e.g., an encryption key on the output of a serial port. Alternatively, fault detection can be accomplished by executing instructions that periodically test the control and functional units of the μP. In this paper, we propose a fault detection method that utilizes an ’engineered’ executable program combined with a small set of strategically placed counters in pursuit of a hardware Periodic Built-In-Self-Test (PBIST). We analyze two distinct methods for generating such a binary; the first uses an Automatic Test Generation Pattern (ATPG)-based methodology, and the second uses a process whereby existing counter-based node-monitoring infrastructure is utilized. We show that complete fault coverage of all leakage faults is possible using relatively small binaries with low latency to fault detection and by utilizing only a few strategically placed counters in the μP.
In this paper we extend the DGiT multirate framework, developed in Connors and Sockwell (2022) for scalar transmission problems, to a solid–solid interaction (SSI) problem involving two coupled elastic solids and a coupled air–sea model with the rotating, thermal shallow water equations. In so doing we aim to demonstrate the broad applicability of the mathematical theory and governing principles established in Connors and Sockwell (2022) to coupled problems characterized by subproblems evolving at different temporal scales. Multirate time integration algorithms employing different time steps, optimized for the dynamics of each subproblem, can significantly improve simulation efficiency for such coupled problems. However, development of multirate algorithms is a highly non-trivial task due to the coupling, which can impact accuracy, stability or other desired properties such as preservation of system invariants. DGiT provides a general template for multirate time integration that can achieve these properties. To elucidate the manner in which DGiT accomplishes this task, we fully detail each step in the application of the framework to the SSI and air–sea coupled problems. Numerical examples illustrate key properties of the resulting multirate schemes for both problems.
Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.
