Publications

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The azimuth of an incoming acoustic wave cannot be determined using microbarometers on a free floating balloon. A single observation of infrasound-induced acceleration on a large zero pressure balloon suggested that a motion sensing "aeroseismometer" could fill this gap. Here, a flight test of prototype balloon-borne aeroseismometers is presented. Two balloons, each carrying accelerometers and IMUs, recorded three sets of chemical explosions. The resulting balloon motion time series allows the explosive source to be geolocated. The future of this technology is discussed, along with a planned publication. Finally, recommendations and lessons learned from the campaign are discussed.

A phenomenological model of cavitation is presented, based on observations that both large relative negative pressures and large negative time derivatives of pressure are required for cavitation onset. We simulated two cavitation experiments to generate cavitation scaling parameters for relative pressure drop and rate of pressure drop. Our results show the model, while simple, is effective at reproducing results from laboratory experiments of cavitation. The parameters were then used in conjunction with a human surrogate computational model to predict, at any position within the head, the probability of intracranial cavitation caused by exposure to a blast event. The results suggest that the magnitude of blast overpressure observed in field data is sufficient to cause intracranial cavitation. Our analysis indicates that the helmeted head, when compared to the unhelmeted head configuration, results in a decrease but not elimination of cavitation exposure. When density functions of cavitation probability versus cumulative brain volume are combined with an injury severity model, the results show helmet efficacy at low and moderate risk levels. However, the convergence of unhelmeted and helmeted probability density functions at high-to-excessive risk thresholds indicates the helmet offers diminishing protection at elevated exposure levels, relative to the unhelmeted baseline. Future investigation and collaboration with neuroscience subject matter experts are needed to contextualize the current computational results. While the present work contributes specific and quantified predictions of intracranial cavitation location and severity, more research is required to apply our results to clinical settings with population-based brain injury subjects and controls. The relationship between our intracranial cavitation predictions with their anticipated clinical sequelae remains a topic in need of exploration.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In October of 2020, dust samples were collected from the surface of spent nuclear fuel (SNF) dry storage canisters during an inspection at an inland Independent Spent Fuel Storage Installation, the second inland site at which surface deposits have been sampled. The purpose of the sampling was to assess the composition and abundance of the soluble salts present on the canister surface, information which provides a metric for potential corrosion risks. 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 texture and mineralogy of the dust and salts. The results of those analyses are presented in this report.

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.

Abstract not provided.

This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.

Abstract not provided.

Abstract not provided.

The Statement of Work, Deliverables and Costs for the project are provided.

Abstract not provided.

A stripmap Synthetic Aperture Radar (SAR) image is a long SAR image along some centerline, and formed from multiple synthetic apertures. At issue is that the centerline in the image actually corresponds to an arc on a round earth, and multiple strategies exist for fitting the image centerline to the round earth. Some of those strategies involve Rhumb lines, great circle paths, and great ellipse paths. Some are better than others in polar regions. Notions of parallel flight paths for the radar during data collection also require careful consideration of the geometry of a round earth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 14959and 18019. Revision 14959 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 18019. Along with the newly updated MELCOR Users Guide and Reference Manual, users are aware and able to assess the new capabilities for their modeling and analysis applications.

The motivation for this report is to discuss and present some realistic tree models employed in computational electromagnetics (EM) simulations to study foliage penetration (FOPEN) at Ku-band. The detail obtained in these trees is unprecedented in FOPEN modeling since many studies in this area focus on lower frequencies where precise tree parameters are not required due to the associated large wavelengths relative to the tree dimensions. The focus of this study is in the Ku-band range where the wavelength is notably smaller and the details of the trees have more of an influence on EM waves (i.e. scattering, attenuating, reflecting, diffracting etc.). Therefore, explicit tree parameters are modeled. Also, moderate foliage is of most interest because with less dense foliage t here is a higher percentage of Ku-band transmission. The EM wave and foliage interaction s are simulated with the computational electromagnetics (CEM) Altair FEKO software. The realistic tree model s implemented for simulations are created in the computer-aided design (CAD) software Arbaro and the module CADFEKO that is offered in FEKO. Details of these tree models are provided, and EM simulation results will be discussed in a follow-on report

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report describes the design, development, and testing of a prototype 100 kWt particle-to-supercritical CO₂ (sCO2) heat exchanger. An analytic hierarchy process was implemented to compare and evaluate alternative heat-exchanger designs (fluidized bed, shell-and-plate moving packed bed, and shell-and-tube moving packed bed) that could meet the high pressure (≥ 20 MPa) and high temperature (≥ 700 °C) operational requirements associated with sCO2 power cycles. Cost, heat-transfer coefficient, structural reliability, manufacturability, parasitics and heat losses, scalability, compatibility, erosion and corrosion, transient operation, and inspection ease were considered in the evaluation. A 100 kWt shell-and-plate design was selected for construction and integration with Sandia’s falling particle receiver system that heats the particles using concentrated sunlight. Sandia worked with industry to design and construct the moving packed-bed shell-and-plate heat exchanger. Tests were performed to evaluate its performance using both electrical heating and concentrated sunlight to heat the particles. Overall heat transfer coefficients at off-design conditions (reduced operating temperatures and only three stainless steel banks in the counter-crossflow heat exchanger) were measured to be approximately ~25 - 70 W/m²-K, significantly lower than simulated values of >100 W/m²-K. Tests using the falling particle receiver to heat the particles with concentrated sunlight yielded overall heat transfer coefficients of ~35 – 80 W/m²-K with four banks (including a nickel-alloy bank above the three stainless steel banks). The overall heat transfer coefficient was observed to decrease with increasing particle inlet temperatures, which contrasted the results of simulations that showed an increase in heat transfer coefficient with temperature due to increased effective particle-bed thermal conductivity from radiation. The likely cause of the discrepancy was particle-flow maldistributions and funnel flow within the heat exchanger caused by internal ledges and cross-bracing, which could have been exacerbated by increased particle-wall friction at higher temperatures. Additional heat loss at higher temperatures may also contribute to a lower overall heat-transfer coefficient. Design challenges including pressure drop, particle and sCO2 flow maldistribution, and reduced heat transfer coefficient are discussed with approaches for mitigation in future designs. Lessons learned regarding instrumentation, performance characterization, and operation of particle components and sCO2 flow loops are also discussed. Finally, a 200 MWt commercial-scale shell-and-plate heat-exchanger design based on the concepts investigated in this report is proposed.

Abstract not provided.

Abstract not provided.

This paper presents a multifidelity uncertainty quantification framework called MFNets. We seek to address three existing challenges that arise when experimental and simulation data from different sources are used to enhance statistical estimation and prediction with quantified uncertainty. Specifically, we demonstrate that MFNets can (1) fuse heterogeneous data sources arising from simulations with different parameterizations, e.g simulation models with different uncertain parameters or data sets collected under different environmental conditions; (2) encode known relationships among data sources to reduce data requirements; and (3) improve the robustness of existing multi-fidelity approaches to corrupted data. MFNets construct a network of latent variables (LVs) to facilitate the fusion of data from an ensemble of sources of varying credibility and cost. These LVs are posited as explanatory variables that provide the source of correlation in the observed data. Furthermore, MFNets provide a way to encode prior physical knowledge to enable efficient estimation of statistics and/or construction of surrogates via conditional independence relations on the LVs. We highlight the utility of our framework with a number of theoretical results which assess the quality of the posterior mean as a frequentist estimator and compare it to standard sampling approaches that use single fidelity, multilevel, and control variate Monte Carlo estimators. We also use the proposed framework to derive the Monte Carlo-based control variate estimator entirely from the use of Bayes rule and linear-Gaussian models -- to our knowledge the first such derivation. Finally, we demonstrate the ability to work with different uncertain parameters across different models.

Sandia National Laboratories’ (Sandia’s) Electrical Sciences Group (1350) provides a spectrum of solutions for electromagnetic (EM) environment effects on electrical systems to advance physical understanding including experiments, analytical and numerical model development leveraging first-principles physics, and high-performance computational modeling and simulation (M&S). The six departments in the Electrical Sciences Group provide analysis, design guidance, and experiments in support of nuclear weapons qualification in normal, abnormal, and hostile environments, as well as research and development for advanced electrical systems that can operate through these environments. This document is intended to provide a quick summary of Sandia’s EM experimental facilities, most of which provide unique national capabilities.

By 2030 about half of all spent nuclear fuel (SNF) arising from the current fleet of commercial power plants will be in dual-purpose canisters (DPCs), which are designed for storage and transportation but not for disposal. As an alternative to complete repackaging of the fuel for disposal, considerable cost savings and lower worker dose could be realized by directly disposing of this SNF in DPCs. The principal technical consideration is criticality control in a geologic repository, because the DPCs are large and depend on neutron absorbing basket components for criticality control. Neutron absorbing materials are generally aluminum-based, and under disposal conditions can degrade after a few hundred years contact with ground water. Simple modifications to the SNF assemblies or the DPC baskets could help to achieve direct disposal, and this is one of the approaches being studied to address the possibility of disposal criticality (SNL 2020a). Five fuel/basket modification concepts have been proposed (SNL 2020b) and a virtual workshop was conducted to solicit review and feedback on these concepts. The proposed solutions are: 1) zone loading of DPCs to limit reactivity, 2) replacing absorber plates with advanced neutron absorbing (ANA) material, 3) adding disposal control rods to pressurized water reactor (PWR) assemblies, 4) rechanneling boiling water reactor (BWR) assemblies with ANA material, and 5) basket insert plates (chevron inserts) made from ANA material. The presentations from the workshop are provided in this report, and the workshop discussions are summarized. This information includes prioritization of the proposed fuel/basket modification solutions, and prioritization of the associated model development, validation testing, and quality assurance activities. Information documented in this report will help to steer research and development efforts at Sandia National Laboratories, Oak Ridge National Laboratory, and Idaho National Laboratory that support the U.S. Department of Energy, Office of Nuclear Energy, Spent Fuel and Waste Science and Technology program

Tools are now available that enable measurement electromagnetic radiation (EMR) from active electronics in an item. This radiation may be intended WIFI or cellular network links, for example or unintended such as the switching noise generated by DC-to-DC converters. It would be extremely valuable to have the capability to discriminate between the low-voltage DC-to-DC converters or other digital noise prevalent in most modern electronics, versus the high-voltage DC-to-DC converters used in utility firesets. Previous work performed under a Sandia Laboratory Directed Research and Development (LDRD) project on Charge State Detection using a deep neural network has been continued in this effort. A state-of-the-art supervised machine learning algorithm has not only been extended to discriminate between low and high voltage converters but has been validated in determining a converters make and model.

The vertical borehole array at Farnsworth Unit, TX is used to monitor microseismic activity in the subsurface around the Carbon Capture and Sequestration (CCS) reservoir. The array consists of 16 3-component seismometers spaced vertically in a single borehole. Tube or borehole waves traveling up or down the borehole can corrupt signals of interest, such as microseismic events. A denoising convolutional neural network (DCNN) was trained to remove borehole waves from seismic waveforms of microseismic events for the purpose of reducing unwanted signal detections and better characterizing events of interest. This R&D leverages the work of Sandia colleague Rigo Tibi, who used a DCNN developed by Greg Beroza at Stanford University to improve the signal-to-noise ratio (SNR) of teleseismic events detected by the International Monitoring System.

This report summarizes the key contributions and lessons learned from SNL experience in technical reviews of Controls awardees in the DOE SPA program from 2013 - 2020. The purpose of this report is to provide observations and technical suggestions that are likely to be beneficial to the WEC industry as a whole. Over the course of the SPA FOA program, SNL has engaged in technical review for a total of 5 different Controls awardees. The awardees represent a diversity of WEC devices and the application of different control design approaches. The report begins with a summary of key performance metrics results reported by the 5 Controls awardees. This is followed by a summary of observations and lessons learned distilled from the technical reviews of the awardees . The report concludes with a list of general technical suggestions for future WEC controls projects.

Natural and induced fractures are potential preferential pathways for migration of radioactive gases to earths surface from underground nuclear explosions (UNEs). This report documents X-ray computed tomography (XRCT) imaging on 26 samples of rock core that was collected to support the Underground Nuclear Explosion Signatures Experiment (UNESE) program. The XRCT datasets are intended to help fill a data gap on the three-dimensional (3D) characteristics of natural and/or induced fractures at the centimeter and smaller scale, which may strongly influence multiphase fluid flow and transport properties of preferential flow paths and interaction with the matrix of the surrounding host rock. Pre- and post-UNE rock samples were carefully chosen to enable comparison of fractures as a function of lithologic and petrophysical properties, as well as distance to the past UNEs. This report serves as documentation for the data, including an introduction with the research motivation, a methods and materials section, descriptions of the XRCT datasets without post-processing, and recommendations for 3D quantification via image analysis and digital rock physics.

Vaccination and the alternative behavior, vaccine refusal, are a classic example of manifesting behaviors driven by social norms and norm violations. Establishing how norms emerge, and under what circumstances people choose to violate them are key issues to understand in modeling epidemics. Interactions between individuals can lead to large-scale patterning of behavior (emergent phenomena). As norm violations are revealed through human behavior, drawing on psychological theory and principles to predict those violations is a viable approach for more human-constrained epidemiological models. As an example of the implications at scale, vaccine refusal is correlated with the spread of mis/disinformation about vaccine side-effects. Considering the complexities of network dynamics, the downstream effects means that if even a small group within a population are persuaded against vaccination, there is a reservoir from which disease and disease outbreaks can propagate. This work will attempt to identify those psychological indicators, to define circumstances that predict health behaviors, and identify potentially modifiable antecedents of health behavior, and factors that influence changes toward health protective behaviors.

This report summarizes design optimization and performance evaluation studies for a new prototype CONFIDANTE (CONfirmation using a Fast-neutron Imaging Detector with Anti-image Null-positive Time Encoding) warhead confirmation system. It was found that a spherical mask geometry and a 2” diameter cylindrical central detector is expected to best discriminate between the parametrically varied source distributions that were evaluated. The optimized design as fabricated and its performance was evaluated in a series of laboratory measurements. The performance was in good agreement with the design studies, with demonstrated discrimination between objects with differences in scale on the order of 5 centimeters or better at 1 meter stand-off

Dispensers are the top cause of maintenance events and down-time at hydrogen fueling stations. In an effort to help characterize and enable improvements in dispenser reliability, an extensive accelerated lifetime testing set-up was designed and built at NREL involving components typically part of dispensing operations at fueling stations. Device Under Test (DUTs) included different components such as normally open valves, normally closed valves, fueling nozzles, breakaways devices and filters. Conditions of testing included pressures, and flow rates similar to light duty fuel cell electric vehicles fueling at -40°C, and -20°C for thousands of cycles in hydrogen. Tested components (failed and non-failed) were disassembled at SNL and polymeric O-rings were carefully retrieved and cataloged for chemical and physical characterization. Data collected was compared to similar O-rings from unexposed or non-tested components for hydrogen effects, and failure modes. Degradation analyses, based on select polymer chemistries common across all component types, their location within components, visual assessment of damage coupled with strong hydrogen effects from chemical characterization, was completed and presented to NREL and DOE. Overall, the failure rate amongst the components was not as high as expected for the test conditions. Among the component types tested, breakaways were the most susceptible to damage under these test conditions, with fueling nozzles a close second. The proper combination of selection of the right polymer and optimum component design was found to make a strong difference in component reliability under severe dispenser operating conditions. Physical degradation of polymers, rather than chemical changes due to low temperature hydrogen exposure, is more prevalent as failure mode for these test conditions. The nature and the extent of the degradation was much less at -20°C as compared to -40°C. The damage and failure rates were higher at lower temperatures than at higher test temperatures. As expected, increasing the number of cycles at the lowest test temperature (-40°C) increased damage. This indicates that cycling at the low temperature of -40°C required by SAE J2601 can reduce component life in fuel dispensing operations

Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Microscopically, how does the hydration/dehydration process occur?

Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.

The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO₂ has been captured in a tamped Richtmyer–Meshkov instability (RMI) environment at a pressure regime of 4–10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker–Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.

Abstract not provided.

Grid-level energy storage systems are needed to enable intermittent renewables. Li-ion, Pb-acid battery systems have been implemented but pose safety and environmental risks. Successful grid storage must be safe, reliable, and low-cost.

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

Abstract not provided.

Grid-level energy storage systems are needed to enable intermittent renewables. Li-ion, Pb-acid battery systems have been implemented but pose safety and environmental risks. Successful grid storage must be safe, reliable, and low-cost.

Abstract not provided.

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

Cryogenic environments make superconducting computing possible by reducing thermal noise, electrical resistance and heat dissipation. Heat generated by the electronics and thermal conductivity of electrical transmission lines to the outside world constitute two main sources of thermal load in such systems. As a result, higher data rates require additional transmission lines which come at an increasingly higher cooling power cost. Hybrid or monolithic integration of silicon photonics with the electronics can be the key to higher data rates and lower power costs in these systems. We present a 4-channel wavelength division multiplexing photonic integrated circuit (PIC) built from modulators in the AIM Photonics process development kit (PDK) that operate at 25 Gbps at room temperature and 10 Gbps at 40 K. We further demonstrate 2-channel operation for 20 Gbps aggregate data rate at 40 K using two different modulators/wavelengths, with the potential for higher aggregate bit rates by utilizing additional channels.

This report provides an analysis of the clad barrier function associated with the direct disposal of dual purpose canisters (DPCs) under hypothetical conditions in a shale repository and in an alluvial repository, including the effect of a postulated criticality event inside a disposed DPC. Should a postulated criticality event occur in a hypothetical shale repository, cladding will primarily degrade by general corrosion. Stress corrosion cracking, hydride cracking, creep failure, pitting and crevice corrosion, rod pressurization, and clad unzipping are calculated to have little impact on cladding persistence. At the higher temperature expected during a postulated criticality event in a saturated shale repository, general corrosion of cladding would be rapid - on the order of 0.034 microns/yr. A few hundred years after onset of a postulated criticality event in a shale repository complete general corrosion of fuel assembly grid spacer walls and guide tubes will likely result in settling of fuel rods upon each other. This rod consolidation should displace the water moderator and possibly terminate a postulated criticality. The primary potential degradation pathway for cladding in a hypothetical alluvial repository is localized corrosion by fluoride, which cannot occur in a shale repository. Fluoride-enhanced corrosion of cladding would be accelerated under the slightly higher (< 100°C) temperatures associated with a postulated criticality event. The impact of criticality in both cases (shale and alluvial) would be to increase the amount of failed cladding. But it would require very specialized transport pathways.

This report provides a high-level test plan for deploying three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUHOMS Advanced Horizontal Storage Modules (AHSM) from Orano (formerly Transnuclear Inc.). The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or adapt in response to conflicting requirements. The goal of the testing is to collect highly defensible and detailed surface deposition measurements from the surface of dry storage systems in a marine coastal environment to guide chloride-induced stress corrosion crack (CISCC) research. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the thermal-hydraulic environment. Instrumentation throughout the canister, storage module, and environment will provide an extensive amount of information for the use of model validation. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer a high-fidelity quantification of the conditions experienced in a harsh yet realistic environment.

Sandia National Laboratories Radiation Protection department did not successfully address the causes leading to the issue NTS-SS-SNL-NMSITE-2017-0004, “Contamination Monitor Left in Service After Inadequate Calibration”, nor identify through the verification and validation process that corrective actions were not fully implemented and effective at addressing the issue(s). This report is the third causal in four years associated with the iPCM12 instrument and the Radiation Protection Instrumentation (RPI) organization within Radiation Protection (00628).

In this paper, the magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics simulations coupled to molecular dynamics, using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M (T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M (T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering. Analytical approximations to the magnetization as a function of temperature and size are proposed.

This study reveals the unique role on Janus particles of the solid-solid interface at the boundary in determining particle interactions and assembly. In an aqueous ionic liquid (IL) solution, Janus spheres adopt intriguing orientations with their boundaries pinned on the glass substrate. It was further discovered that the orientation was affected by the particle amphiphilicity as well as the chemical structure and concentration of the IL. Further characterization suggests that the adsorption on the hydrophilic side is due to both an electrostatic interaction and hydrogen bonding, while adsorption on the hydrophobic side is due to hydrophobic attraction. Through the concerted interplay of all these interactions, the amphiphilic boundary may attract an excessive amount of IL cations, which guide the unique orientations of the Janus spheres. The results highlight the importance of the Janus boundary that has not been recognized previously. Adsorption at the solid-solid interfaces may inspire new applications in areas such as separation and catalysis.

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

On June 5th, 2020 Sandia National Laboratories signed and recommended for closure the evidence package for the Noncompliance Tracking System NTS—SS-SNL-NMSITE-2017-0004, Contamination Monitor Left in Service After Inadequate Calibration. On July 16th, 2020 the Sandia Field Office sent a follow-up question list to the Radiation Protection Instrumentation Program Lead to clarify some points while reviewing the evidence package. The follow-up questions readily identified deficiencies with all three calibrations of the iPCM12 systems. On July 30th, 2020 Radiation Protection hosted a meeting to better understand the exact nature of the concerns, determine if the iPCM12 machines were out of calibration and therefore would have to be taken out of service. One result of this meeting was to determine if SNL could withdraw the evidence package. SNL transmitted a letter on August 26th, 2020 recalling the evidence package and committed to performing a thorough causal analysis and invoke the required issues management processes to ensure that corrective actions are sustained.

A new scatter calculation algorithm has been implemented in the Gamma Detector Response and Analysis Software (GADRAS) package that accounts for spectral effects of scattering materials not in the line of sight of the detector and the source. Previously, GADRAS would only apply scattering effects due to materials that fall between the source and detector. This new routine will allow better modeling of various scenarios including gamma imagers, collimated detectors, or traditional gamma detectors where scattering materials are present.

Producing and distributing COVID-19 vaccine during the pandemic is a major logistical challenge requiring careful planning and efficient execution. This report presents information on logistical, policy and technical issues relevant to rapidly fielding a COVID-19 vaccination program. For this study we (a) conducted literature review and subject matter expert elicitation to understand current vaccine manufacturing and distribution capabilities and vaccine allocation strategies, (b) designed a baseline vaccine distribution strategy and modeling strategy to provide insight into the potential for targeted distribution of limited initial vaccine supplies, and (c) developed parametric interfaces to enable vaccine distribution scenarios to be analyzed in depth with Sandias Adaptive Recovery Model that will allow us evaluate the additional sub- populations and alternative distribution scenarios from a public health benefit and associated economic disruption Principal issues, challenges, and complexities that complicate COVID-19 vaccine delivery identified in our literature and subject matter expert investigation include these items: The United States has not mounted an urgent nationwide vaccination campaign in recent history. The existing global manufacturing and distribution infrastructure are not able to produce enough vaccine for the population immediately. Vaccines, once available will be scarce resources. Prioritization for vaccine allocation will be built on existing distribution networks. Vaccine distribution may not have a universal impact on disease transmission and morbidity because of scarcity, priority population demographics, and underlying disease transmission rates. Considerations for designing a vaccine distribution strategy are discussed. A baseline distribution strategy is designed and tested using the Adaptive Recovery Model, which couples a deterministic compartmental epidemiological model and a stochastic network model. We show the impact of this vaccine distribution strategy on hospitalizations, mortality, and contact tracing requirements. This model can be used to quantitatively evaluate alternative distribution scenarios, guiding policy decisions as vaccine candidates are narrowed down.

I am a staff scientist at Sandia National Laboratories (SNL). I work on multiple fundamental-science projects and lead modeling/data-analysis for the stellar opacity experiments. I also work with theorists to refine plasma material-property calculations to make NNSA simulations more predictable. These are challenging and important problems for national security. In fact, our experiments raised questions about “opacity”, the property of matter that controls energy transport inside stars. This work was published by Nature due to its serious implication over broad applications, and our team received NNSA Defense Program Award of Excellence, both in 2015. Work at SNL is collegial and full of learning through interactions with the world’s finest experimentalists and theorists. When I started my Ph.D. program, working at a national lab was a dream position. Today, I collaborate with top scientists to solve mission-critical problems. The Stewardship Science Academic Programs (SSAP) and National Laser Users’ Facility (NLUF) played a vital role in the efficient learning and seamless transition to SNL.

The use of multiple types of precision in mathematical software has the potential to increase its performance on new heterogeneous architectures. The xSDK project focuses both on the investigation and development of multiprecision algorithms as well as their inclusion into xSDK member libraries. This report summarizes current efforts on including and/or using mixed precision capabilities in the math libraries Ginkgo, heFFTe, hypre, MAGMA, PETSc/TAO, SLATE, SuperLU, and Trilinos, including KokkosKernels. It contains both numerical results from libraries that already provide mixed precision capabilities, as well as descriptions of the strategies to incorporate multiprecision into established libraries.

Abstract not provided.

The equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50000 K, for densities between -ü=3 and 15g/cm3, with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

The thirty-year non-stationary historical trends in the wave energy climate for United States coastal waters between 1980 and 2009 are investigated using spectral partitioned wave data generated from a WaveWatch III® (version 5.05) hindcast. In addition to historical trends in the omni-directional wave power, frequency and directionally resolved wave power, frequency and directional spreading, and seasonal variability, are examined for the first time, including their geographical distribution. These historical wave energy climate trends are linked to changes to the dominant wave systems and commensurate trends in the historical wind climate. Total wave power trends are consistent with other studies, but the present study identifies regions, and specific frequency and direction bands, where significant wave energy climate changes have. For some regions, changes to wave climate parameters exceeded one-percent annually, more than thirty-percent over the study period. Non-stationary trends of this magnitude have significant implications for ocean and coastal engineering projects designed assuming wave climates are stationary and warrant their consideration in design practices.

Neuromorphic computing is a critical future technology for the computing industry, but it has yet to achieve its promise and has struggled to establish a cohesive research community. A large part of the challenge is that full realization of the potential of brain inspiration requires advances in both device hardware, computing architectures, and algorithms. This simultaneous development across technology scales is unprecedented in the computing field. This article presents a strategy, framed by market and policy pressures, for moving past these current technological and cultural hurdles to realize its full impact across technology. Achieving the full potential of brain-derived algorithms as well as post-complementary metal-oxide-semiconductor (CMOS) scaling neuromorphic hardware requires appropriately balancing the near-term opportunities of deep learning applications with the long-term potential of less understood opportunities in neural computing.

This work presents a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.

High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO), hydroperoxyalkyl peroxy (OOQOOH), and ketohydroperoxide (HOOPO, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. With the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO, OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO, OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O2 concentration, and temperature on the competition among radical decomposition and second O2 addition pathways.

Steady-state photocapacitance (SSPC) was conducted on nonpolar m-plane GaN n-type Schottky diodes to evaluate the defects induced by inductively coupled plasma (ICP) dry etching in etched-and-regrown unipolar structures. An ∼10× increase in the near-midgap Ec - 1.9 eV level compared to an as-grown material was observed. Defect levels associated with regrowth without an etch were also investigated. The defects in the regrown structure (without an etch) are highly spatially localized to the regrowth interface. Subsequently, by depth profiling an etched-and-regrown sample, we show that the intensities of the defect-related SSPC features associated with dry etching depend strongly on the depth away from the regrowth interface, which is also reported previously [Nedy et al., Semicond. Sci. Technol. 30, 085019 (2015); Fang et al., Jpn. J. Appl. Phys. 42, 4207-4212 (2003); and Cao et al., IEEE Trans. Electron Devices 47, 1320-1324 (2000)]. A photoelectrochemical etching (PEC) method and a wet AZ400K treatment are also introduced to reduce the etch-induced deep levels. A significant reduction in the density of deep levels is observed in the sample that was treated with PEC etching after dry etching and prior to regrowth. An ∼2× reduction in the density of Ec - 1.9 eV level compared to a reference etched-and-regrown structure was observed upon the application of PEC etching treatment prior to the regrowth. The PEC etching method is promising for reducing defects in selective-area doping for vertical power switching structures with complex geometries [Meyers et al., J. Electron. Mater. 49, 3481-3489 (2020)].

Machine learning (ML) techniques are being used to detect increasing amounts of malware and variants. Despite successful applications of ML, we hypothesize that the full potential of ML is not realized in malware analysis (MA) due to a semantic gap between the ML and MA communities-as demonstrated in the data that is used. Due in part to the available data, ML has primarily focused on detection whereas MA is also interested in identifying behaviors. We review existing open-source malware datasets used in ML and find a lack of behavioral information that could facilitate stronger impact by ML in MA. As a first step in bridging this gap, we label existing data with behavioral information using open-source MA reports-1) altering the analysis from identifying malware to identifying behaviors, 2)~aligning ML better with MA, and 3)~allowing ML models to generalize to novel malware in a zero/few-shot learning manner. We classify the behavior of a malware family not seen during training using transfer learning from a state-of-the-art model for malware family classification and achieve 57%-84% accuracy on behavioral identification but fail to outperform the baseline set by a majority class predictor. This highlights opportunities for improvement on this task related to the data representation, the need for malware specific ML techniques, and a larger training set of malware samples labeled with behaviors.

To combat dynamic, cyber-physical disturbances in the electric grid, online and adaptive remedial action schemes (RASs) are needed to achieve fast and effective response. However, a major challenge lies in reducing the computational burden of analyses needed to inform selection of appropriate controls. This paper proposes the use of a role and interaction discovery (RID) algorithm that leverages control sensitivities to gain insight into the controller roles and support groups. Using these results, a procedure is developed to reduce the control search space to reduce computation time while achieving effective control response. A case study is presented that considers corrective line switching to mitigate geomagnetically induced current (GIC) -saturated reactive power losses in a 20-bus test system. Results demonstrated both significant reduction of both the control search space and reactive power losses using the RID approach.

The swelling of clay at high temperature and pressure is important for applications including nuclear waste storage but is not well understood. A molecular dynamics study of the swelling of Na montmorillonite in water at several temperatures (T = 298, 400, and 500 K) and water environment pressures (P_e = 5 and 100 MPa) is reported here. Adopting a rarely used setup that enables swelling pressure to be resolved with an accuracy of ~1 MPa, the swelling pressure was computed at basal spacings 1.6–2.6 nm (or 2–5 water layers between neighboring clay sheets), which has not been widely studied before. At T = 298 K and P_e = 5 MPa, swelling pressure P_s oscillates at d-spacing d smaller than 2.2 nm and decays monotonically as d increases. Increasing T to 500 K but keeping P_e at 5 MPa, P_s remains oscillatory at small d, but its repulsive peak at d = 2.2 nm shifts to ~2.0 nm and P_s at different d-spacings can grow more attractive or repulsive. At d > 2.0 nm, P_s is weakened greatly. Keeping T at 500 K and increasing P_e to 100 MPa, P_s recovers toward that at T = 298 K and Pe = 5 MPa, however, the repulsive peak at d = 2.0 nm remains the same. The opposite effects of increasing temperature and pressure on the density and dielectric screening of water, which control ion correlations and thus double layer repulsion, are essential for understanding the observed swelling pressure at elevated temperatures and its response to environment pressures.

Alkaline zinc-manganese dioxide (Zn-MnO2) batteries are well suited for grid storage applications because of their inherently safe, aqueous electrolyte and established materials supply chain, resulting in low production costs. With recent advances in the development of Cu/Bi-stabilized birnessite cathodes capable of the full 2-electron capacity equivalent of MnO2 (617 mA h/g), there is a need for selective separators that prevent zincate (Zn(OH)4)2- transport from the anode to the cathode during cycling, as this electrode system fails in the presence of dissolved zinc. Herein, we present the synthesis of N-butylimidazolium-functionalized polysulfone (NBI-PSU)-based separators and evaluate their ability to selectively transport hydroxide over zincate. We then examine their impact on the cycling of high depth of discharge Zn/(Cu/Bi-MnO2) batteries when inserted in between the cathode and anode. Initially, we establish our membranes' selectivity by performing zincate and hydroxide diffusion tests, showing a marked improvement in zincate-blocking (DZn (cm2/min): 0.17 ± 0.04 × 10-6 for 50-PSU, our most selective separator vs 2.0 ± 0.8 × 10-6 for Cellophane 350P00 and 5.7 ± 0.8 × 10-6 for Celgard 3501), while maintaining similar crossover rates for hydroxide (DOH (cm2/min): 9.4 ± 0.1 × 10-6 for 50-PSU vs 17 ± 0.5 × 10-6 for Cellophane 350P00 and 6.7 ± 0.6 × 10-6 for Celgard 3501). We then implement our membranes into cells and observe an improvement in cycle life over control cells containing only the commercial separators (cell lifetime extended from 21 to 79 cycles).

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. Here we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5‐km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100‐kHz, continuous seismic monitoring system deployed in six 60‐m length monitoring boreholes surrounding the experimental domain in 3‐D. The achieved location uncertainty was on the order of 1 m and limited by the signal‐to‐noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear‐reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well‐behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

Passive silicon photonic waveguides are exposed to gamma radiation to understand how the performance of silicon photonic integrated circuits is affected in harsh environments such as space or high energy physics experiments. The propagation loss and group index of the mode guided by these waveguides is characterized by implementing a phase sensitive swept-wavelength interferometric method. We find that the propagation loss associated with each waveguide geometry explored in this study slightly increases at absorbed doses of up to 100 krad (Si). The measured change in group index associated with the same waveguide geometries is negligibly changed after exposure. Additionally, we show that the post-exposure degradation of these waveguides can be improved through heat treatment.

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

Impact problems of plate-like parts by punch-like objects with relatively large mass moving at slow speeds of a few feet per second constitute a subset of impact problems of interest at Sandia. This is in contrast to small objects moving in the range of hundreds or thousands of feet per second or higher. The objective of this work is to develop a simple formula that can be used to estimate a lower bound for the puncture energy of metal plates impacted by cylindrical, essentially rigid punches of circular cross-section and at nose. Such geometry is used as a basis in the design of puncture mitigation barriers or procedures. This was accomplished by deriving an expression using non-dimensional analysis and then calibrating it based on tests results in the range of speeds of interest. Lower bounds can then be determined based on confidence intervals or factors of safety.

Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.

Researchers at Sandia National Laboratories have integrated the GRANTA materials database with the MatCal calibration engine to calibrate material models from characterization data. GRANTA is gaining acceptance across the NNSA Tri-lab complex and is being populated with weapons-specific test data by Sandia experimentalists. To use that data to create material models for use by weapons systems analysts, MatCal has been enabled import calibration data and test conditions from GRANTA to quickly and reproducibly produce a calibrated set of parameters for a given constitutive model. The team is currently working to store the parameters characterizing material behavior in GRANTA to make them accessible by all weapons analysts.

As the US electrifies the transportation sector, cyber attacks targeting vehicle charging could bring consequences to electrical system infrastructure. This is a growing area of concern as charging stations increase power delivery and must communicate to a range of entities to authorize charging, sequence the charging process, and manage load (grid operators, vehicles, OEM vendors, charging network operators, etc.). The research challenges are numerous and are complicated because there are many end users, stakeholders, and software and equipment vendors interests involved. Poorly implemented electric vehicle supply equipment (EVSE), electric vehicle (EV), or grid communication system cybersecurity could be a significant risk to EV adoption because the political, social, and financial impact of cyberattacks - or public perception of such - ripples across the industry and has lasting and devastating effects. Unfortunately, there is no comprehensive EVSE cybersecurity approach and limited best practices have been adopted by the EV/EVSE industry. There is an incomplete industry understanding of the attack surface, interconnected assets, and unsecured interfaces. Thus, comprehensive cybersecurity recommendations founded on sound research are necessary to secure EV charging infrastructure. This project is providing the power, security, and automotive industry with a strong technical basis for securing this infrastructure by developing threat models, determining technology gaps, and identifying or developing effective countermeasures. Specifically, the team is creating a cybersecurity threat model and performing a technical risk assessment of EVSE assets, so that automotive, charging, and utility stakeholders can better protect customers, vehicles, and power systems in the face of new cyber threats.

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of inertial confinement fusion (ICF) experiments. The Grüneisen parameter increases with density below ∼3.1g/cm3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. The discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

The incorporation of nanostructured and amorphous metals into modern applications is reliant on the understanding of deformation and failure modes in constrained conditions. To study this, a 105 nm crystalline Cu/160 nm amorphous Cu45Zr55 (at.%) multilayer structure was fabricated with the two crystalline layers sputter deposited between the top-middle-bottom amorphous layers and prepared to electron transparency. The multilayer was then in situ indented either under a single load to a depth of ~ 100 nm (max load of ~ 100 μN) or held at 20 μN and then repeatedly indented with an additional 5 μN up to 20,000 cycles in a transmission electron microscope to compare the deformation responses in the nanolaminate. For the single indentation test, the multilayer showed serrated load-displacement behavior upon initial indentation inductive of shear banding. At an indentation depth of ~ 32 nm, the multilayer exhibited perfect plastic behavior and no strain hardening. Both indented and fatigue-indented films revealed diffraction contrast changes with deformation. Subsequent Automated Crystal Orientation Mapping (ACOM) measurements confirmed and quantified global texture changes in the crystalline layers with specifically identified grains revealing rotation. Using a finite element model, the in-plane displacement vectors under the indent mapped conditions where ACOM determined grain rotation was observed, indicating the stress flow induced grain rotation. The single indented Cu layers also exhibited evidence of deformation induced grain growth, which was not evident in the fatigue-indented Cu based multilayer. Finally, the single indented multilayer retained a significant plastic crater in the upper most amorphous layer that directly contacted the indenter; a negligible crater impression in the same region was observed in the fatigued tested multilayer. These differences are explained by the different loading methods, applied load, and deformation mechanisms experienced in the multilayers.

The reactivity of liquid hydrazine (N2H4) with respect to H-, Cl-, and Br-terminated Si(100) surfaces was investigated to uncover the principles of nitrogen incorporation into the interface. This process has important implications in a wide variety of applications, including semiconductor surface passivation and functionalization, nitride growth, and many others. The use of hydrazine as a precursor allows for reactions that exclude carbon and oxygen, the primary sources of contamination in processing. In this work, the reactivity of N2H4 with H- and Cl-terminated surfaces prepared by traditional solvent-based methods and with a Br-terminated Si(100) prepared in ultrahigh vacuum was compared. The reactions were studied with X-ray photoelectron spectroscopy, atomic force microscopy, and scanning tunneling microscopy, and the observations were supported by computational investigations. The H-terminated surface led to the highest level of nitrogen incorporation; however, the process proceeds with increasing surface roughness, suggesting possible etching or replacement reactions. In the case of Cl-terminated (predominantly dichloride) and Br-terminated (monobromide) surfaces, the amount of nitrogen incorporation on both surfaces after the reaction with hydrazine was very similar despite the differences in preparation, initial structure, and chemical composition. Density functional theory was used to propose the possible surface structures and to analyze surface reactivity.

Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.

The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe_0.70Ni_0.11Cr_0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.

Here, we describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.

Shock-induced detonation is a key property of energetic materials (EM) that remains empirically understood. One proposed mechanism of shock-initiation in EM is “phonon up-pumping” to initiate chemical reactions, where excitation of lattice phonon modes rapidly transfers energy into intramolecular vibrations, ultimately resulting in the breaking of chemical bonds. We are developing novel ultrafast laser spectroscopy techniques to study vibrational energy transfer from phonon modes to intramolecular vibrations (phonon up-pumping), as well as competing energy transfer pathways from intramolecular vibrations to phonon modes (vibrational cooling). Through combinations of plasma-generated supercontinuum infrared, tunable near- and mid-infrared, and terahertz pulses in pump-probe spectroscopy, supplemented with ab inito simulations, we can explore the energy transfer processes on a sub-picosecond time scale to elucidate vibrational energy transfer pathways and lifetimes in EM. Herein we highlight recent progress, including the spectral and temporal characteristics of the infrared and THz sources as well as preliminary results on select EM.

Energetic materials with different properties can be mixed or layered to control performance. However, reactions at material interfaces are poorly understood and performance may be highly dependent on the degree of mixing. In this work, we use vapor-deposited explosive multilayers as a model system to investigate shock interactions between different explosive materials with precisely controlled spacings. Samples consisted of alternating pentaerythritol tetranitrate (PETN) and hexanitrostilbene (HNS) layers, materials that have substantial differences in detonation velocity, with individual layer thicknesses in the vicinity of the critical thickness for detonation propagation of each material (~100 - 200 μm). Additional experiments on PETN/HNS bilayer samples were conducted to elucidate the role of non-ideal interfaces on detonation propagation. Preliminary hydrocode simulations were employed to simulate detonation performance, using an Arrhenius reactive burn model that was parameterized from detonation velocity and failure data from vapor-deposited films of each constituent material. Measured detonation velocities in the multilayer samples were significantly lower than expected, given that the individual PETN layer thicknesses were larger than the critical thickness for detonation propagation. The bilayer experiments highlight the role of non-ideal interfaces in contributing to this result.

Within the energetics community, considerable effort is being put forth to find a robust scale-bridging link between unreacted material microstructures and the observed material responses, e.g. detonation and sub-detonative phenomena. Specifically, one area where this scale-bridging capability is needed is mesoscale modeling of explosives initiation (MMEI); here, material microstructures are imported directly or as statistical reconstructions into a hydrocode. While MMEI is attractive for simulating the shock initiation process with ever-increasing model fidelity, a large gap remains between the data being generated at the mesoscale and the calibration of burn model parameters. In this work, stochastic burn models are explored as a paradigm-shift to address possible scale-bridging schemes. These stochastic, particle-based methods are similar to those used for granular and droplet-laden flows, with Langevin-type equations. Further parallels are drawn to turbulent combustion modeling and preliminary developments using probability density function (pdf) theory by Baer, Gartling, and DesJardin. In order to implement these new scale-bridging schemes, one example of a stochastic burn model is explained in greater detail. Results from the stochastic burn model and MMEI simulations are given to illustrate the proposed approach. Ultimately, the execution of this work will be a community endeavor; to achieve such a capability, research efforts should focus on full-field data mining and pdf evolution, in addition to new numerical techniques for hydrocodes.

We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg2SiO4), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.

An A-and B-site substitutional study of SrFeO3−δ perovskites (A’x A1−x B’y B1−y O3−δ, where A = Sr and B = Fe) was performed for a two-step solar thermochemical air separation cycle. The cycle steps encompass (1) the thermal reduction of A’x Sr1−x B’y Fe1−y O3−δ driven by concentrated solar irradiation and (2) the oxidation of A’x Sr1−x B’y Fe1−y O3−δ in air to remove O2, leaving N2 . The oxidized A’x Sr1−x B’y Fe1−y O3−δ is recycled back to the first step to complete the cycle, resulting in the separation of N2 from air and concentrated solar irradiation. A-site substitution fractions between 0 ≤ x ≤ 0.2 were examined for A’ = Ba, Ca, and La. B-site substitution fractions between 0 ≤ y ≤ 0.2 were examined for B’ = Cr, Cu, Co, and Mn. Samples were prepared with a modified Pechini method and characterized with X-ray diffractometry. The mass changes and deviations from stoichiometry were evaluated with thermogravimetry in three screenings with temperature-and O2 pressure-swings between 573 and 1473 K and 20% O2 /Ar and 100% Ar at 1 bar, respectively. A’ = Ba or La and B’ = Co resulted in the most improved redox capacities amongst temperature-and O2 pressure-swing experiments.

The Extended History Variable Reactive Burn model (XHVRB), as proposed by Starkenburg, uses shock capturing rather than current pressure for calculating the pseudo-entropy that is used to model the reaction rate of detonating explosives. In addition to its extended capabilities for modeling explosive desensitization in multi-shock environments, XHVRB's shock capturing offers potential improvement for single shock modeling over the historically used workhorse model HVRB in CTH, an Eulerian shock physics code developed at Sandia National Labs. The detailed transition to detonation of PBX9501, as revealed by embedded gauge data, is compared to models with both HVRB and XHVRB. Improvements to the comparison of model to test data are shown with XHVRB, though not all of the details of the transition are captured by these single-rate models.

Energetic materials with different properties can be mixed or layered to control performance. However, reactions at material interfaces are poorly understood and performance may be highly dependent on the degree of mixing. In this work, we use vapor-deposited explosive multilayers as a model system to investigate shock interactions between different explosive materials with precisely controlled spacings. Samples consisted of alternating pentaerythritol tetranitrate (PETN) and hexanitrostilbene (HNS) layers, materials that have substantial differences in detonation velocity, with individual layer thicknesses in the vicinity of the critical thickness for detonation propagation of each material (~100 - 200 μm). Additional experiments on PETN/HNS bilayer samples were conducted to elucidate the role of non-ideal interfaces on detonation propagation. Preliminary hydrocode simulations were employed to simulate detonation performance, using an Arrhenius reactive burn model that was parameterized from detonation velocity and failure data from vapor-deposited films of each constituent material. Measured detonation velocities in the multilayer samples were significantly lower than expected, given that the individual PETN layer thicknesses were larger than the critical thickness for detonation propagation. The bilayer experiments highlight the role of non-ideal interfaces in contributing to this result.

Signal arrival-time estimation plays a critical role in a variety of downstream seismic analyses, including location estimation and source characterization. Any arrival-time errors propagate through subsequent data-processing results. In this article, we detail a general framework for refining estimated seismic signal arrival times along with full estimation of their associated uncertainty. Using the standard short-term average/long-term average threshold algorithm to identify a search window, we demonstrate how to refine the pick estimate through two different approaches. In both cases, new waveform realizations are generated through bootstrap algorithms to produce full a posteriori estimates of uncertainty of onset arrival time of the seismic signal. The onset arrival uncertainty estimates provide additional data-derived information from the signal and have the potential to influence seismic analysis along several fronts.

Energy storage systems (ESSs) are being deployed widely due to numerous benefits including operational flexibility, high ramping capability, and decreasing costs. This study investigates the economic benefits provided by battery ESSs when they are deployed for market-related applications, considering the battery degradation cost. A comprehensive investment planning framework is presented, which estimates the maximum revenue that the ESS can generate over its lifetime and provides the necessary tools to investors for aiding the decision making process regarding an ESS project. The applications chosen for this study are energy arbitrage and frequency regulation. Lithium-ion batteries are considered due to their wide popularity arising from high efficiency, high energy density, and declining costs. A new degradation cost model based on energy throughput and cycle count is developed for Lithium-ion batteries participating in electricity markets. The lifetime revenue of ESS is calculated considering battery degradation and a cost-benefit analysis is performed to provide investors with an estimate of the net present value, return on investment and payback period. The effect of considering the degradation cost on the estimated revenue is also studied. The proposed approach is demonstrated on the IEEE Reliability Test System and historical data from PJM Interconnection.

Understanding microstructural and strain evolutions induced by noble gas production in the nuclear fuel matrix or plasma-facing materials is crucial for designing next generation nuclear reactors, as they are responsible for volumetric swelling and catastrophic failure. We describe a multimodal approach combining synchrotron-based nanoscale X-ray imaging techniques with atomic-scale electron microscopy techniques for mapping chemical composition, morphology and lattice distortion in a single crystal W induced by Kr irradiation. We report that Kr-irradiated single crystal W undergoes surface deformation, forming Kr containing cavities. Furthermore, positive strain fields are observed in Kr-irradiated regions, which lead to compression of underlying W matrix.

Supported by the U.S. Department of Energy (DOE), the Energy Exascale Earth System Model (E3SM) project aims to optimize the use of DOE resources to address the grand challenge of actionable predictions of Earth system variability and change. This requires sustained advancement to (1) integrate model development with leading-edge computational advances toward ultra-high-resolution modeling; (2) represent the coupled human-Earth system to address energy sector vulnerability to variability and change; and (3) address uncertainty in model simulations and projections. Scientific development of the E3SM modeling system is driven by the simulation requirements in three overarching science areas centering on understanding the Earth's water cycle, biogeochemistry, and cryosphere systems and their future changes. This paper serves as an introduction to the E3SM special collection, which includes 50 papers published in several AGU journals. It provides an overview of the E3SM project, including its goals and science drivers. It also provides a brief history of the development of E3SM version 1 and highlights some key findings from papers included in the special collection.

Work performed under this one-year LDRD was concerned with estimating resource requirements for small quantum test beds that are expected to be available in the near future. This work represents a preliminary demonstration of our ability to leverage quantum hardware for solving small quantum simulation problems in areas of interest to the DOE. The algorithms enabling such studies are hybrid quantum-classical variational algorithms, in particular the widely-used variational quantum eigensolver (VQE). Employing this hybrid algorithm, in which the quantum computer complements the classical one, we implemented an end-to-end application-level toolchain that allows the user to specify a molecule of interest and compute the ground state energy using the VQE approach. We found significant limitations attributable to the classical portion of the hybrid system, including a greater than greater-than-quartic power scaling of the classical memory requirements with the system size. Current VQE approaches would require an exascale machine for solving any molecule with size greater than 150 nuclei. Our findings include several improvements that we implemented into the VQE toolchain, including a new classical optimizer that is decades old but hadn't been considered before in the VQE ecosystem. Our findings suggest limitations to variational hybrid approaches to simulation that further motivate the need for a gate-based fault-tolerant quantum processor that can implement larger problems using the fully digital quantum phase estimation algorithm.

An organic glass scintillator developed by Sandia National Laboratories was characterized in terms of its light output and pulse shape discrimination (PSD) properties and compared to commercial liquid (EJ-309) and plastic (EJ-276) organic scintillators. The electron light output was determined through relative comparison of the 137Cs Compton edge location. The proton light yield was measured using a double time-of-flight technique at the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory. Using a tunable broad-spectrum neutron source and an array of pulse-shape-discriminating observation scintillators, a continuous measurement of the proton light yield was performed for EJ-309 (200 keV–3.2 MeV), EJ-276 (170 keV–4.9 MeV), and the organic glass (50 keV–20 MeV). Finally, the PSD properties of the organic glass, EJ-309, and EJ-276 were evaluated using an AmBe source and compared via a figure-of-merit metric. The organic glass exhibited a higher electron light output than both EJ-309 and EJ-276. Its proton light yield and PSD performance were comparable to EJ-309 and superior to that of EJ-276. With these performance characteristics, the organic glass scintillator is well poised to replace current state-of-the-art PSD-capable scintillators in a range of fast neutron detection applications.

Conathane EN-7 (referred to as EN-7) has been used for decades to pot electrical connectors, providing mechanical support for solder joints in cables. Unfortunately, the EN-7 formulation contains a suspect carcinogen and chemical sensitizer, toluene diisocyanate (TDI). Because of this, various groups have been formulating replacement materials, but all have come up short in final properties or in processing. We propose Arathane 5753 HVB as a replacement for EN-7. The properties compare very well with EN-7 and the processing has both advantages and disadvantages over EN-7 as discussed in detail below.

This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported classical phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.

Atomic Force Microscopy (AFM), in conjunction with Peak Force Kelvin Probe Force Microscopy (PF-KPFM) and Peak Force Scanning Spreading Resistance Microscopy (PF-SSRM), was used to assess changes on thin metal films that underwent accelerated aging. The AFM technique provides a relatively easy, non-destructive methodology that does not require high-vacuum facilities to obtain nanometer-scale spatial resolution of surface chemistry changes. Surface morphology, roughness, contact potential difference, and spreading resistance were monitored to qualitatively identify effects of aging-morphology changes and oxidation of Au, Al, Cu thin film standards as well as diffusion of CuAu and AlAu thin film stacks at 65C under dried nitrogen flow conditions. AFM PF-KPFM and PF-SSRM modes have been exercised, refined and have proven to be viable and necessary early aging detection tools.

Secondary networks are used to supply power to loads that require very high reliability and are in use around the world, particularly in dense urban downtown areas. The protection of secondary network systems poses unique challenges. The addition of distributed energy resources (DERs) to secondary networks can compound these challenges, and deployment of microgrids on secondary networks will create a new set of challenges and opportunities. This report discusses secondary networks and their protection, the challenges associated with interconnecting DERs to a secondary network, issues expected to be associated with creating microgrids on secondary networks, standards that deal with these challenges and issues, and suggestions for research and development foci that would yield new means for addressing these challenges.

There are numerous self-shielded research irradiators used in various facilities throughout the United States. The irradiators employ radioactive sources containing either ¹³⁷Cs or ⁶⁰Co and the irradiators are used for a variety of radiobiological investigations involving cellular and/or animal models. A report from the National Academy of Sciences described security issues associated with particular radiation sources and the desire for their replacement with suitable X-ray irradiators. One possible replacement would be a 320 kV X-ray irradiator. The participants in this research successfully performed in vivo radiobiological studies involving mice exposed to filtered (HVL ≈ 4 mm Cu) 320 kV X rays. Two publications (one published and one submitted at the publishing of this report) documenting key findings are provided in Appendices A and B of this report. The 320 kV X rays were found suitable for in vivo (in mice) cell survival studies and are expected to be suitable for bone marrow transplantation studies using mice but this needs to be experimentally validated.

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4–6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron–ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. Here, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments.

Wide-area time-synchronized measurements have recently revealed troublesome forced oscillations (FOs) within modern synchronized power grids. In some cases, these FOs represent a dangerous hazard to the system. Recent research has focused on locating the source of FOs to provide operators with knowledge for mitigating their impact locally. This paper presents a complementary mitigation strategy, which is to purposely induce a second oscillation into the grid which cancels the impact of the FO. Such a strategy is complementary in that it may provide valuable time to operators attempting to locate the FOs source and to determine how to rectify it. This paper presents a suppression control strategy which modulates controllable devices to automatically cancel the impact of the FO without the need for locating the source of the original FO. The strategy is based upon tuned feedback control. The approach is demonstrated on a simulation system via modulation of inverter-based generation.

Lean operation of Spark-Ignition engines can provide higher thermal efficiency compared to standard stoichiometric operation. However, for a homogeneous lean mixture, the associated reduction of flame speeds becomes an important issue from the perspective of robust ignition and fast flame spread throughout the charge. This study is focused on the use of a lean partial fuel stratification strategy that can stabilize the deflagration, while sufficiently fast combustion is ensured via the use of end-gas autoignition. The engine has a spray-guided Direct-Injection Spark-Ignition combustion system and was fueled with either a high-octane certification gasoline or E85. Partial fuel stratification was achieved using several fuel injections during the intake stroke in combination with a small pilot-injection concurrent with the Spark-Ignition. The results reveal that partial fuel stratification enables very stable combustion, offering higher thermal efficiency for parts of the load range in comparison to well-mixed lean and stoichiometric combustion. The heat release and flame imaging demonstrate that the combustion often has three distinct stages. The combustion of the pilot-injected fuel, ignited by the normal spark, acts as a “super igniter,” ensuring a very repeatable initiation of combustion, and flame incandescence reveals locally rich conditions. The second stage is mainly composed of blue flame propagation in a well-mixed lean mixture. The third stage is the compression autoignition of a well-mixed and typically very lean end-gas. The end-gas autoignition is critical for achieving high combustion efficiency, high thermal efficiency, and stable combustion. Partial fuel stratification enables very effective combustion-phasing control, which is critical for controlling the occurrence and intensity of end-gas autoignition. Comparing the gasoline and E85 fuels, it is noted that achieving end-gas autoignition for the higher octane E85 requires a more aggressive compression of the end-gas via the use of a more advanced combustion phasing or higher intake-air temperature.

We develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

Graph coloring is often used in parallelizing scientific computations that run in distributed and multi-GPU environments; it identifies sets of independent data that can be updated in parallel. Many algorithms exist for graph coloring on a single GPU or in distributed memory, but hybrid MPI+GPU algorithms have been unexplored until this work, to the best of our knowledge. We present several MPI+GPU coloring approaches that use implementations of the distributed coloring algorithms of Gebremedhin et al. and the shared-memory algorithms of Deveci et al. The on-node parallel coloring uses implementations in KokkosKernels, which provide parallelization for both multicore CPUs and GPUs. We further extend our approaches to solve for distance-2 coloring, giving the first known distributed and multi-GPU algorithm for this problem. In addition, we propose novel methods to reduce communication in distributed graph coloring. Our experiments show that our approaches operate efficiently on inputs too large to fit on a single GPU and scale up to graphs with 76.7 billion edges running on 128 GPUs.

Benefits of local recovery (restarting only a failed process or task) have been previously demonstrated in parallel solvers. Local recovery has a reduced impact on application performance due to masking of failure delays (for message-passing codes) or dynamic load balancing (for asynchronous many-task codes). In this paper, we implement MPI-process-local checkpointing and recovery of data (as an extension of the Fenix library) in combination with an existing method for local detection of silent errors in partial-differential-equation solvers, to show a path for incorporating lightweight silent-error resilience. In addition, we demonstrate how asynchrony introduced by maximizing computation-communication overlap can halt the propagation of delays. For a prototype stencil solver (including an iterative-solver-like variant) with injected memory bit flips, results show greatly reduced overhead under weak scaling compared to global recovery, and high failure-masking efficiency. The approach is expected to be generalizable to other MPI-based solvers.

Exceptions and errors occurring within mission critical applications due to hardware failures have a high cost. With the emerging Next Generation Platforms (NGPs), the rate of hardware failures will likely increase. Therefore, designing our applications to be resilient is a critical concern in order to retain the reliability of results while meeting the constraints on power budgets. In this paper, we discuss software resilience in AMTs at both local and distributed scale. We choose HPX to prototype our resiliency designs. We implement two resiliency APIs that we expose to the application developers, namely task replication and task replay. Task replication repeats a task n-times and executes them asynchronously. Task replay reschedules a task up to n-times until a valid output is returned. Furthermore, we expose algorithm based fault tolerance (ABFT) using user provided predicates (e.g., checksums) to validate the returned results. We benchmark the resiliency scheme for both synthetic and real world applications at local and distributed scale and show that most of the added execution time arises from the replay, replication or data movement of the tasks and not the boilerplate code added to achieve resilience.

We present the execution model of Virtual Transport (VT) a new, Asynchronous Many-Task (AMT) runtime system that provides unprecedented integration and interoperability with MPI. We have developed VT in conjunction with large production applications to provide a highly incremental, high-value path to AMT adoption in the dominant ecosystem of MPI applications, libraries, and developers. Our aim is that the'MPI+X' model of hybrid parallelism can smoothly extend to become'MPI+VT +X'. We illustrate a set of design and implementation techniques that have been useful in building VT. We believe that these ideas and the code embodying them will be useful to others building similar systems, and perhaps provide insight to how MPI might evolve to better support them. We motivate our approach with two applications that are adopting VT and have begun to benefit from increased asynchrony and dynamic load balancing.

As network speeds increase, the overhead of processing incoming messages is becoming onerous enough that many manufacturers now provide network interface cards (NICs) with offload capabilities to handle these overheads. This increase in NIC capabilities creates an opportunity to enable computation on data in-situ on the NIC. These enhanced NICs can be classified into several different categories of SmartNICs. SmartNICs present an interesting opportunity for future runtime software designs. Designing runtime software to be located in the network as opposed to the host level leads to new radical distributed runtime possibilities that were not practical prior to SmartNICs. In the process of transitioning to a radically different runtime software design for SmartNICs there are intermediary steps of migrating current runtime software to be offloaded onto a SmartNIC that also present interesting possibilities. This paper will describe SmartNIC design and how SmartNICs can be leveraged to offload current generation runtime software and lead to future radically different in-network distributed runtime systems.

Exceptions and errors occurring within mission critical applications due to hardware failures have a high cost. With the emerging Next Generation Platforms (NGPs), the rate of hardware failures will likely increase. Therefore, designing our applications to be resilient is a critical concern in order to retain the reliability of results while meeting the constraints on power budgets. In this paper, we discuss software resilience in AMTs at both local and distributed scale. We choose HPX to prototype our resiliency designs. We implement two resiliency APIs that we expose to the application developers, namely task replication and task replay. Task replication repeats a task n-times and executes them asynchronously. Task replay reschedules a task up to n-times until a valid output is returned. Furthermore, we expose algorithm based fault tolerance (ABFT) using user provided predicates (e.g., checksums) to validate the returned results. We benchmark the resiliency scheme for both synthetic and real world applications at local and distributed scale and show that most of the added execution time arises from the replay, replication or data movement of the tasks and not the boilerplate code added to achieve resilience.

Multithreaded MPI applications are gaining popularity in scientific and high-performance computing. While the combination of programming models is suited to support current parallel hardware, it moves threading models and their interaction with MPI into focus. With the advent of new threading libraries, the flexibility to select threading implementations of choice is becoming an important usability feature. Open MPI has traditionally avoided componentizing its threading model, relying on code inlining and static initialization to minimize potential impacts on runtime fast paths and synchronization. This paper describes the implementation of a generic threading runtime support in Open MPI using the Opal Modular Component Architecture. This architecture allows the programmer to select a threading library at compile-or run-time, providing both static initialization of threading primitives as well as dynamic instantiation of threading objects. In this work, we present the implementation, define required interfaces, and discuss trade-offs of dynamic and static initialization.

Understanding the fundamental limits of gas deliverable capacity in porous materials is of critical importance as it informs whether technical targets (e.g., for on-board vehicular storage) are feasible. High-throughput screening studies of rigid materials, for example, have shown they are not able to achieve the original ARPA-E methane storage targets, yet an interesting question remains: what is the upper limit of deliverable capacity in flexible materials? In this work we develop a statistical adsorption model that specifically probes the limit of deliverable capacity in intrinsically flexible materials. The resulting adsorption thermodynamics indicate that a perfectly designed, intrinsically flexible nanoporous material could achieve higher methane deliverable capacity than the best benchmark systems known to date with little to no total volume change. Density functional theory and grand canonical Monte Carlo simulations identify a known metal-organic framework (MOF) that validates key features of the model. Therefore, this work (1) motivates a continued, extensive effort to rationally design a porous material analogous to the adsorption model and (2) calls for continued discovery of additional high deliverable capacity materials that remain hidden from rigid structure screening studies due to nominal non-porosity.

Exceptions and errors occurring within mission critical applications due to hardware failures have a high cost. With the emerging Next Generation Platforms (NGPs), the rate of hardware failures will likely increase. Therefore, designing our applications to be resilient is a critical concern in order to retain the reliability of results while meeting the constraints on power budgets. In this paper, we discuss software resilience in AMTs at both local and distributed scale. We choose HPX to prototype our resiliency designs. We implement two resiliency APIs that we expose to the application developers, namely task replication and task replay. Task replication repeats a task n-times and executes them asynchronously. Task replay reschedules a task up to n-times until a valid output is returned. Furthermore, we expose algorithm based fault tolerance (ABFT) using user provided predicates (e.g., checksums) to validate the returned results. We benchmark the resiliency scheme for both synthetic and real world applications at local and distributed scale and show that most of the added execution time arises from the replay, replication or data movement of the tasks and not the boilerplate code added to achieve resilience.

Benefits of local recovery (restarting only a failed process or task) have been previously demonstrated in parallel solvers. Local recovery has a reduced impact on application performance due to masking of failure delays (for message-passing codes) or dynamic load balancing (for asynchronous many-task codes). In this paper, we implement MPI-process-local checkpointing and recovery of data (as an extension of the Fenix library) in combination with an existing method for local detection of silent errors in partial-differential-equation solvers, to show a path for incorporating lightweight silent-error resilience. In addition, we demonstrate how asynchrony introduced by maximizing computation-communication overlap can halt the propagation of delays. For a prototype stencil solver (including an iterative-solver-like variant) with injected memory bit flips, results show greatly reduced overhead under weak scaling compared to global recovery, and high failure-masking efficiency. The approach is expected to be generalizable to other MPI-based solvers.

Many-core systems are beginning to feature novel large, high-bandwidth intermediate memory as a visible part of the memory hierarchy. This paper discusses how to make use of intermediate memory when composing matrix multiply with transpose to compute $A$ * AT. We re-purpose the cache-oblivious approach developed by Frigo et al. and apply it to the composition of a bandwidth-bound kernel (transpose) with a compute-bound kernel (matrix multiply). Particular focus is on regions of matrix shapes far from square that are not usually considered. Our codes are simpler than optimized codes, but reasonably close in performance. Also, perhaps of more importance is developing a paradigm for how to construct other codes using intermediate memories.

A system’s response to disturbances in an internal or external driving signal can be characterized as performing an implicit computation, where the dynamics of the system are a manifestation of its new state holding some memory about those disturbances. Identifying small disturbances in the response signal requires detailed information about the dynamics of the inputs, which can be challenging. This paper presents a new method called the Information Impulse Function (IIF) for detecting and time-localizing small disturbances in system response data. The novelty of IIF is its ability to measure relative information content without using Boltzmann’s equation by modeling signal transmission as a series of dissipative steps. Since a detailed expression of the informational structure in the signal is achieved with IIF, it is ideal for detecting disturbances in the response signal, i.e., the system dynamics. Those findings are based on numerical studies of the topological structure of the dynamics of a nonlinear system due to perturbated driving signals. The IIF is compared to both the Permutation entropy and Shannon entropy to demonstrate its entropy-like relationship with system state and its degree of sensitivity to perturbations in a driving signal.

In this synthesis, we assess present research and anticipate future development needs in modeling water quality in watersheds. We first discuss areas of potential improvement in the representation of freshwater systems pertaining to water quality, including representation of environmental interfaces, in-stream water quality and process interactions, soil health and land management, and (peri-)urban areas. In addition, we provide insights into the contemporary challenges in the practices of watershed water quality modeling, including quality control of monitoring data, model parameterization and calibration, uncertainty management, scale mismatches, and provisioning of modeling tools. Finally, we make three recommendations to provide a path forward for improving watershed water quality modeling science, infrastructure, and practices. These include building stronger collaborations between experimentalists and modelers, bridging gaps between modelers and stakeholders, and cultivating and applying procedural knowledge to better govern and support water quality modeling processes within organizations.

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.

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.

Abstract not provided.

Abstract not provided.

A High-Altitude Electromagnetic Pulse (HEMP) is a potential threat to the power grid. HEMP can couple to transmission lines and cables, causing significant overvoltages which can be harmful to line connected equipment. The effects of overvoltages on various types of power systems components need to be understood. HEMP effects on trip coils were tested and presented in this report. A high voltage pulser was built to replicate the induced voltage waveform from a HEMP. The pulser was used to test breaker trip coils with increasing pulse magnitudes ranging from 20 kV to 80 kV. The State-of-Health of each trip coils was measured via mechanical operation and impedance measurements before and after each insult to identify any damage or degradation to the trip coils. Dielectric breakdown was observed at the conductor leads during testing, causing the HEMP insult to be diverted to the grounded casing. However, the dielectric breakdown did not cause interference with regular device operation.

Abstract not provided.

Abstract not provided.

Research shows that individuals often overestimate their knowledge and performance without realizing they have done so, which can lead to faulty technical outcomes. This phenomenon is known as the Dunning-Kruger effect (Kruger & Dunning, 1999). This research sought to determine if some individuals were more prone to overestimating their performance due to underlying personality and cognitive characteristics. To test our hypothesis, we first collected individual difference measures. Next, we asked participants to estimate their performance on three performance tasks to assess the likelihood of overestimation. We found that some individuals may be more prone to overestimating their performance than others, and that faulty problem-solving abilities and low skill may be to blame. Encouraging individuals to think critically through all options and to consult with others before making a high-consequence decision may reduce overestimation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A useful and popular waveform for high-performance radar systems is the Linear Frequency Modulated (LFM) chirp. The chirp may have a positive frequency slope with time (up-chirp) or a negative frequency slope with time (down-chirp). There is no inherent advantage to one with respect to the other, except that the receiver needs to be matched to the proper waveform. However, if up-chirps and down-chirps are employed on different pulses in the same Coherent Processing Interval (CPI), then care must be taken to maintain coherence in the range-compressed echo signals. We present the mathematics for doing so, for both correlation processing and stretch processing.

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.

We present an implementation that can keep a coldatom ensemble within a sub-millimeter diameter hole in a transparent membrane. Based on the effective beam diameter of the magneto-optical trap (MOT), d = 400 mm-hole diameter, we measure the atom number that is 105 times higher than the predicted value using the conventional d6 scaling rule. Atoms trapped by the membrane MOT are cooled down to 10 mK with sub- Doppler cooling process and can be potentially coupled to the photonic/electronic integrated circuits that can be fabricated in the membrane device by taking a step toward the atom trap integrated platform.

Abstract not provided.

Future nuclear fuel cycle facilities will see a significant benefit from considering materials accountancy requirements early in the design process. The Material Protection, Accounting, and Control Technologies (MPACT) working group is demonstrating Safeguards and Security by Design (SSBD) for a notional electrochemical reprocessing facility as part of a 2020 Milestone. The idea behind SSBD is to consider regulatory requirements early in the design process to provide more optimized systems and avoid costly retrofits later in the design process. Safeguards modeling, using single analyst tools, allows the designer to efficiently consider materials accountancy approaches that meet regulatory requirements. However, safeguards modeling also allows the facility designer to go beyond current regulations and work toward accountancy designs with rapid response and lower thresholds for detection of anomalies. This type of modeling enables new safeguards approaches and may inform future regulatory changes. The Separation and Safeguards Performance Model (SSPM) has been used for materials accountancy system design and analysis. This paper steps through the process of designing a Material Control and Accountancy (MC&A) system, presents the baseline system design for an electrochemical reprocessing facility, and provides performance metrics from the modeling analysis. The most critical measurements in the electrochemical facility are the spent fuel input, electrorefiner salt, and U/TRU product output measurements. Finally, material loss scenario analysis found that measurement uncertainties (relative standard deviations) for Pu would need to be at 1% (random and systematic error components) or better in order to meet domestic detection goals or as high as 3% in order to meet international detection goals, based on a 100 metric ton per year plant size.

Abstract not provided.

The Sodium-Cooled Fast Reactor (SFR) system was identified during the Generation IV Technology Roadmap as a promising technology to perform the actinide management mission and, if enhanced economics for the system could be realized, also the electricity and heat production missions. The main characteristics of the SFR that make it especially suitable for the actinide management mission are: Consumption of transuranics in a closed fuel cycle, thus reducing the radiotoxicity and heat load which facilitates waste disposal and geologic isolation; Enhanced utilization of uranium resources through efficient management of fissile materials and multi-recycle; and, High level of safety achieved through inherent and passive means that accommodate transients and bounding events with significant safety margins.

Abstract not provided.

Abstract not provided.

Abstract not provided.