This document is aimed at providing guidance to the National Nuclear Security Administration’s (NNSA) Office of International Nuclear Security’s (INS) country and regional teams for implementing effective physical protection systems (PPSs) for nuclear power plants (NPPs) to prevent the radiological consequences of sabotage. This recommendation document includes input from the Physical Protection Functional Team (PPFT), the Response Functional Team (RFT), and the Sabotage Functional Team (SFT) under INS. Specifically, this document provides insights into increasing and sustaining physical protection capabilities at INS partner countries’ NPP sites. Nuclear power plants should consider that the intent of this document is to provide a historical context as well as technologies and methodologies that may be applied to improve physical protection capabilities. It also refers to relevant guidance from the International Atomic Energy Agency (IAEA) and the U.S. Nuclear Regulatory Commission (NRC).
Optimization of automated data collection is gaining increased interest for the purposes of enabling closed-loop self-correcting systems that inherently maximize operational efficiencies and reduce waste. Many data collection systems have several variables which influence data accuracy or consistency and which can require frequent user interaction to be monitored and maintained. Operating upon a Robo-MET.3D™ automated mechanical serial-sectioning system, a run-to-run control algorithm has been developed to accelerate data collection and reduce data inconsistency. Using historical data amassed over a decade of experiments, a linear regression model of the deterministic system dynamics is created and used to employ a run-to-run control algorithm that optimizes selected system inputs to reduce operator intervention and increase efficacy while reducing variance of system output.
We show how to sample uniformly within the three-sided region bounded by a circle, a radial ray, and a tangent, called a “chock.” By dividing a 2D planar rectangle into a background grid, and subtracting Poisson disks from grid squares, we are able to represent the available region for samples exactly using triangles and chocks. Uniform random samples are generated from chock areas precisely without rejection sampling. This provides the first implemented algorithm for precise maximal Poisson-disk sampling in deterministic linear time. We prove O(n · M(b) log b), where n is the number of samples, b is the bits of numerical precision and M is the cost of multiplication. Prior methods have higher time complexity, take expected time, are non-maximal, and/or are not Poisson-disk distributions in the most precise mathematical sense. We fill this theoretical lacuna.
Jang, Taejin; Mishra, Lubhani; Roberts, Scott A.; Planden, Brady; Subramaniam, Akshay; Uppaluri, Maitri; Linder, David; Gururajan, Mogadalai P.; Zhang, Ji G.; Subramanian, Venkat R.
Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. The article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.
In June of 2022, dust samples were collected from the surface of an in-service spent nuclear fuel dry storage canister during an inspection at an Independent Spent Fuel Storage Installation. The site is anonymous but is a near-marine or brackish water east coast location referred to here as "Site C". The purpose of the sampling was to assess the composition and abundance of the soluble salts present on the canister surface, information that provides a metric for potential corrosion risks. Following collection, the samples were delivered to Sandia National Laboratories for analysis. At Sandia, the soluble salts were leached from the dust and quantified by ion chromatography. In addition, subsamples of the dust were taken for scanning electron microscopy to determine the particle sizes, morphology, and mineralogy of the dust and salts. The results of those analyses are presented in this report.
In this work we employ an encoder–decoder convolutional neural network to predict the failure locations of porous metal tension specimens based only on their initial porosities. The process we model is complex, with a progression from initial void nucleation, to saturation, and ultimately failure. The objective of predicting failure locations presents an extreme case of class imbalance since most of the material in the specimens does not fail. In response to this challenge, we develop and demonstrate the effectiveness of data- and loss-based regularization methods. Since there is considerable sensitivity of the failure location to the particular configuration of voids, we also use variational inference to provide uncertainties for the neural network predictions. We connect the deterministic and Bayesian convolutional neural network formulations to explain how variational inference regularizes the training and predictions. We demonstrate that the resulting predicted variances are effective in ranking the locations that are most likely to fail in any given specimen.
Thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio-posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens-is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300-430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.
This report is a companion document to a series of six white papers, prepared jointly by the Proliferation Resistance and Physical Protection Working Group (PRPPWG) and the six System Steering Committees (SSCs) and provisional System Steering Committees (pSSCs). This publication is an update to a similar series published in 2011 presenting crosscutting Proliferation Resistance & Physical Protection (PR&PP) characteristics for the six systems selected by the Generation IV International Forum (GIF) for further research and development, namely: the Lead-cooled Fast Reactor (LFR), the Sodium-cooled fast Reactor (SFR), the Very high temperature reactor (VHTR), the gas-cooled fast reactor (GFR), the Molten salt reactor (MSR) and the Supercritical water–cooled reactor (SCWR).
For decades, Arctic temperatures have increased twice as fast as average global temperatures. As a first step toward quantifying parametric uncertainty in Arctic climate, we performed a variance-based global sensitivity analysis (GSA) using a fully coupled, ultra-low resolution (ULR) configuration of version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SMv1). Specifically, we quantified the sensitivity of six quantities of interests (QOIs), which characterize changes in Arctic climate over a 75 year period, to uncertainties in nine model parameters spanning the sea ice, atmosphere, and ocean components of E3SMv1. Sensitivity indices for each QOI were computed with a Gaussian process emulator using 139 random realizations of the random parameters and fixed preindustrial forcing. Uncertainties in the atmospheric parameters in the Cloud Layers Unified by Binormals (CLUBB) scheme were found to have the most impact on sea ice status and the larger Arctic climate. Our results demonstrate the importance of conducting sensitivity analyses with fully coupled climate models. The ULR configuration makes such studies computationally feasible today due to its low computational cost. When advances in computational power and modeling algorithms enable the tractable use of higher-resolution models, our results will provide a baseline that can quantify the impact of model resolution on the accuracy of sensitivity indices. Moreover, the confidence intervals provided by our study, which we used to quantify the impact of the number of model evaluations on the accuracy of sensitivity estimates, have the potential to inform the computational resources needed for future sensitivity studies.
We present a polynomial preconditioner for solving large systems of linear equations. The polynomial is derived from the minimum residual polynomial (the GMRES polynomial) and is more straightforward to compute and implement than many previous polynomial preconditioners. Our current implementation of this polynomial using its roots is naturally more stable than previous methods of computing the same polynomial. We implement further stability control using added roots, and this allows for high degree polynomials. We discuss the effectiveness and challenges of root-adding and give an additional check for stability. In this article, we study the polynomial preconditioner applied to GMRES; however it could be used with any Krylov solver. This polynomial preconditioning algorithm can dramatically improve convergence for some problems, especially for difficult problems, and can reduce dot products by an even greater margin.
The U.S. Department of Energy has a broad mission to ensure the nation’s security by addressing ongoing environmental challenges, developing novel energy production technologies, and mitigating nuclear security concerns. These efforts benefit from developing more effective semiconductor materials and devices. The research discussed in this report aims to further the development of these new electronics by improving implementations of vertical device structures, which offer theoretically better performance compared to lateral structures. Specifically, it reviews the background and development of a new and easy-to-use program for data processing and analysis from experiments run on vertical gallium nitride power devices. Current vertical structures fall short in critical performance metrics, such as On-Resistance and Breakdown Voltage, due to poor management of the device’s electric field. Junction termination extensions, or JTEs, are a field management technique that increases a device’s resilience to expected failure modes by controlling its surface electric fields. Each JTE design requires significant experimental validation, consisting of hundreds of tested devices, each outputting an enormous set of data points. The new software includes multiple methods of filtering large sets of device testing data to identify and remove flawed devices, automate device analyses, and provide statistical breakdowns for four metrics used to define the effectiveness of the JTE and the current carrying capacity of the device. These results can then be related back to fabrication techniques and device design, and in the future, be combined with processing simulations to further improve our understanding of the JTE mechanisms.
We present an elementary and self-contained proof that there are exactly four symmetry classes of the elasticity tensor in two dimensions: oblique, rectangular, square, and isotropic. In two dimensions, orthogonal transformations are either reflections or rotations. The proof is based on identification of constraints imposed by reflections and rotations on the elasticity tensor, and it simply employs elementary tools from trigonometry, making the proof accessible to a broad audience. For completeness, we identify the sets of transformations (rotations and reflections) for each symmetry class and report the corresponding equations of motions in classical linear elasticity.
Thin-film coatings can be found everywhere in modern technological applications due to desirable electrical, mechanical, chemical, and optical properties. These properties directly depend upon the thin-film's microstructural features, which are themselves influenced by the materials and vapor-deposition processing conditions used for fabrication. As such, understanding processing-microstructure relationships is essential to designing thin-films with optimized properties, and discovering new processing conditions that allow for novel thin-films with multifunctional microstructures. Here, a short review is presented on recent developments that utilize the phase-field method to simultaneously model the vapor-deposition process and corresponding microstructure formation at the mesoscale. Phase-field-based vapor-deposition models that simulate thin-film growth of immiscible alloy and polycrystalline systems are highlighted in addition to machine-learning-based surrogate models that can facilitate accelerated high-fidelity simulations along with materials design and exploration studies.
Qu, Haozheng J.; Srinivasan, Jayendran; Zhao, Yangyang; Mao, Keyou S.; Taylor, Jason M.; Marino, Gabriella; Montoya, Timothy M.; Johnson, Kyle; Locke, Jenifer S.; Schaller, Rebecca S.; Schindelholz, Eric; Wharry, Janelle P.
The chloride-induced stress corrosion cracking (CISCC) mechanism of cold spray (CS) coating on a galvanically similar substrate is investigated. Arc welded 304L stainless steel (SS) specimens are loaded into four-point bend fixtures, cold sprayed with 304L SS, then immersed in boiling MgCl2. Interconnected porosity forms through crevice corrosion along CS splat boundaries, allowing corrosive species to penetrate through the CS layer. Nevertheless, the substrate is resistant to CISCC likely because of compressive stress introduced by peening during CS particle impacts. These findings underscore the importance of residual stress in the environmental degradation of CS coatings or repairs of engineering structures.
Fault detection and classification in photovoltaic (PV) systems through real-time monitoring is a fundamental task that ensures quality of operation and significantly improves the performance and reliability of operating systems. Different statistical and comparative approaches have already been proposed in the literature for fault detection; however, accurate classification of fault and loss incidents based on PV performance time series remains a key challenge. Failure diagnosis and trend-based performance loss routines were developed in this work for detecting PV underperformance and accurately identifying the different fault types and loss mechanisms. The proposed routines focus mainly on the differentiation of failures (e.g., inverter faults) from irreversible (e.g., degradation) and reversible (e.g., snow and soiling) performance loss factors based on statistical analysis. The proposed routines were benchmarked using historical inverter data obtained from a 1.8 MWp PV power plant. The results demonstrated the effectiveness of the routines for detecting failures and loss mechanisms and the capability of the pipeline for distinguishing underperformance issues using anomaly detection and change-point (CP) models. Finally, a CP model was used to extract significant changes in time series data, to detect soiling and cleaning events and to estimate both the performance loss and degradation rates of fielded PV systems.
Kinsey, Geoffrey S.; Riedel-Lyngskaer, Nicholas C.; Miguel, Alonso A.; Boyd, Matthew; Braga, Marilia; Shou, Chunhui; Cordero, Raul R.; Duck, Benjamin C.; Fell, Christopher J.; Feron, Sarah; Georghiou, George E.; Habryl, Nicholas; John, Jim J.; Ketjoy, Nipon; Lopez, Gabriel; Louwen, Atse; Maweza, Elijah L.; Mittal, Ankit; Molto, Cecile; Garrido, Gustavo N.; Norton, Matthew; Paudyal, Basant R.; Pereira, Enio B.; Poissant, Yves; Pratt, Lawrence; Shen, Qu; Reindl, Thomas; Rennhofer, Marcus; Rodriguez-Gallegos, Carlos D.; Ruther, Ricardo; Van Sark, Wilfried; Sevillano-Bendezu, Miguel A.; Seigneur, Hubert; Tejero, Jorge A.; Theristis, Marios; Tofflinger, Jan A.; Vilela, Waldeir A.; Xia, Xiangao; Yamasoe, Marcia A.
In photovoltaic power ratings, a single solar spectrum, AM1.5, is the de facto standard for record laboratory efficiencies, commercial module specifications, and performance ratios of solar power plants. More detailed energy analysis that accounts for local spectral irradiance, along with temperature and broadband irradiance, reduces forecast errors to expand the grid utility of solar energy. Here, ground-level measurements of spectral irradiance collected worldwide have been pooled to provide a sampling of geographic, seasonal, and diurnal variation. Applied to nine solar cell types, the resulting divergence in solar cell efficiencies illustrates that a single spectrum is insufficient for comparisons of cells with different spectral responses. Cells with two or more junctions tend to have efficiencies below that under the standard spectrum. Silicon exhibits the least spectral sensitivity: relative weekly site variation ranges from 1% in Lima, Peru to 14% in Edmonton, Canada.
Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. A computationally led, material design and synthesis methodology is presented here that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.
For decades, Arctic temperatures have increased twice as fast as average global temperatures. As a first step toward quantifying parametric uncertainty in Arctic climate, we performed a variance-based global sensitivity analysis (GSA) using a fully coupled, ultra-low resolution (ULR) configuration of version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SMv1). Specifically, we quantified the sensitivity of six quantities of interests (QOIs), which characterize changes in Arctic climate over a 75 year period, to uncertainties in nine model parameters spanning the sea ice, atmosphere, and ocean components of E3SMv1. Sensitivity indices for each QOI were computed with a Gaussian process emulator using 139 random realizations of the random parameters and fixed preindustrial forcing. Uncertainties in the atmospheric parameters in the Cloud Layers Unified by Binormals (CLUBB) scheme were found to have the most impact on sea ice status and the larger Arctic climate. Our results demonstrate the importance of conducting sensitivity analyses with fully coupled climate models. The ULR configuration makes such studies computationally feasible today due to its low computational cost. When advances in computational power and modeling algorithms enable the tractable use of higher-resolution models, our results will provide a baseline that can quantify the impact of model resolution on the accuracy of sensitivity indices. Moreover, the confidence intervals provided by our study, which we used to quantify the impact of the number of model evaluations on the accuracy of sensitivity estimates, have the potential to inform the computational resources needed for future sensitivity studies.
Role-based access control (RBAC) is adopted in the information and communication technology domain for authentication purposes. However, due to a very large number of entities within organizational access control (AC) systems, static RBAC management can be inefficient, costly, and can lead to cybersecurity threats. In this article, a novel hybrid RBAC model is proposed, based on the principles of offline deep reinforcement learning (RL) and Bayesian belief networks. The considered framework utilizes a fully offline RL agent, which models the behavioral history of users as a Bayesian belief-based trust indicator. Thus, the initial static RBAC policy is improved in a dynamic manner through off-policy learning while guaranteeing compliance of the internal users with the security rules of the system. By deploying our implementation within the smart grid domain and specifically within a Distributed Energy Resources (DER) ecosystem, we provide an end-To-end proof of concept of our model. Finally, detailed analysis and evaluation regarding the offline training phase of the RL agent are provided, while the online deployment of the hybrid RL-based RBAC model into the DER ecosystem highlights its key operation features and salient benefits over traditional RBAC models.
The surfaces of most metals immersed in aqueous electrolytes have a several-nanometer-thick oxide/hydroxide surface layer. This gives rise to the existence of both metal∣oxide and oxide∣liquid electrotlyte interfaces, and makes it challenging to correlate atomic length-scale structures with electrochemical properties such the potential-of-zero-charge (PZC). The PZC has been shown to be correlated the onset potential for pitting corrosion. In this work, we conduct large-scale Density Functional Theory and ab initio molecular dynamics to calculate the PZC of a Al(111)∣γ-Al2O3(110)∣ water double-interface model within the context of aluminum corrosion. By partitioning the multiple interfaces involved into binary components with additive contributions to the overall work function and voltage, we predict the PZC to be −1.53 V vs SHE for this model. We also calculate the orbital energy levels of defects like oxygen vacancies in the oxide, which are critical parameters in theories associated with pitting corrosion. We predict that the Fermi level at the PZC lies above the impurity defect levels of the oxygen vacancies, which are therefore uncharged at the PZC. From the PZC estimate, we predict the voltage needed to create oxygen vacancies with net postive charges within a flatband approximation.
An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. We show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child-Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.