Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.
Coupled Mode Theory (CMT) is a classic model that addresses many diverse problems in physics and electrical engineering. Although CMT is well-established in areas such as directional coupler design, its significance in antenna engineering is less well-recognized. Recently, Characteristic Mode Analysis (CMA) has shown how CMT quantitatively models a wide variety of multi-mode patch structures. In this paper, we extend these ideas to a double-tuned slot antenna utilizing some unique features of the method of moments code FEKO. The modeled slot antenna has a stable radiation pattern throughout its 50% impedance bandwidth. CMA decomposition indicates the antenna consists of two resonators: the slot and a lumped LC resonator that represents the feed probe. The two are electromagnetically coupled to produce broad impedance bandwidth. An equivalent circuit whose values are derived directly from the CMA results is given.
MTTKRP is the bottleneck operation in algorithms used to compute the CP tensor decomposition. For sparse tensors, utilizing the compressed sparse fibers (CSF) storage format and the CSF-oriented MTTKRP algorithms is important for both memory and computational efficiency on distributed-memory architectures. Existing intelligent tensor partitioning models assume the computational cost of MTTKRP to be proportional to the total number of nonzeros in the tensor. However, this is not the case for the CSF-oriented MTTKRP on distributed-memory architectures. We outline two deficiencies of nonzero-based intelligent partitioning models when CSF-oriented MTTKRP operations are performed locally: failure to encode processors' computational loads and increase in total computation due to fiber fragmentation. We focus on existing fine-grain hypergraph model and propose a novel vertex weighting scheme that enables this model encode correct computational loads of processors. We also propose to augment the fine-grain model by fiber nets for reducing the increase in total computational load via minimizing fiber fragmentation. In this way, the proposed model encodes minimizing the load of the bottleneck processor. Parallel experiments with real-world sparse tensors on up to 1024 processors prove the validity of the outlined deficiencies and demonstrate the merit of our proposed improvements in terms of parallel runtimes.
We report on progress in developing macroscopic balance equations for combustion and electrochemistry systems. A steady state solution capability is described for the macroscopic reactor network, with an associated steady state continuation method and solution storage capability added in. An example is provided of continuation of a hydrogen flame versus the equivalence ratio. The reactor modeling capability is extended to charged fluid systems, with a description of the new ChargedFluidReactor, SubstrateElement, and MetalCurrentElement reactor classes and novel setup of unknowns within these reactors that preserve charge neutrality. Zuzax's setup for electrochemistry is explained including the specification of the electron chemical potential and the adherence to the SHE Reference electrode specification. The description of the different ways to enter electrochemical reaction rates are described, contrasted, and their derivations with respect to one another are derived. An example of using the ChargedFluidReactor within corrosion problems is provided. We present a description of calculations to understand the phenomena of corrosion of copper from a micron sized droplet of NaCl water droplet, where secondary spreading occurs. An analysis of the discrepancies with experiment is carried out, demonstrating that macroscopic balances can be an important tool for understanding what major factors need to be addressed for a better understanding of a physical system.
The purpose of this work is to fit a previously developed empirical equation for puncture energy to simulation data. The conservative puncture energy equation could be used to expedite the process of performing calculations in the development of safety measures, avoiding the need to create complex finite element models for specific puncture scenarios. A total of 108 simulations are developed by varying coupon thickness, coupon material, probe shape, and probe diameter. The simulations are comprised of a low-velocity probe puncturing the coupons, from which the probe kinetic energy change is calculated. The empirical equation is fit to the dimensions, material properties, and energy results using a non-linear least-squares regression method within Python, which determines the two constant parameters for each fit. More statistically significant fit results are achieved by separating the data by probe shape and coupon material.
We discuss thinned InAsSb resonant infrared detectors that are designed to enable high quantum efficiency by using interleaved nanoantennas to read out two wavelengths from each pixel simultaneously.
The structures that surround and support optical components play a key role in the performance of the overall optical system. For aerospace applications, creating an opto-mechanical structure that is athermal, lightweight, robust, and can be quickly developed from concept through to hardware is challenging. This project demonstrates a design and fabrication method for optical structures using origami-style folded, photo-etched sheetmetal pieces that are micro-welded to each other or to 3d printed metal components. Thin flexures, critical for athermal mounting of optics, can be thinner with sheetmetal than from standard machining, which leads to more compact designs and the ability to mount smaller optics. Building a structure by starting with the thinnest features, then folding that thin material to make the ''thicker'' sections is the opposite of standard machining (cutting thin features from thicker blocks). A design method is shown with mass savings of >90%, and stiffness to weight ratio improvements of 5x to 10x compared to standard methods for space systems hardware. Designs and processes for small, flexured, actively aligned systems are demonstrated as are methods for producing lightweight, structural, Miura-core sandwich panels in both flat and curved configurations. Concepts for deployable panels and component hinges are explored, as is a lens subcell with tunable piston movement with temperature change and an ultralight sunshade.
Based on the latest DOE (Department of Energy) milestones, Sandia needs to convert to IPv6 (Internet Protocol version 6)-only networks over the next 5 years. Our original IPv6 migration plan did not include migrating to IPv6-only networks at any point within the next 10 years, so it must necessarily change. To be successful in this endeavor, we need to evaluate technologies that will enable us to deploy IPv6-only networks early without creating system stability or security issues. We have set up a test environment using technology representative of our production network where we configured and evaluated industry standard translation technologies and techniques. Based on our results, bidirectional translation between IPv4 (Internet Protocol version 4) and IPv6 is achievable with our current equipment, but due to the complexity of the configuration, may not scale well to our production environment.
Several Department of Energy (DOE) facilities have materials stored in robust, welded, stainless - steel containers with presumed fire - induced pressure response behaviors. Lack of test data related to fire exposure requires conservative safety analysis assumptions for container response at these facilities. This conservatism can in turn result in the implementation of challenging operational restrictions with costly nuclear safety controls. To help address this issue for sites that store DOE 3013 stainless - steel containers, a series of ten tests were undertaken at Sandia National Laboratories. The goal of this test series was to obtain the response behavior for various configurations of DOE 3013 containers with various payload compositions when exposed to one of two ASTM fire conditions. Key parameters measured in the test series included identification of failure - specific characteristics such as pressure, temperature, and whether or not a vessel was breached during a test . Numerous failure - specific characteristics were identified from the ten tests. This report describes the implementation and execution of the test series performed to identify these failure - specific characteristics. Discussions on the test configurations, payload compositions, thermal insults, and experimental setups are presented. Test results in terms of pressurization and vessel breach (or no - breach) are presented along with corresponding discussions for each test.
Transition metal dichalcogenides, such as MoS2, are known to undergo a structural phase transformation as well as a change in the electronic conductivity upon Li intercalation. These properties make them candidates for charge-tunable ion-insertion materials that could be used in electrochemical devices. In this paper, we study the phase stability and electronic structure of H and T' Li-intercalated MoX2 bilayers with X = S, Se, or Te. Using first-principles calculations in combination with classical and machine learning modeling approaches, we find that the H phase is more stable at low Li concentration for all X and the critical Li concentration at which the T'→H transition occurs decreases with increasing mass of X. Furthermore, the relative free energy of the two phases becomes less sensitive to Li insertion with increasing atomic mass of the chalcogen atom X. While the electronic conductivity increases with increasing ion concentration at low concentrations, we do not observe a (positive) conductivity jump at the phase transition from H to T'.
Independent gamma spectroscopy data analysis of a plutonium oxide sample was requested on July 23, 2021. The primary request was to assess the Pu-239 activity/mass of the sample using previously collected gamma spectral data. Using the provided gamma spectral analysis report and spectral files, an independent evaluation of the data was conducted without any prior knowledge of the isotopic activity/mass of the sample.
A new Boltzmann-CSD solver has been developed within the SCEPTRE radiation-transport code, based on the 1st-order form of the transport equation, using discontinuous finite elements in space and energy and discrete ordinates in angle. The Boltzmann-CSD solver has been validated against experimental data for electron energy deposition distributions and for electron emission spectra. Comparison of the calculated results with experimental data shows excellent agreement for many of the test configurations and reasonable agreement for other test configurations. The tests have also been modeled with the ITS Monte Carlo code, which also shows excellent to reasonable agreement with the SCEPTRE results and experimental data. The SCEPTRE Boltzmann-CSD solver relies on electron cross sections generated by the legacy CEPXS code, which currently is limited to electron-only Boltzmann-CSD cross sections. Performing full electron-photon radiation transport with the Boltzmann-CSD solver will require further development in the cross section generating code. For the energy- deposition calculations, neglecting photon transport results in at most about 5% overprediction of the energy deposition for high-energy electrons on high-Z targets, and relatively insignificant difference for the other test configurations.
With the growing number of applications designed for heterogeneous HPC devices, application programmers and users are finding it challenging to compose scalable workflows as ensembles of these applications, that are portable, performant and resilient. The Kokkos C++ library has been designed to simplify this cumbersome procedure by providing an intra-application uniform programming model and portable performance. However, assembling multiple Kokkos-enabled applications into a complex workflow is still a challenge. Although Kokkos enables a uniform programming model, the inter-application data exchange still remains a challenge from both performance and software development cost perspectives. In order to address this issue, we propose Kokkos data staging memory space, an extension of Kokkos' data abstraction (memory space) for heterogeneous computing systems. This new abstraction allows to express data on a virtual shared-space for multiple Kokkos applications, thus extending Kokkos to support inter-application data exchange to build an efficient application workflow. Additionally, we study the effectiveness of asynchronous data layout conversions for applications requiring different memory access patterns for the shared data. Our preliminary evaluation with a synthetic benchmark indicate the effectiveness of this conversion adapted to three different scenarios representing access frequency and use patterns of the shared data.
Sub-channel codes are one of the the modeling and simulation tools used for thermal-hydraulic analysis of nuclear reactors. A few examples of such sub-channel codes are the COolant Boiling in Rod Arrays (COBRA) family of codes. The approximations that are used to simplify the fluid conservation equations into sub-channel form, mainly that of axially-dominated flow, lead to noticeable limitations on sub-channels solvers for problems with significant flow in lateral directions. In this report, a two-dimensional Cartesian solver is developed and implemented within CTF-R, which is the residual solver in the North Carolina State University version of COBRA-TF (CTF). The new solver will enable CTF to simulate flow that is not axially-dominated. The appropriate Cartesian forms of the conservation equations are derived and implemented in the solver. Once the conservation equations are established, the process of constructing the matrix system was altered to solve a two-dimensional staggered grid system. A simple case was used to test that the two-dimensional Cartesian solver is accurate. The test problem does not include any source terms or flow in the lateral direction. The results show that the solver was able to run the simple case and converge to a steady-state solution. Future work will focus on testing existing capabilities by using test cases that include transients and equation cross-terms. Future work will also include adding additional capabilities such as enabling the solver to include cases with source terms and three dimensional cases.
For a wide range of PDEs, the discontinuous Petrov–Galerkin (DPG) methodology of Demkowicz and Gopalakrishnan provides discrete stability starting from a coarse mesh and minimization of the residual in a user-controlled norm, among other appealing features. Research on DPG for transient problems has mainly focused on spacetime discretizations, which has theoretical advantages, but practical costs for computations and software implementations. The sole examination of time-stepping DPG formulations was performed by Führer, Heuer, and Gupta, who applied Rothe's method to an ultraweak formulation of the heat equation to develop an implicit time-stepping scheme; their work emphasized theoretical results, including error estimates in time and space. In the present work, we follow Führer, Heuer, and Gupta in examining the heat equation; our focus is on numerical experiments, examining the stability and accuracy of several formulations, including primal as well as ultraweak, and explicit as well as implicit and Crank–Nicolson time-stepping schemes. We are additionally interested in communication-avoiding algorithms, and we therefore include a highly experimental formulation that places all the trace terms on the right-hand side of the equation.
We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.
High-temperature particle receivers are being pursued to enable next-generation concentrating solar thermal power (CSP) systems that can achieve higher temperatures (> 700 °C) to enable more efficient power cycles, lower overall system costs, and emerging CSP-based process-heat applications. The objective of this work was to develop characterization methods to quantify the particle and heat losses from the open aperture of the particle receiver. Novel camera- based imaging methods were developed and applied to both laboratory-scale and larger 1 MWt on-sun tests at the National Solar Thermal Test Facility in Albuquerque, New Mexico. Validation of the imaging methods was performed using gravimetric and calorimetric methods. In addition, conventional particle-sampling methods using volumetric particle-air samplers were applied to the on-sun tests to compare particle emission rates with regulatory standards for worker safety and pollution. Novel particle sampling methods using 3-D printed tipping buckets and tethered balloons were also developed and applied to the on-sun particle-receiver tests. Finally, models were developed to simulate the impact of particle size and wind on particle emissions and concentrations as a function of location. Results showed that particle emissions and concentrations were well below regulatory standards for worker safety and pollution. In addition, estimated particle temperatures and advective heat losses from the camera-based imaging methods correlated well with measured values during the on-sun tests.
A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. The model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.
This document is a guide to setting the system and processing configuration for the Geophysical Monitoring System (GMS) Station State-of-Health (SOH) Monitoring and Interactive Analysis (IAN) applications.
Stress corrosion cracking (SCC) is an important failure degradation mechanism for storage of spent nuclear fuel. Since 2014, Sandia National Laboratories has been developing a probabilistic methodology for predicting SCC. The model is intended to provide qualitative assessment of data needs, model sensitivities, and future model development. In fiscal year 2021, improvement of the SCC model focused on the salt deposition, maximum pit size, and crack growth rate models.
Fungi produce and excrete various proteins, enzymes, and polysaccharides, which may be used for the synthesis of nanoparticles. This study investigated the effect an anion species on the synthesis of ceramic nanoparticles using fungal filtrates. In this work, ceramic zinc oxide (ZnO) nanoparticles ranging between 1 nm and 1000 nm were successfully synthesized using three different filamentous fungi: Aspergillus sp., Penicillium sp., and Paecilomyces variotti. Each fungus was cultured, and the filtrate was extracted and individually exposed to zinc nitrate, zinc sulfate, or zinc chloride. The formation of nanoparticles was characterized using UV-visible spectrophotometry (UV-Vis), fluorescence microscopy, and with transmission electron microscopy (TEM) analyses. UV-Vis spectra exhibited broad increases in the absorption across the range of 200 nm - 800 nm, which corresponded to the formation of ZnO nanoparticles under various conditions. Nanoparticle formation was confirmed with fluorescence microscopy and TEM analysis and determined to form particles with an irregular spherical shape. To date, our work demonstrates that the ability of fungi to synthesize ZnO nanoparticles is not genus/species-specific but is dependent on the starting composition of a given metal salt.
Twin boundaries play an important role in the thermodynamics, stability, and mechanical properties of nanocrystalline metals. Understanding their structure and chemistry at the atomic scale is key to guide strategies for fabricating nanocrystalline materials with improved properties. We report an unusual segregation phenomenon at gold-doped platinum twin boundaries, which is arbitrated by the presence of disconnections, a type of interfacial line defect. By using atomistic simulations, we show that disconnections containing a stacking fault can induce an unexpected transition in the interfacial-segregation structure at the atomic scale, from a bilayer, alternating-segregation structure to a trilayer, segregation-only structure. This behavior is found for faulted disconnections of varying step heights and dislocation characters. Supported by a structural analysis and the classical Langmuir-McLean segregation model, we reveal that this phenomenon is driven by a structurally induced drop of the local pressure across the faulted disconnection accompanied by an increase in the segregation volume.
The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using relatively high backfill pressures (up to approximately 800 kPa) in the canister to enhance internal natural convection. This pressure differential offers a comparatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.89 mm (0.350 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO2, were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to demonstrate a new capability to characterize SCCs under well-controlled boundary conditions. Modeling efforts were also initiated that evaluate the depletion of aerosols in a commercial dry storage canister. These preliminary modeling and ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.
Well integrity is one of the major concerns in long-term geologic storage sites due to its potential risk for well leakage and groundwater contamination. Evaluating changes in electrical responses due to energized steel-cased wells has the potential to quantify and predict possible wellbore failures, as any kind of breakage or corrosion along highly-conductive well casings will have an impact on the distribution of subsurface electrical potential. However, realistic wellbore-geoelectrical models that can fully capture fine scale details of well completion design and the state of well damage at the field scale require extensive computational e.ort, or can even be intractable to simulate. To overcome this computational burden while still keeping the model realistic, we use the hierarchical finite element method which represents electrical conductivity at each dimensional component (1-D edges, 2-D planes and 3-D cells) of a tetrahedra mesh. This allows well completion designs with real-life geometric scales and well systems with realistic, detailed, progressive corrosion and damage in our models. Here, we present a comparison of possible discretization approaches of a multi-casing completion design in the finite-element model. The e.ects of the surface casing length and the coupling between concentric well casings, as well as the e.ects of the degree and the location of well damage on the electrical responses are also examined. Finally, we analyze real surface electric field data to detect wellbore integrity failure associated with damage.
An approach to increase the value of carbon fiber for wind turbines blades, and other compressive strength driven designs, is to identify pathways to increase its cost-specific compressive strength. A finite element model has been developed to evaluate the predictiveness of current finite element methods and to lay groundwork for future studies that focus on improving the cost-specific compressive strength. Parametric studies are conducted to understand which uncertainties in the model inputs have the greatest impact on compressive strength predictions. A statistical approach is also presented that enables the micromechanical model, which is deterministic, to efficiently account for statistical variability in the fiber misalignment present in composite materials; especially if the results from the hexagonal and square pack models are averaged. The model was found to agree well with experimental results for a Zoltek PX-35 pultrusion. The sensitivity studies suggest that the fiber packing and the interface shear strength have the greatest impact on compressive strength prediction for the fiber reinforced polymer studied here. Based on the performance of the modeling approach presented in this work, it is deemed sufficient for future work which will seek to identify carbon fiber composites with improved cost-specific compressive strength.
The commercial software package Barracuda, developed by CPFD Software for simulating particle-laden fluid flows, is evaluated as a means to simulate the motion of bubbles in vibrating liquid-filled containers. Demonstration simulations of bubbles rising due to buoyancy forces in a cylinder filled with silicone oil and angled at 0, 30, 45, and 60 degrees from the vertical were performed by CPFD Software. The results of these simulations are discussed, and the capabilities of Barracuda for simulating bubble motion are assessed. It was determined that at present Barracuda does not meet the needs of the desired application. Further developments that would enable its use for this application are highlighted.
Esteves, Giovanni E.; Lundh, James S.; Coleman, Kathleen; Song, Yiwen; Griffin, Benjamin A.; Douglas, Erica A.; Edstrand, Adam E.; Badescu, Stefan C.; Leach, Jacob H.; Moody, Baxter; Trolier-Mckinstry, Susan; Choi, Sukwon; Moore, Elizabeth A.
In this study, the Raman biaxial stress coefficients KII and strain-free phonon frequencies ω0 have been determined for the E2 (low), E2 (high), and A1 (LO) phonon modes of aluminum nitride, AlN, using both experimental and theoretical approaches. The E2 (high) mode of AlN is recommended for the residual stress analysis of AlN due to its high sensitivity and the largest signal-to-noise ratio among the studied modes. The E2 (high) Raman biaxial stress coefficient of -3.8 cm-1/GPa and strain-free phonon frequency of 656.68 cm-1 were then applied to perform both macroscopic and microscopic stress mappings. For macroscopic stress evaluation, the spatial variation of residual stress was measured across an AlN-on-Si wafer prepared by sputter deposition. A cross-wafer variation in residual stress of ∼150 MPa was observed regardless of the average stress state of the film. Microscopic stress evaluation was performed on AlN piezoelectric micromachined ultrasonic transducers (pMUTs) with submicrometer spatial resolution. These measurements were used to assess the effect of device fabrication on residual stress distribution in an individual pMUT and the effect of residual stress on the resonance frequency. At ∼20 μm directly outside the outer edge of the pMUT electrode, a large lateral spatial variation in residual stress of ∼100 MPa was measured, highlighting the impact of metallization structures on residual stress in the AlN film.
Alkali metals, such as lithium, sodium, potassium, etc., are highly reactive elements. While researchers generally handle these metals with caution, less caution is taken when these elements have been “reacted”. In this work, a recent incident is examined in which a pair of researchers ignited a lithium silicide alloy sample that was assumed to be fully hydrated to lithium hydroxide and, thereby, no longer water-reactive. However, variations in the original chemical composition of the lithium compounds examined resulted in select mixtures failing to hydrate and react completely to lithium hydroxide in the time frame allowed. This gave rise to residual unreacted, water-sensitive lithium silicide which resulted in a violent exothermic reaction with water and autoignition of the produced hydrogen gas. This Article describes this incident and improvements that can be implemented to prevent similar incidents from occurring.
The electric power grid is one of the most critical national infrastructures, and determining the susceptibility of power grid elements to external factors is of significant importance for ensuring grid resilience. Reliable energy is vital to the safety and security of society. One potential threat to the power grid comes in the form of strong electromagnetic field transients arising from high-altitude nuclear weapon detonation. The radiated EM fields from these can affect the operation of electronic components via direct field exposure or from the conducted transients that arise from coupling onto long cables. Vulnerability to these pulses for many electrical components on the grid is unknown. This research focuses on conducted pulse testing of digital protective relays in a power substation and their associated high-voltage circuit breaker circuit and instrumentation transformer circuits. The relays, yard cables, power supplies, and components representing yard equipment were assembled in a manner consistent with installation in a substation to represent the pulse's propagation in the components and wiring. Equipment was tested using pulsed injection into the yard cable. The results showed no equipment damage or undesired operations for insult levels below 180 kV peak open circuit voltage, which is significantly higher than the anticipated coupling to substation yard cables.
Resurrecting a battery chemistry thought to be only primary, we demonstrate the first example of a rechargeable alkaline zinc/copper oxide battery. With the incorporation of a Bi2O3additive to stabilize the copper oxide-based conversion cathode, Zn/(CuO-Bi2O3) cells are capable of cycling over 100 times at >124 W h/L, with capacities from 674 mA h/g (cycle 1) to 362 mA h/g (cycle 150). The crucial role of Bi2O3in facilitating the electrochemical reversibility of Cu2O, Cu(OH)2, and Cuowas supported by scanning and transmission electrochemical microscopy, cyclic voltammetry, and rotating ring-disc electrode voltammetry and monitoredvia operandoenergy-dispersive X-ray diffraction measurements. Bismuth was identified as serving two roles, decreasing the cell resistance and promoting Cu(I) and Cu(II) reduction. To mitigate the capacity losses of long-term cycling CuO cells, we demonstrate two limited depth of discharge (DOD) strategies. First, a 30% DOD (202 mA h/g) retains 99.9% capacity over 250 cycles. Second, the modification of the CuO cathode by the inclusion of additional Cu metal enables performance at very high areal capacities of ∼40 mA h/cm2and unprecedented energy densities of ∼260 W h/L, with near 100% Coulombic efficiency. This work revitalizes a historically primary battery chemistry and opens opportunity to future works in developing copper-based conversion cathode chemistries for the realization of low-cost, safe, and energy-dense secondary batteries.
Classification of features in a scene typically requires conversion of the incoming photonic field into the electronic domain. Recently, an alternative approach has emerged whereby passive structured materials can perform classification tasks by directly using free-space propagation and diffraction of light. In this manuscript, we present a theoretical and computational study of such systems and establish the basic features that govern their performance. We show that system architecture, material structure, and input light field are intertwined and need to be co-designed to maximize classification accuracy. Our simulations show that a single layer metasurface can achieve classification accuracy better than conventional linear classifiers, with an order of magnitude fewer diffractive features than previously reported. For a wavelength λ, single layer metasurfaces of size 100λ × 100λ with an aperture density λ-2 achieve ∼96% testing accuracy on the MNIST data set, for an optimized distance ∼100λ to the output plane. This is enabled by an intrinsic nonlinearity in photodetection, despite the use of linear optical metamaterials. Furthermore, we find that once the system is optimized, the number of diffractive features is the main determinant of classification performance. The slow asymptotic scaling with the number of apertures suggests a reason why such systems may benefit from multiple layer designs. Finally, we show a trade-off between the number of apertures and fabrication noise.
Despite its promise as a safe, reliable system for grid-scale electrical energy storage, traditional molten sodium (Na) battery deployment remains limited by cost-inflating high-temperature operation. Here, we describe a high-performance sodium iodide-gallium chloride (NaI-GaCl3) molten salt catholyte that enables a dramatic reduction in molten Na battery operating temperature from near 300°C to 110°C. We demonstrate stable, high-performance electrochemical cycling in a high-voltage (3.65 V) Na-NaI battery for >8 months at 110°C. Supporting this demonstration, characterization of the catholyte physical and electrochemical properties identifies critical composition, voltage, and state of charge boundaries associated with this enabling inorganic molten salt electrolyte. Symmetric and full cell testing show that the catholyte salt can support practical current densities in a low-temperature system. Collectively, these studies describe the critical catholyte properties that may lead to the realization of a new class of low-temperature molten Na batteries.
In polymer nanoparticle composites (PNCs) with attractive interactions between nanoparticles (NPs) and polymers, a bound layer of the polymer forms on the NP surface, with significant effects on the macroscopic properties of the PNCs. The adsorption and wetting behaviors of polymer solutions in the presence of a solid surface are critical to the fabrication process of PNCs. In this study, we use both classical density functional theory (cDFT) and molecular dynamics (MD) simulations to study dilute and semi-dilute solutions of short polymer chains near a solid surface. Using cDFT, we calculate the equilibrium properties of polymer solutions near a flat surface while varying the solvent quality, surface-fluid interactions, and the polymer chain lengths to investigate their effects on the polymer adsorption and wetting transitions. Using MD simulations, we simulate polymer solutions near solid surfaces with three different curvatures (a flat surface and NPs with two radii) to study the static conformation of the polymer bound layer near the surface and the dynamic chain adsorption process. We find that the bulk polymer concentration at which the wetting transition in the poor solvent system occurs is not affected by the difference in surface-fluid interactions; however, a threshold value of surface-fluid interaction is needed to observe the wetting transition. We also find that with good solvent, increasing the chain length or the difference in the surface-polymer interaction relative to the surface-solvent interaction increases the surface coverage of polymer segments and independent chains for all surface curvatures. Finally, we demonstrate that the polymer segmental adsorption times are heavily influenced only by the surface-fluid interactions, although polymers desorb more quickly from highly curved surfaces.
Fine control over the thermal expansion and contraction behavior of polymer materials is challenging. Most polymers have large coefficients of thermal expansion (CTEs), which preclude long performance lifetimes of composite materials. Herein, we report the design and synthesis of epoxy thermosets with low CTE values below their Tg and large contraction behavior above Tg by incorporating thermally contractile dibenzocyclooctane (DBCO) motifs within the thermoset network. This atypical thermomechanical behavior was rationalized in terms of a twist-boat to chair conformational equilibrium of the DBCO linkages. We anticipate these findings to be generally useful in the preparation of materials with designed CTE values.
Phosphor thermometry has been successfully applied within several challenging environments. Typically, the thermographic phosphors are excited by an ultraviolet light source, and the temperature-dependent spectral or temporal response is measured. However, this is challenging or impossible in optically thick environments. In addition, emission from other sources (e.g., a flame) may interfere with the optical phosphor emission. A temperature dependent x-ray excitation/emission could alleviate these issues as x-rays could penetrate obscurants with no interference from flame luminosity. In addition, x-ray emission could allow for thermometry within solids while simultaneously x-ray imaging the structural evolution. In this study, select thermographic phosphors were excited via x-ray radiation, and their x-ray emission characteristics were measured at various temperatures. Several of the phosphors showed varying levels of temperature dependence with the strongest sensitivity occurring for YAG:Dy and ZnGa2O4:Mn. This approach opens a path for less intrusive temperature measurements, particularly in optically opaque multiphase and solid phase combustion environments.
The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX
Villalba, Gara; Whelan, Mary; Montzka, Stephen A.; Cameron-Smith, Philip J.; Fischer, Marc; Zumkehr, Andrew; Hilton, Tim; Stinecipher, James; Baker, Ian; Bambha, Ray B.; Lafranchi, Brian W.; Estruch, Carme; Campbell, Elliott
Cities are implementing additional urban green as a means to capture CO2 and become more carbon neutral. However, cities are complex systems where anthropogenic and natural components of the CO2 budget interact with each other, and the ability to measure the efficacy of such measures is still not properly addressed. There is still a high degree of uncertainty in determining the contribution of the vegetation signal, which furthermore confounds the use of CO2 mole fraction measurements for inferring anthropogenic emissions of CO2. Carbonyl sulfide (OCS) is a tracer of photosynthesis which can aid in constraining the biosphere signal. This study explores the potential of using OCS to track the urban biosphere signal. We used the Sulfur Transport and dEposition Model (STEM) to simulate the OCS concentrations and the Carnegie Ames Stanford Approach ecosystem model to simulate global CO2 fluxes over the Bay Area of San Francisco during March 2015. Two observation towers provided measurements of OCS and CO2: The Sutro tower in San Francisco (upwind from the area of study providing background observations), and a tower located at Sandia National Laboratories in Livermore (downwind of the highly urbanized San Francisco region). Our results show that the STEM model works better under stable marine influence, and that the boundary layer height and entrainment are driving the diurnal changes in OCS and CO2 at the downwind Sandia site. However, the STEM model needs to better represent the transport and boundary layer variability, and improved estimates of gross primary productivity for characterizing the urban biosphere signal are needed.
Potts, Alexander M.; Bajaj, Sanyam; Daughton, David R.; Allerman, A.A.; Armstrong, Andrew A.; Razzak, Towhidur; Sohel, Shahadat H.; Rajan, Siddharth
In this work, ultrawide bandgap Al0.7Ga0.3N MESFETs with refractory Tungsten Schottky and Ohmic contacts are studied in 300–675 K environments. Variable-temperature dc electrical transport reveals large ON-state drain current densities for an AlGaN device: 209 mA/mm at 300 K and 156 mA/mm at 675 K in the ON-state (25% reduction). Drain and gate currents are only weakly temperature-dependent, suggesting potential for engineering temperature invariant operation. The ON-/ OFF-ratio is limited by OFF-state leakage through the gate, which is attributed to damage from sputter deposition. Future work using refractory metals with larger work functions that are deposited by electron beam deposition is proposed.
The addition of a common amino acid, phenylalanine, to a Layer-by-Layer (LbL) deposited polyelectrolyte (PE) film on a nanoporous membrane can increase its ionic selectivity over a PE film without the added amino acid. The addition of phenylalanine is inspired by detailed knowledge of the structure of the channelrhodopsins family of protein ion channels, where phenylalanine plays an instrumental role in facilitating sodium ion transport. The normally deposited and crosslinked PE films increase the cationic selectivity of a support membrane in a controllable manner where higher selectivity is achieved with thicker PE coatings, which in turn also increases the ionic resistance of the membrane. The increased ionic selectivity is desired while the increased resistance is not. We show that through incorporation of phenylalanine during the LbL deposition process, in solutions of NaCl with concentrations ranging from 0.1 to 100 mM, the ionic selectivity can be increased independently of the membrane resistance. Specifically, the addition is shown to increase the cationic transference of the PE films from 81.4% to 86.4%, an increase on par with PE films that are nearly triple the thickness while exhibiting much lower resistance compared to the thicker coatings, where the phenylalanine incorporated PE films display an area specific resistance of 1.81 Ω cm2in 100 mM NaCl while much thicker PE membranes show a higher resistance of 2.75 Ω cm2in the same 100 mM NaCl solution.