Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

U.S. nuclear power facilities face increasing challenges in meeting dynamic security requirements caused by evolving and expanding threats while keeping costs reasonable to make nuclear energy competitive. The past approach has often included implementing security features after a facility has been designed and without attention to optimization, which can lead to cost overruns. Incorporating security in the design process can provide robust, economical, and effective physical protection systems (PPS). The purpose of this work is both to develop a framework for the integration of security into the design phase of a molten salt reactor (MSR) and show how to effectively design a PPS with a reduced staffing headcount. Specifically, this work focuses on integrating PPS design features into a developed facility layout by making minor modifications to building structures. A suite of tools, including Scribe3D©, PathTrace©, and Blender©, were used to model a hypothetical, generic domestic MSR facility. Physical protection elements such as sensors, cameras, barriers, and responders were added into the model based on defending the hypothetical MSR facility against a hypothetical design basis threat (DBT). Multiple outsider sabotage scenarios were examined, with adversary team sizes ranging from 4–8 to determine security system effectiveness. The results of this work will influence PPS designs and facility designs for U.S. domestic MSRs. This work will also demonstrate how a series of experimental and modeling capabilities across the Department of Energy (DOE) complex can impact the design and completion of security-by-design (SeBD) for small modular reactors (SMRs). The conclusions and recommendations in this document may be applicable to all SMR designs.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Describes some of the non-lithium-ion battery technologies that CEC is evaluating for long duration energy storage in CA.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract As modern neuroscience tools acquire more details about the brain, the need to move towards biological-scale neural simulations continues to grow. However, effective simulations at scale remain a challenge. Beyond just the tooling required to enable parallel execution, there is also the unique structure of the synaptic interconnectivity, which is globally sparse but has relatively high connection density and non-local interactions per neuron. There are also various practicalities to consider in high performance computing applications, such as the need for serializing neural networks to support potentially long-running simulations that require checkpoint-restart. Although acceleration on neuromorphic hardware is also a possibility, development in this space can be difficult as hardware support tends to vary between platforms and software support for larger scale models also tends to be limited. In this paper, we focus our attention on Simulation Tool for Asynchronous Cortical Streams (STACS), a spiking neural network simulator that leverages the Charm++ parallel programming framework, with the goal of supporting biological-scale simulations as well as interoperability between platforms. Central to these goals is the implementation of scalable data structures suitable for efficiently distributing a network across parallel partitions. Here, we discuss a straightforward extension of a parallel data format with a history of use in graph partitioners, which also serves as a portable intermediate representation for different neuromorphic backends. We perform scaling studies on the Summit supercomputer, examining the capabilities of STACS in terms of network build and storage, partitioning, and execution. We highlight how a suitably partitioned, spatially dependent synaptic structure introduces a communication workload well-suited to the multicast communication supported by Charm++. We evaluate the strong and weak scaling behavior for networks on the order of millions of neurons and billions of synapses, and show that STACS achieves competitive levels of parallel efficiency.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

An analysis was performed to determine whether a hydrogen jet flame impinging on a tunnel ceiling composed of multiple prestressed steel reinforced concrete box beams could result in permanent damage to the tunnel. The lower layer of the concrete box beam was modeled to determine whether heat reaches the steel reinforcing bars and whether spalling could occur. Heat transfer analysis shows that the temperature remains constant at the location of the steel reinforcing bars after 1.3 minutes of impingement and reaches a maximum of 130°C after 5 minutes. However, assuming a constant impingement for 5 minutes is an over estimation due the existing fire model which includes conservative assumptions. Explosive spalling may occur at a thin layer (~0.05 in. at 50 seconds, 0.1 in. at 5 minutes) at the bottom surface of the concrete box beam, but the steel reinforcing bars will not be exposed to the hydrogen flame.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. We rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator's statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

Shock-driven implosions with 100% deuterium (D2) gas fill compared to implosions with 50:50 nitrogen-deuterium (N2D2) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Z fill gas on implosion performance. Ion temperature (Tion) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N2D2 implosions than for the D2-only case, while the DD-neutron yield from the D2-only implosion is 7.2±0.5 times higher than from the N2D2 gas fill. The Tion enhancement for N2D2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D2-only versus N2D2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. This picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N2D2 implosions relative to the pure D2-fill implosion is determined to result from the reduced amount of D2 in the fill (fourfold effect on yield) combined with a lower fraction of the D2 fuel being hot enough to burn in the N2D2 case. The experimental yield and Tion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D2 in the N2D2 relative to the D2-only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D xrage simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased Tion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N2D2 versus D2-only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

Abstract not provided.

ARMA conference poster for a year-round grad intern from CSM

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.