Publications

Most multivalent secondary batteries have employed electrolytes composed of cyclic ether solvents such as tetrahydrofuran or linear glycol ether solvents (glymes) such as 1,2-dimethoxyethane (G1). A robust understanding of multivalent cation solvation tendencies in these classes of solvents provides insight into corresponding structure-property relationships which, in turn, promotes the design and discovery of improved electrolytes. In this work, our goal is to systematically address how electrolyte constituent properties, namely, ether solvent structure and dication size, direct the solvation interactions of divalent electrolytes and their resultant properties. This study utilizes pulsed-field gradient (PFG) nuclear magnetic resonance (NMR) spectroscopy in conjunction with Raman spectroscopy and ionic conductivity measurements to elucidate the preferential interactions between multivalent cations, anions, and solvent molecules along with their correlated ion dynamics. These investigations incorporate two representative divalent cations (Ca2+ and Zn2+) as well as two ethereal solvent representatives from both the cyclic ether and glyme structural classes. The results reveal that anions coordinate more readily with divalent cations in cyclic ethers than in glymes. Furthermore, the coordination of the anions with Ca2+, i.e., contact-ion pair (CIP) formation is more pronounced than with Zn2+ in a glyme solvent of limited chain length (G1), providing insight into cation size effects that are important for translating solvation behavior across various multivalent electrolytes. Importantly, we find that specific anion coordination is more strongly controlled by solvent structure than by salt concentration in the practical range of 0.1-0.5 M. However, simply reducing these inner-sphere inter-ionic interactions by changing solvent structure does not necessarily de-correlate ionic motion. Instead, concentration-dependent changes in molar ionic conductivity suggest that second-shell interactions, i.e., solvent separated ion pairs (SSIPs), are prevalent in these electrolytes and that the solution dielectric constant, which is increased by the presence of dipolar ion pairs, is critical for controlling these interactions. These findings thus provide a basis for understanding the physical chemistry of multivalent battery electrolytes.

Two events of magnitude (m_b) 3.6–3.8 occurred in southern North Korea (NK) on 27 June 2019 and 11 May 2020. Although these events were located ~330–400 km from the known nuclear test site, the fact that they occurred within the territory of NK, a country with a recent history of underground nuclear tests, made them events of interest for the monitoring community. In this work, we used P/Lg ratios from regional stations to categorize seismic events that occurred in NK from 2006 to May 2020, including these two recent events, the six declared NK nuclear tests, and the cavity collapse and triggered earthquakes that followed the 3 September 2017 nuclear explosion. We were able to separate the cavity collapse from the population of nuclear explosions. However, based on P/Lg ratios, the distinction between the earthquakes and the cavity collapse is ambiguous. The performed discriminant analyses suggest that combining Pg/Lg and Pn/Lg ratios results in improved discriminant power compared with any of the ratio types alone. We used the two ratio types jointly in a quadratic discriminant function and successfully classified the six declared nuclear tests and the triggered earthquakes that followed the September 2017 explosion. Our analyses also confirm that the recent southern events of June 2019 and May 2020 are both tectonic earthquakes that occurred naturally.

We present a methodology based on the Néel model to build a classical spin-lattice Hamiltonian for cubic crystals capable of describing magnetic properties induced by the spin-orbit coupling like magnetocrystalline anisotropy and anisotropic magnetostriction, as well as exchange magnetostriction. Taking advantage of the analytical solutions of the Néel model, we derive theoretical expressions for the parametrization of the exchange integrals and Néel dipole and quadrupole terms that link them to the magnetic properties of the material. This approach allows us to build accurate spin-lattice models with the desired magnetoelastic properties. We also explore a possible way to model the volume dependence of magnetic moment based on the Landau energy. This feature allows us to consider the effects of hydrostatic pressure on the saturation magnetization. We apply this method to develop a spin-lattice model for BCC Fe and FCC Ni, and we show that it accurately reproduces the experimental elastic tensor, magnetocrystalline anisotropy under pressure, anisotropic magnetostrictive coefficients, volume magnetostriction, and saturation magnetization under pressure at zero temperature. This work could constitute a step towards large-scale modeling of magnetoelastic phenomena.

Polymer foams are an essential class of lightweight materials used to protect assets against mechanical insults, such as shock and vibration. Two features are important to enhance their energy absorption characteristics: the foam structure and the matrix phase mechanical behavior. This study investigates novel approaches to control both of these features to enhance the energy absorption capability of flexible lattice foams. First, we consider 3D printing via digital light processing (DLP) as a method to control the foam mesostructure across a suite of periodic unit cells. Second, we introduce an additional energy dissipation mechanism in the solid matrix phase material by 3D printing the lattice foams with polydomain liquid crystal elastomer (LCE), which undergo a mechanically induced phase transition under large strains. This phase transition is associated with LC mesogen rotation and alignment and provides a second mechanism for mechanical energy dissipation in addition to the viscoelastic relaxation of the polymer network. We contrast the 3D printed LCE lattices with conventional, thermomechanically near-equivalent elastomer lattice foams to quantify the energy-absorbing enhancement the LCE matrix phase provides. Under cyclic quasi-static uniaxial compression conditions, the LCE lattices show dramatically enhanced energy dissipation in uniaxial compression compared to the non-LCE equivalent foams printed with a commercially available photocurable elastomer resin. The lattice geometry also plays a prominent role in determining the energy dissipation ratio between the LCE and non-LCE foams. We show that when increasing the lattice connectivity, the foam deformation transitions from bending-dominated to stretching-dominated deformations, which generates higher axial strains in the struts and higher energy dissipation in the lattice foam, as stretching allows greater mesogen rotation than bending. The LCE foams demonstrate superior energy absorption during the repeated dynamic loading during drop testing compared with the non-LCE equivalent foams, demonstrating the potential of LCEs to enhance physical protection systems against mechanical impact.

Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.

Ionizable block copolymers with distinctive block characteristics display the diversity crucial for the design of macromolecules for targeted applications. In contrast to van der Waals copolymers, their interfaces, which are critical to their function, consist of nanodomains, each of a different nature and thus unique interfacial behavior. Here, the interfacial response of a symmetric block copolymer with a sulfonated polystyrene polyelectrolyte center, tethered to polyethylene-r-propylene and terminated by poly(t-butyl styrene) is probed as polymer films are exposed to three polar solvents, water, propanol, and tetrahydrofuran (THF), using molecular dynamics simulations. Each of the solvents captures a distinctive interaction with the individual blocks. We find that at the film boundary, the interfacial response is initially dominated by that of the hydrophobic blocks to all solvents. At later times, the solvent distribution among the blocks, where water molecules associate predominantly with the sulfonated groups and propanol and THF reside at multiple different sites, determines the chemical composition and the polymer conformation at the interface. Overall, these simulations provide the first direct molecular insight into the interfacial response of ionizable copolymers.

Steady-state shear viscosity (γ˙) of unconcatenated ring polymer melts as a function of the shear rate γ˙ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity ∼γ˙-0.56 ± 0.02, independent of chain length, for Weissenberg numbers up to about 102. Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with ∼γ˙-0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high γ˙. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ∼104, the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime ∼γ˙-0.57, we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large γ˙. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/γ˙. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox-Merz rule for ring polymers.

Electrochemical systems function via interconversion of electric charge and chemical species and represent promising technologies for our cleaner, more sustainable future. However, their development time is fundamentally limited by our ability to identify new materials and understand their electrochemical response. To shorten this time frame, we need to switch from the trial-and-error approach of finding useful materials to a more selective process by leveraging model predictions. Machine learning (ML) offers data-driven predictions and can be helpful. Herein we ask if ML can revolutionize the development cycle from decades to a few years. We outline the necessary characteristics of such ML implementations. Instead of enumerating various ML algorithms, we discuss scientific questions about the electrochemical systems to which ML can contribute.

We present the design, fabrication, and initial characterization of a CMOS compatible, ultra-high bandwidth, bulk-micro machined, optomechanical accelerometer. Displacement detection is achieved via a SiN integrated photonics Mach-Zehnder interferometer (MZI) fabricated on the surface of the device that is optomechanically coupled to acceleration-induced deformation of the accelerometer's proof mass tethers. The device is designed to measure vibrations at microsecond timescales with high dynamic range for the characterization of shock dynamics.

While there is no known fundamental lower limit to the thermal conductivity of a material, the lowest thermal conductivities are typically found in amorphous and strongly disordered materials, not highly crystalline materials. Here, we demonstrate a surprising nanostructuring route to ultralow thermal conductivity in a large-unit-cell oxide crystal (Fe₃O₄) containing close-packed nanoscale pores. The electrical conductivity of this material reduces by a factor of 5 relative to dense v, independent of pore size. In contrast, thermal conductivity has a strong dependence on pore size with a factor of 40 of suppression relative to dense Fe₃O₄ for 40 nm pores vs a factor of 5 for 500 nm pores. The matrix thermal conductivity of Fe₃O₄ containing 40 nm pores falls below the predicted minimum thermal conductivity by a factor of 3. Finally, we attribute this to strong acoustic phonon scattering and intrinsically limited contributions to thermal conductivity from optical phonons with small dispersion.

The flash pyrolysis oftrans3-pentenenitrile (3-PN, CH3-CH-CH-CH2-CN) was studied by combining the results of VUV photoionization mass spectra with broadband microwave spectra recorded as a function of the temperature of the pyrolysis tube. The two separated functional groups (vinyl and nitrile) open up isomerization as an initial step in competition with unimolecular dissociation. Primary products were detected by keeping the 3-PN concentration low and limiting reaction times to the traversal time of the gas in the pyrolysis tube (100 μs). The reaction is quenched and products are cooled by expansion into vacuum before interrogation over the 8-18 GHz region using chirped-pulse broadband methods. 118 nm VUV photoionization of the same reaction mixture provides a means of detecting all products with ionization potentials below 10.5 eV with minimal fragmentation. These results are combined with a detailed computational investigation of the C5H7N and related potential energy surfaces, leading to a consistent picture of the unimolecular decomposition of 3-PN. Loss of two H-atoms to form a 79 amu product is proven from its microwave transitions to containtrans-Z-2,4-pentadienenitrile, while no pyridine is observed. Methyl loss, HCN loss, and breaking the central C(2)-C(3) bond all occur following isomerization of the position of the double bond, thereby opening up low-energy pathways to these decomposition channels.

We demonstrate all-epitaxial structures capable of supporting short- and long-range surface plasmon polariton (SRSPP and LRSPP) modes in the long-wave infrared region of the electromagnetic spectrum. The SRSPP and LRSPP modes are bound to the interfaces of a buried heavily doped (n + +) semiconductor layer and surrounding quantum-engineered type-II superlattice (T2SL) materials. The surrounding T2SLs are designed to allow optical transitions across the frequency dispersion of the SPP modes. We map the SPP dispersion in our structure using grating-coupled angle- and polarization-dependent reflection and photoluminescence spectroscopy. The epitaxial structures are analytically described using a simplified three-layer system (T 2 SL / n + + / T 2 SL) and modeled using rigorous coupled wave analysis with excellent agreement to our experimental results. The presented structures offer the potential to serve as long-range interconnects or waveguides in all-epitaxial plasmonic/optoelectronic systems operating in the long-wave infrared.

This report describes plans for dust sampling and analysis for the multi-year Canister Deposition Field Demonstration. The demonstration will use three commercial 32PTH2 NUHOMS welded stainless steel storage canisters, which will be stored at an ISFSI site in Advanced Horizontal Storage Modules. One canister will be unheated; the other two will have heaters to achieve canister surface temperatures that match, to the degree possible, spent nuclear fuel (SNF) loaded canisters with heat loads of 10 kW and 40 kW. Surface sampling campaigns will take place on a yearly or bi-yearly basis. The goal of the planned dust sampling and analysis is to determine important environmental parameters that impact the potential occurrence of stress corrosion cracking on SNF dry storage canisters. Specifically, the size, morphology, and composition of the deposited dust and salt particles will be quantified, as well as the soluble salt load per unit area and the rate of deposition, as a function of canister surface temperature, location, time, and orientation. Sampling locations on the canister surface will nominally include 25 locations, corresponding to 5 circumferential locations at each of the 5 longitudinal locations. At each sampling location, a 2x2 sampling grid (containing 4 sample cells) will be painted onto the metal surface. During each sampling campaign, two samples at each sampling location will be collected, in a specific routine to measure both periodic (yearly or bi-yearly) and cumulative deposition rates. For each sample, a wet and a dry sample will be collected. Wet samples will be analyzed to determine the composition of the soluble salt fraction and to estimate salt loading per unit area. Dry samples will be analyzed to assess particle size, morphology, mineralogy, and identity (e.g. for floral/faunal fragments). The data generated by this proposed sampling plan will provide detailed information on dust and salt aerosol deposits on spent nuclear fuel canister surfaces. The anticipated results include information regarding particle compositions, size distributions, and morphologies, in addition to particle deposition rates as a function of canister surface location, orientation, time, and temperature. The information gathered during the Canister Deposition Field Demonstration is critical for ongoing efforts to develop a detailed understanding of the potential for stress corrosion cracking on SNF dry storage canisters

High and low salinity water flooding are common oil recovery processes performed in the oil fields for extracting crude oil from the reservoir. These processes are often performed sequentially, naturally establishing non-uniform salinity in the porous subsurface. In this article, we investigate oil transport in porous media induced by salinity change upon flooding with high and low salinity water. As we observe a large number of impervious dead-ends from three-dimensional imaging of the actual reservoir, we identify that these areas play an important role in oil recovery where the oil transport is governed by the salinity change rather than hydrodynamics. The salinity gradients induced upon high salinity water flooding provide pathways to enhance the transport of oil drops trapped in the dead-end regions via non-equilibrium effects. However, above a critical salinity, we observe a rapid aggregation of drops that lead to the complete blockage of the pore space, thereby inhibiting oil recovery. We also find that, at an intermediate salinity where the drop aggregation is modest, the aggregation rather promotes the oil recovery. Our observations suggest that there exist optimal salinity conditions for maximizing oil recovery during chemical flooding.

The fillers research and development (R&D) program, mostly experimental, is part of a broader R&D program that includes new process modeling and performance assessment of criticality effects and the overall importance of criticality to repository performance (consequence screening). A literature research and consultation effort with experts by Hardin and Brady (2018) identified several potentially effective and workable filler materials including cements (primarily phosphate based), moltenmetal alloys, and low-temperature glasses. Filler attributes were defined, and the preliminary lists were compared qualitatively. Further comparative analysis will be done (e.g., cost estimates) after experimental screening has narrowed the list of alternatives. The following cement filler compositions were selected for experimental development work and accelerated testing in FY20: Aluminum phosphate cements (APCs); more specifically aluminum oxide / aluminum phosphate (Al₂O₃ / AlPO₄) cements in which Al₂O₃ serves as the filler material bound by an AlPO₄ binder formed by the reaction of Al₂O₃ with H₃PO₄; Calcium phosphate cements (CPCs); more specifically composed of pure or nearly pure hydroxyapatite or HAP (Ca₅(PO₄)₃(OH)); Wollastonite phosphate cements (WPC), specifically wollastonite and aluminum or calcium aluminum phosphates in which CaSiO₃ serves as the filler material and the phosphate serves as the binder. The FY20 effort focused on the optimization of compositions and subsequent processing of these three materials to achieve dense and well-consolidated monolithic samples with 30 to 40% porosity and permeabilities of 1 millidarcy. At the close of this progress report the aluminum phosphate cements (APCs) and the wollastonite phosphate cements (WPCs) appear to show the most promise for continued development. Less progress has been made with the calcium phosphate cements (CPCs); their slurry viscosities are high (and difficult to measure) and they exhibit relatively short cure times of 2 to 3 hours with concomitant and excessive volatile (e.g. CO₂) generation.

We experimentally show that the thermal conductance across confined solid-solution crystalline thin films between parent materials does not necessarily lead to an increase in thermal resistances across the thin-film geometries with increasing film thicknesses, which is counterintuitive to the notion that adding a material serves to increase the total thermal resistance. Confined thin epitaxial Ca0.5Sr0.5TiO3 solid-solution films with systematically varying thicknesses in between two parent perovskite materials of calcium titanate and (001)-oriented strontium titanate are grown, and thermoreflectance techniques are used to accurately measure the thermal boundary conductance across the confined solid-solution films, showing that the thermal resistance does not substantially increase with the addition of solid-solution films with increasing thicknesses from μ1 to μ10 nm. Contrary to the macroscopic understanding of thermal transport where adding more material along the heat propagation direction leads to larger thermal resistances, our results potentially offer experimental support to the computationally predicted concept of vibrational matching across interfaces. This concept is based on the fact that a better match in the available heat-carrying vibrations due to an interfacial layer can lead to lower thermal boundary resistances, thus leading to an enhancement in thermal boundary conductance across interfaces driven by the addition of a thin "vibrational bridge"layer between two solids.

The fillers research and development (R&D) program, mostly experimental, is part of a broader R&D program that includes new process modeling and performance assessment of criticality effects and the overall importance of criticality to repository performance (consequence screening). A literature research and consultation effort with experts by Hardin and Brady (2018) identified several potentially effective and workable filler materials including cements (primarily phosphate based), moltenmetal alloys, and low-temperature glasses. Filler attributes were defined, and the preliminary lists were compared qualitatively. Further comparative analysis will be done (e.g., cost estimates) after experimental screening has narrowed the list of alternatives. The research presented here is focused Sandia’s efforts for the development of phosphate-based cement fillers. Molten metal filler research is an ongoing activity at Oak Ridge National Laboratories and is not discussed herein. After the completion of the FY20 research effort the following cement filler compositions were selected for further experimental development work and advanced testing in FY21: 1. Aluminum phosphate cements (APCs); more specifically aluminum oxide / aluminum phosphate (Al₂O₃ / AlPO₄) cements in which Al₂O₃ serves as the filler material bound by an AlPO₄ binder formed by the reaction of Al₂O₃ with various phosphate sources; 2. Wollastonite phosphate cements (WPCs), specifically wollastonite and aluminum or calcium aluminum phosphates in which CaSiO₃ serves as the filler material bound by a calcium phosphate that serves as the binder; and 3. Calcium aluminate phosphate cements (CAPCs) specifically grossite (CaAl₄O₇) and hibonite (CaAl₁₁O₁₈) fillers bound by an aluminum phosphate that serves as the binder. This effort focused on the optimization and subsequent processing of these three cements to achieve dense and well-consolidated monolithic samples. Upon completion of the FY21 effort the aluminum phosphate cements (APCs) and the calcium aluminate phosphate cements (CAPCs) show the most promise for advanced testing and scale up. We will begin the work in FY22 focused on testing the performance of these two cements in small scale DPCs as well as advanced materials testing to evaluate cement performance under expected radiation doses and representative post-closure geochemical environments.

The natural convection boundary layer (${\delta }_{nc}$) and its influence on cathodic current in a galvanic couple under varying electrolytes as a function of concentration (1 - 5.3 M NaCl) and temperature (25 °C-45 °C) were understood. Polarization scans were obtained under quiescent conditions and at defined boundary layer thicknesses using a rotating disk electrode on platinum and stainless steel 304L (SS304L); these were combined to determine ${\delta }_{nc}.$ With increasing chloride concentration and temperature, ${\delta }_{nc}$ decreased. Increased mass transport (Sherwood number) results in a decrease in ${\delta }_{nc},$ providing a means to predict this important boundary. Using Finite Element Modeling, the cathodic current was calculated for an aluminum alloy/SS304L galvanic couple as a function of water layer (WL) thickness and cathode length. Electrolyte domains were delineated, describing (i) dominance of ohmic resistance over mass transport under thin WL, (ii) the transition from thin film to bulk conditions at ${\delta }_{nc},$ and (iii) dominance of mass transport under thick WL. With increasing chloride concentration, cathodic current decreased due to decreases in mass transport. With increasing temperature, increased cathodic current was related to increases in mass transport and solution conductivity. This study has implications for sample sizing and corrosion prediction under changing environments.

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.0 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and $\textit{J}$-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

This is an addendum to the Sierra/SolidMechanics 5.0 User’s Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State’s International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 5.0 User’s Guide should be referenced for most general descriptions of code capability and use.

Presented in this document are the theoretical aspects of capabilities contained in the Sierra/SM code. This manuscript serves as an ideal starting point for understanding the theoretical foundations of the code. For a comprehensive study of these capabilities, the reader is encouraged to explore the many references to scientific articles and textbooks contained in this manual. It is important to point out that some capabilities are still in development and may not be presented in this document. Further updates to this manuscript will be made as these capabilities come closer to production level.

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra / SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra / SM test suite that are not included in this manual.

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra / SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra / SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 5.0 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

Thermally activated batteries undergo a series of coupled physical changes during activation that influence battery performance. These processes include energetic material burning, heat transfer, electrolyte phase change, capillary-driven two-phase porous flow, ion transport, electrochemical reactions, and electrical transport. Several of these processes are strongly coupled and have a significant effect on battery performance, but others have minimal impact or may be suitably represented by reduced-order models. Assessing the relative importance of these phenomena must be based on comparisons to a high-fidelity model including all known processes. In this work, we first present and demonstrate a high-fidelity, multi-physics model of electrochemical performance. This novel multi-physics model enables predictions of how competing physical processes affect battery performance and provides unique insights into the difficult-to-measure processes that happen during battery activation. We introduce four categories of model fidelity that include different physical simplifications, assumptions, and reduced-order models to decouple or remove costly elements of the simulation. Using this approach, we show an order-of-magnitude reduction in computational cost while preserving all design-relevant quantities of interest within 5 percent. The validity of this approach and these model reductions is demonstrated by comparison between results from the full fidelity model and the different reduced models.

Lanthanide elements have well-documented similarities in their chemical behavior, which make the valuable trivalent lanthanide cations (Ln3+) particularly difficult to separate from each other in water. In this work, we applyab initiomolecular dynamics simulations to compare the free energies (ΔGads) associated with the adsorption of lanthanide cations to silica surfaces at a pH condition where SiO−groups are present. The predicted ΔGadsfor lutetium (Lu3+) and europium (Eu3+) are similar within statistical uncertainties; this is in qualitative agreement with our batch adsorption measurements on silica. This finding is remarkable because the two cations exhibit hydration free energies (ΔGhyd) that differ by >2 eV, different hydration numbers, and different hydrolysis behavior far from silica surfaces. We observe that the similarity in Lu3+and Eu3+ΔGadsis the result of a delicate cancellation between the difference in Eu3+and Lu3+hydration (ΔGhyd), and their difference in binding energies to silica. We propose that disrupting this cancellation at the two end points, either for adsorbed or completely desorbed lanthanides (e.g.,viananoconfinment or mixed solvents), will lead to effective Ln3+separation.

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approachis a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D2O confined inhydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on thesilica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. Thisappears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated withtheweaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.

The detection of special nuclear materials (SNM) requires the understanding of nuclear signatures that allow the discrimination against background. In particular, understanding neutron background characteristics such as count rates and energies and their correlations with environmental conditions and surroundings of measurement locations is important in enhancing SNM detection capabilities. The Mobile Imager of Neutrons for Emergency Responders (MINER) was deployed for 8 weeks in downtown San Francisco (CA) to study such neutron background characteristics in an urban environment. Of specific interest was the investigation of the impact of surrounding buildings on the neutron background count rates and to answer the question whether buildings act as absorber of neutrons or as sources via the so-called ship effect. MINER consists of 16 liquid scintillator detector elements and can be operated as a neutron spectrometer, as a neutron imager, or simply as a counter of fast neutrons. As expected, the neutron background rate was found to be inversely proportional to the atmospheric pressure. In the energy range where MINER is most sensitive, approximately 1–10 MeV, it was found that the shape of the detected background spectrum is similar to that of a detected fission spectrum, indicating the limited discrimination power of the neutron energy. The similarities between the detected background neutron spectrum and fission sources makes it difficult to discriminate SNM from background based solely on the energies observed. The images produced using maximum likelihood expectation maximization revealed that neutrons preferentially are coming from areas in the environment that have open sky, indicating that the surrounding buildings act as absorbers of neutrons rather than sources as expected by the ship effect. The inherent properties of a neutron scatter camera limit the achievable image quality and the effective deployment to systematically map neutron background signatures due to the low count rate.

While magnetized turbulence is ubiquitous in many astrophysical and terrestrial systems, our understanding of even the simplest physical description of this phenomena, ideal magnetohydrodynamic (MHD) turbulence, remains substantially incomplete. As such, in this work we highlight the shortcomings of existing theoretical and phenomenological descriptions of MHD turbulence that focus on the joint (kinetic and magnetic) energy fluxes and spectra by demonstrating that treating these quantities separately enables fundamental insights into the dynamics of MHD turbulence. This is accomplished through the analysis of the scale-wise energy transfer over time within an implicit large eddy simulation of subsonic, super-Alfvénic MHD turbulence. Our key finding is that the kinetic energy spectrum develops a scaling of approximately k^–4/3 in the stationary regime as magnetic tension mediates large-scale kinetic to magnetic energy conversion and significantly suppresses the kinetic energy cascade. This motivates a reevaluation of existing MHD turbulence theories with respect to a more differentiated modeling of the energy fluxes.

The III-nitride semiconductors have many attractive properties for field-emission vacuum electronics, including high thermal and chemical stability, low electron affinity, and high breakdown fields. Here, we report top-down fabricated gallium nitride (GaN)-based nanoscale vacuum electron diodes operable in air, with record ultralow turn-on voltages down to ∼0.24 V and stable high field-emission currents, tested up to several microamps for single-emitter devices. We leverage a scalable, top-down GaN nanofabrication method leading to damage-free and smooth surfaces. Gap-dependent and pressure-dependent studies provide new insights into the design of future, integrated nanogap vacuum electron devices. The results show promise for a new class of high-performance and robust, on-chip, III-nitride-based vacuum nanoelectronics operable in air or reduced vacuum.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User’s Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user’s notes and of course the material in the open literature.

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

Graphite fluoride (CFx) cathodes coupled with lithium anodes yield one of the highest theoretical specific capacities (>860 mAh/g) among primary batteries. In practice, the observed discharge voltage (∼2.5 V) is significantly lower than thermodynamic limits (>4.5 V), the discharge rate is low, and so far Li/CFx has only been used in primary batteries. Understanding the discharge mechanism at atomic length scales will improve practical CFx energy density, rate capability, and rechargeability. So far, purely experimental techniques have not identified the correct discharge mechanism or explained the discharge voltage. We apply density functional theory calculations to demonstrate that a CFx-edge propagation discharge mechanism based on lithium insertion at the CF/C boundary in partially discharged CFx exhibits a voltage range of 2.5 to 2.9 V - depending on whether solvent molecules are involved. The voltages and solvent dependence agree with our discharge and galvanostatic intermittent titration technique measurements. The predicted discharge kinetics are consistent with CFx operations. Finally, we predict some Li/CFx rechargeability under the application of high potentials, along a charging pathway different from that of discharge. Our work represents a general, quasi-kinetic framework to understand the discharge of conversion cathodes, circumventing the widely used phase diagram approach which most likely does not apply to Li/CFx because equilibrium conditions are not attained in this system.

The theoretical understanding of plasmon behavior is crucial for an accurate interpretation of inelastic scattering diagnostics in many experiments. We highlight the utility of linear response time-dependent density functional theory (LR-TDDFT) as a first-principles framework for consistently modeling plasmon properties. We provide a comprehensive analysis of plasmons in aluminum from ambient to warm dense matter conditions and assess typical properties such as the dynamical structure factor, the plasmon dispersion, and the plasmon lifetime. We compare our results with scattering measurements and with other TDDFT results as well as models such as the random phase approximation, the Mermin approach, and the dielectric function obtained using static local field corrections of the uniform electron gas parametrized from path-integral Monte Carlo simulations. We conclude that results for the plasmon dispersion and lifetime are inconsistent between experiment and theories and that the common practice of extracting and studying plasmon dispersion relations is an insufficient procedure to capture the complicated physics contained in the dynamic structure factor in its full breadth.

The elastic moduli of amorphous and crystalline atomic layer-deposited Hf1-xZrxO2 (HZO, x = 0, 0.31, 0.46, 0.79, 1) films prepared with TaN electrodes on silicon substrates were investigated using picosecond acoustic measurements. The moduli of the amorphous films were observed to increase between 211 ± 6 GPa for pure HfO2 and 302 ± 9 GPa for pure ZrO2. In the crystalline films, it was found that the moduli increased upon increasing the zirconium composition from 248 ± 6 GPa for monoclinic HfO2 to 267 ± 9 GPa for tetragonal ZrO2. Positive deviations from this increase were observed for the Hf0.69Zr0.31O2 and Hf0.54Zr0.46O2 compositions, which were measured to have moduli of 264 ± 8 GPa and 274 ± 8 GPa, respectively. These two compositions contained the largest fractions of the ferroelectric orthorhombic phase, as assessed from polarization and diffraction data. The biaxial stress states of the crystalline films were characterized through sin2(ψ) x-ray diffraction analysis. The in-plane stresses were all found to be tensile and observed to increase with the increasing zirconium composition, between 2.54 ± 0.6 GPa for pure HfO2 and 5.22 ± 0.5 GPa for pure ZrO2. The stresses are consistent with large thermal expansion mismatches between the HZO films and silicon substrates. These results demonstrate a device-scale means to quantify biaxial stress for investigation on its effect on the ferroelectric properties of hafnia-based materials.

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r₀. A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r₀. As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

Invited for this month's cover is the joint redox flow battery team from Sandia and Los Alamos National Laboratories. The cover image shows the stylized components of a redox flow battery (RFB) in the foreground, with renewable sources of energy generation in the background. The Review itself is available at 10.1002/cssc.202002354.

This article presents a global prioritization methodology that evaluates the relative risks of non-state actor acquisition of materials that could be used in chemical, biological, radiological, nuclear and high explosive Weapons of Mass Destruction (WMD) from the country’s relevant infrastructure. Prioritization is based on three domains: 1. Assessing relative scale of materials in each country, 2. The country’s corresponding security posture, and 3. The presence of threat actors. The output is a list of countries prioritized from greatest risk to least. Rather than providing an overall 1 to N ranking, however, the results are placed into tiers based upon their natural groupings within the three domains. The countries in the highest tiers are flagged as potential US national security concern; those scoring in the middle and at the bottom are flagged as posing lower US national security concern. A systematic approach assesses each country by leveraging many disciplines, such as risk and decision analysis, as well as expert judgement. A quantitative value model based on Multi-Attribute Value Theory (MAVT) organizes the objectives scoring criteria into a value tree using lessons learned from previous studies, published literature, and expert judgement. The article presents the prioritization categories and corresponding value model scoring criteria to include measurement type, weight, range, and value preference. Country names and data are notional in order to share the details on the underlying methodology and model without identification of actual security risks. A deliberative process addresses factors external to the model and scrutinizes inputs, methodology, model, and results.

An increase in Electric Vehicles (EV) will result in higher demands on the distribution electric power systems (EPS) which may result in thermal line overloading and low voltage violations. To understand the impact, this work simulates two EV charging scenarios (home-and work-dominant) under potential 2030 EV adoption levels on 10 actual distribution feeders that support residential, commercial, and industrial loads. The simulations include actual driving patterns of existing (non-EV) vehicles taken from global positioning system (GPS) data. The GPS driving behaviors, which explain the spatial and temporal EV charging demands, provide information on each vehicles travel distance, dwell locations, and dwell durations. Then, the EPS simulations incorporate the EV charging demands to calculate the power flow across the feeder. Simulation results show that voltage impacts are modest (less than 0.01 p.u.), likely due to robust feeder designs and the models only represent the high-voltage (“primary”) system components. Line loading impacts are more noticeable, with a maximum increase of about 15%. Additionally, the feeder peak load times experience a slight shift for residential and mixed feeders (≈1 h), not at all for the industrial, and 8 h for the commercial feeder.

Laser powder bed fusion (LPBF) additive manufacturing (AM) offers a variety of advantages over traditional manufacturing, however its usefulness for manufacturing of high-performance components is currently hampered by internal defects (porosity) created during the LPBF process that have an unknown impact on global mechanical performance. By inducing porosity distributions through variations in print energy density and inspecting the resulting tensile samples using computed tomography, nearly 50,000 pores across 75 samples were identified. Porosity characteristics were quantitatively extracted from inspection data and compared with mechanical properties to understand the strength of relationships between porosity and global tensile performance. Useful porosity characteristics were identified for prediction of part performance. Results indicate that ductility and strain at ultimate tensile strength are the global tensile properties most significantly impacted by porosity and can be predicted with reasonable accuracy using simple porosity shape descriptors such as volume, diameter, and surface area. Moreover, it was found that the largest pores influenced behavior most significantly. Specifically, pores in excess of 125 µm in diameter were found to be a sufficient threshold for property estimation. These results establish an initial understanding of the complex defect-performance relationship in AM 316L stainless steel and can be leveraged to develop certification standards and improve confidence in part quality and reliability for the broader set of engineering alloys.

Joint and marginal distributions in the frequency, direction, and time domain are employed to demonstrate their value for wave energy resource characterization when full spectra are available. Insights gained through analysis of these distributions support wave energy converter concept design, operation and maintenance. Spatial trends in the wave energy resource and contributing wave energy systems along the continental shelf of the West Coast of the United States are investigated using the most recent two-dimensional wave spectra measurements at four buoys over an eleven year period (2008 to 2018). Resource hot spots and dominant resolved energy resource bands in the frequency-direction-time domain are delineated. Resource attributes, including frequency and directional spreading, and seasonal variability, are characterized using joint distributions and marginal distributions of wave power spectra. North Pacific westerly swells in the winter season, augmented by Aleutian low-pressure southwesterly swells, are the principal suppliers of the dominant resource and main drivers influencing resource attributes. The modification of these systems southward, especially the North Pacific westerly swells, explains the observed spatial resource trends. The dominant resource wave period shifts two seconds to higher wave periods, thirty degrees in the dominant direction band to a more northward orientation, and forward by one month.

The inductively driven transmission line (IDTL) is a miniature current-carrying device that passively couples to fringe magnetic fields in the final power feed on the Z Pulsed Power Facility. The IDTL redirects a small amount of Z's magnetic energy along a secondary path to ground, thereby enabling pulsed power diagnostics to be driven in parallel with the primary load for the first time. IDTL experiments and modeling presented here indicate that IDTLs operate non-perturbatively on Z and that they can draw in excess of 150 kA of secondary current, which is enough to drive an X-pinch backlighter. Additional experiments show that IDTLs are also capable of making cleaner, higher-fidelity measurements of the current flowing in the final feed.

We present a fully discrete approximation technique for the compressible Navier–Stokes equations that is second-order accurate in time and space, semi-implicit, and guaranteed to be invariant domain preserving. The restriction on the time step is the standard hyperbolic CFL condition, i.e. τ≲O(h)∕V where V is some reference velocity scale and h the typical meshsize.

We have commissioned a new time-resolved, x-ray imaging diagnostic for the Z facility. The primary intended application is for diagnosing the stagnation behavior of Magnetized Liner Inertial Fusion (MagLIF) and similar targets. We have a variety of imaging systems at Z, both time-integrated and time-resolved, that provide valuable x-ray imaging information, but no system at Z up to this time provides a combined high-resolution imaging with multi-frame time resolution; this new diagnostic, called TRICXI for Time Resolved In-Chamber X-ray Imager, is meant to provide time-resolved spatial imaging with high resolution. The multi-frame camera consists of a microchannel plate camera. A key component to achieving the design goals is to place the instrument inside the Z vacuum chamber within 2 m of the load, which necessitates a considerable amount of x-ray shielding as well as a specially designed, independent vacuum system. A demonstration of the imaging capability for a series of MagLIF shots is presented. Predictions are given for resolution and relative image irradiance to guide experimenters in choosing the desired configuration for their experiments.

In this work, scratch and nanoindentation testing was used to determine hardness, fracture toughness, strain rate sensitivity, and activation volumes on additively manufactured graded and uniform Ni-Nb bulk specimens. Characterization showed the presence of a two phase system consisting of Ni₃Nb and Ni₆Nb₇ intermetallics. Intermetallics were multimodal in nature, having grain and cell sizes spanning from a few nanometers to 10s of micrometers. The unique microstructure resulted in impressively high hardness, up to 20 GPa in the case of the compositionally graded sample. AM methods with surface deformation techniques are a useful way to rapidly probe material properties and alloy composition space.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn-Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼ 20 %.

Abstract not provided.

While studies of urban acoustics are typically restricted to the audio range, anthropogenic activity also generates infrasound (<20 Hz, roughly at the lower end of the range of human hearing). Shutdowns related to the COVID-19 pandemic unintentionally created ideal conditions for the study of urban infrasound and low frequency audio (20-500 Hz), as closures reduced human-generated ambient noise, while natural signals remained relatively unaffected. An array of infrasound sensors deployed in Las Vegas, NV, provides data for a case study in monitoring human activity during the pandemic through urban acoustics. The array records a sharp decline in acoustic power following the temporary shutdown of businesses deemed nonessential by the state of Nevada. This decline varies spatially across the array, with stations close to McCarran International Airport generally recording the greatest declines in acoustic power. Further, declines in acoustic power fluctuate with the time of day. As only signals associated with anthropogenic activity are expected to decline, this gives a rough indication of periodicities in urban acoustics throughout Las Vegas. The results of this study reflect the city's response to the pandemic and suggest spatiotemporal trends in acoustics outside of shutdowns.

Three M-MOF-74 (M = Co, Mg, Ni) metal-organic framework (MOF) thin film membranes have been synthesized through a sensor functionalization method for the direct electrical detection of NO2. The two-step surface functionalization procedure on the glass/Pt interdigitated electrodes resulted in a terminal carboxylate group, with both steps confirmed through infrared spectroscopic analysis. This surface functionalization allowed the MOF materials to grow largely in a uniform manner over the surface of the electrode forming a thin film membrane over the Pt sensing elec-trodes. The growth of each membrane was confirmed through scanning electron microscopy (SEM) and X-ray diffraction analysis. The Ni and Mg MOFs grew as a continuous but non-defect free membrane with overlapping polycrystallites across the glass surface, whereas the Co-MOF-74 grew dis-continuously. To demonstrate the use of these MOF membranes as an NO2 gas sensor, Ni-MOF-74 was chosen as it was consistently fabricated as the best thin and homogenous membrane, as confirmed by SEM. The membrane was exposed to 5 ppm NO2 and the impedance magnitude was observed to decrease 123× in 4 h, with a larger change in impedance and a faster response than the bulk material. Importantly, the use of these membranes as a sensor for NO2 does not require them to be defect-free, but solely continuous and overlapping growth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

NaI-AlBr3 is a very appealing low melting temperature (<100 C), salt system for use as an electrochemically-active electrolyte. This system was investigated for its electrochemical and physical properties with focus to energy storage considerations. A simple phase diagram was generated; at >100 C, lower NaI concentrations had two partially miscible liquid phases, while higher NaI concentrations had solid particles. Considering the fully molten regime, electrical conductivities were evaluated over 5-25 mol% NaI and 110 C-140 C. Conductivities of 6.8-38.9 mS cm-1 were observed, increasing with temperature and NaI concentration. Effective diffusion coefficients of the I-/I3- redox species were found to decrease with both increasing NaI concentration and increasing applied potential. Regardless, oxidation current density at 3.6 V vs Na/Na+ was observed to increase with increasing NaI concentration over 5-25 mol%. Finally, the critical interface between the molten salt electrolyte and electrode materials was found to significantly affect reaction kinetics. When carbon was used instead of tungsten, an adsorbed species, most likely I2, blocked surface sites and significantly decreased current densities at high potentials. This study shows the NaI-AlBr3 system offers an attractive, low-temperature molten salt electrolyte that could be useful to many applied systems, though composition and electrode material must be considered.

Learning 3D representations that generalize well to arbitrarily oriented inputs is a challenge of practical importance in applications varying from computer vision to physics and chemistry. We propose a novel multi-resolution convolutional architecture for learning over concentric spherical feature maps, of which the single sphere representation is a special case. Our hierarchical architecture is based on alternatively learning to incorporate both intra-sphere and inter-sphere information. We show the applicability of our method for two different types of 3D inputs, mesh objects, which can be regularly sampled, and point clouds, which are irregularly distributed. We also propose an efficient mapping of point clouds to concentric spherical images, thereby bridging spherical convolutions on grids with general point clouds. We demonstrate the effectiveness of our approach in improving state-of-the-art performance on 3D classification tasks with rotated data.

Nonvolatile, molecular multiferroic devices have now been demonstrated, but it is worth giving some consideration to the issue of whether such devices could be a competitive alternative for solid‐state nonvolatile memory. For the Fe (II) spin crossover complex [Fe{H2B(pz)2}2(bipy)], where pz = tris(pyrazol‐1‐yl)‐borohydride and bipy = 2,2′‐bipyridine, voltage‐controlled isothermal changes in the electronic structure and spin state have been demonstrated and are accompanied by changes in conductance. Higher conductance is seen with [Fe{H2B(pz)2}2(bipy)] in the high spin state, while lower conductance occurs for the low spin state. Plausibly, there is the potential here for low‐cost molecular solid‐state memory because the essential molecular thin films are easily fabricated. However, successful device fabrication does not mean a device that has a practical value. Here, we discuss the progress and challenges yet facing the fabrication of molecular multiferroic devices, which could be considered competitive to silicon.

Abstract not provided.

Abstract not provided.

Electrodialysis (ED) desalination performance of different conventional and laboratoryscale ion exchange membranes (IEMs) has been evaluated by many researchers, but most of these studies used their own sets of experimental parameters such as feed solution compositions and concentrations, superficial velocities of the process streams (diluate, concentrate, and electrode rinse), applied electrical voltages, and types of IEMs. Thus, direct comparison of ED desalination performance of different IEMs is virtually impossible. While the use of different conventional IEMs in ED has been reported, the use of bioinspired ion exchange membrane has not been reported yet. The goal of this study was to evaluate the ED desalination performance differences between novel laboratory-scale bioinspired IEM and conventional IEMs by determining (i) limiting current density, (ii) current density, (iii) current efficiency, (iv) salinity reduction in diluate stream, (v) normalized specific energy consumption, and (vi) water flux by osmosis as a function of (a) initial concentration of NaCl feed solution (diluate and concentrate streams), (b) superficial velocity of feed solution, and (c) applied stack voltage per cell-pair of membranes. A laboratory-scale single stage batchrecycle electrodialysis experimental apparatus was assembled with five cell-pairs of IEMs with an active cross-sectional area of 7.84 cm2. In this study, seven combinations of IEMs (commercial and laboratory-made) were compared: (i) Neosepta AMX/CMX, (ii) PCA PCSA/PCSK, (iii) Fujifilm Type 1 AEM/CEM, (iv) SUEZ AR204SZRA/CR67HMR, (v) Ralex AMH-PES/CMH-PES, (vi) Neosepta AMX/Bare Polycarbonate membrane (Polycarb), and (vii) Neosepta AMX/Sandia novel bioinspired cation exchange membrane (SandiaCEM). ED desalination performance with the Sandia novel bioinspired cation exchange membrane (SandiaCEM) was found to be competitive with commercial Neosepta CMX cation exchange membrane.

High heat flux (>500 kW/m2) ignitions occur in scenarios involving metal fires, propellants, lightning strikes, above ground nuclear weapon use, etc. Data for material response in such environments is primarily limited to experimental programs in the 1950s and 1960s. We have recently obtained new data in this environment using concentrated solar energy. A portion of the experimental data were taken with the objective that the data be useful for model validation. To maximize the utility of the data for validation of predictive codes, additional focus is placed on repeatability of the data, reduction of uncertainties, and characterization of the environment. We illustrate here a portion of the data and methods used to assess environmental and response parameters. The data we present are novel in the flux range and materials tested, and these data constitute progress in the ability to characterize fires from high flux events.

Momentum, advection, diffusion, and turbulence are component physics relating to fire simulation tools like computational fluid dynamics (CFD). Magnetic Resonance Velocimetry and Magnetic Resonance Concentration MRV/MRC techniques can produce heretofore unrivaled detailed measurements of three-component velocity and concentration fields in turbulent flows. This study exhibits 3D flow comparisons between velocity and concentration fields obtained using MRC/MRV and SIERRA/Fuego for an urban geometry based on a section of downtown Oklahoma City. A 1:2500 scale water flow scenario provides 0.8 mm resolution data. Various techniques are employed to quantify the accuracy of the simulation results. The techniques all generally suggest a good comparison between the model and experiments throughout the compared volume. The selected metrics provide benchmark accuracy measures that can be used to indicate quantitative accuracy of the simulations, as well as for targets for future simulation improvements.

We have completed a series of both vented and sealed cookoff experiments of black powder and smokeless powder in our Sandia Instrumented Thermal Ignition (SITI) apparatus at bulk densities of 1078 and 729 kg/m3, respectively. The confining aluminum cylinder was ramped from room temperature to a set point temperature and then held at the setpoint temperature until ignition. The setpoint temperatures varied between 495 to 523 K for the black powder and 401 to 412 K for the more sensitive smokeless powder. The vented experiments show a significant delay in thermal ignition, indicating that the ignition is dependent on pressure. Post experimental debris shows greater violence for our smokeless powder experiments than our black powder experiments. A simplified universal cookoff model (UCM) was calibrated using the black powder and smokeless powder SITI data and used to predict pressurization and thermal ignition. The current work presents the first calibration of the UCM with a double base propellant. This work also presents the first pressure-dependent cookoff model for black powder and smokeless powder.

The impact on the final morphology of ceria (CeO2) nanoparticles made from different precursors (commercial: cerium acetate/nitrate) and in house: cerium tri(methylsilyl)amide (Ce-TMSA)) via a microwave solid state reaction has been determined. In all instances, powder X-ray diffraction indicated that the cubic fluorite CeO2 phase (PDF# 04–004-9150, with the space group Fm-3 m) had formed. Scanning electron microscopy (SEM) images revealed spherical nanoparticles were produced from the Ce-TMSA precursor. The commercial acetate and nitrate precursors produced particles with irregular morphology. The roles of the precursor decomposition and binding energy in the synthesis of the nanocrystals with various morphologies, as well as a possible growth mechanism, were evaluated based on experimental and computational data. The formation of spherical shaped nanoparticles was determined to be due to the preferential single-step decomposition of the Ce-TMSA as well as the low activation energy to overcome decomposition. Due to the complicated decomposition of the commercial precursors and high activation energy the resulting particles adopted an irregular morphology. Highly uniform samarium doped ceria (SmxCe1-xO2-δ) nanospheres were also synthesized from Ce-TMSA and samarium tri(methylsilyl)amide (Sm-TMSA). The effects of reaction time and temperature, on the final morphology were observed through SEM. The rapid single-step decomposition of TMSA-based precursors as observed through thermogravimetric analysis (TGA) and confirmed through the calculation of potential energy surfaces and binding energies from density functional theory (DFT) calculations, indicated that nanoparticle formation follows LaMer’s classical nucleation theory.

The dynamic behavior of an elastic peridynamic material with a nonconvex bond potential is studied. In spite of the material’s inherently unstable nature, initial value problems can be solved using essentially the same techniques as with conventional materials, both analytically and numerically. In a suitably constructed material model, small perturbations grow exponentially over time until the material fails. The time for this growth is computed explicitly for a stretching bar that passes from the stable to the unstable phase of the material model. This time to failure represents an incubation time for the nucleation of a crack. The finiteness of the failure time in effect creates a rate dependence in the failure properties of the material. Thus, the unstable nature of the elastic material leads to a rate effect even though it does not contain any terms that explicitly include a strain rate dependence.

Modern electrochemical energy conversion devices require more advanced proton conductors for their broad applications. Phosphonated polymers have been proposed as anhydrous proton conductors for fuel cells. However, the anhydride formation of phosphonic acid functional groups lowers proton conductivity and this prevents the use of phosphonated polymers in fuel cell applications. Here, we report a poly(2,3,5,6-tetrafluorostyrene-4-phosphonic acid) that does not undergo anhydride formation and thus maintains protonic conductivity above 200 °C. We use the phosphonated polymer in fuel cell electrodes with an ion-pair coordinated membrane in a membrane electrode assembly. This synergistically integrated fuel cell reached peak power densities of 1,130 mW cm−2 at 160 °C and 1,740 mW cm−2 at 240 °C under H2/O2 conditions, substantially outperforming polybenzimidazole- and metal phosphate-based fuel cells. Our result indicates a pathway towards using phosphonated polymers in high-performance fuel cells under hot and dry operating conditions.

Photovoltaic (PV) soiling profiles exhibit a sawtooth shape, where cleaning events and soiling deposition periods alternate. Generally, the rate at which soiling accumulates is assumed to be constant within each deposition period. In reality, changes in rates can occur because of sudden variations in climatic conditions, e.g., dust storms or prolonged periods of rain. The existing models used to extract the soiling profile from the PV performance data might fail to capture the change points and occasionally estimate incorrect soiling profiles. This work analyzes how the introduction of change points can be beneficial for soiling extraction. Data from nine soiling stations and a 1-MW site were analyzed by using piecewise regression and three change point detection algorithms. The results showed that accounting for change points can provide significant benefits to the modeling of soiling even if not all the change point algorithms return the same improvements. Considering change points in historical trends is found to be particularly important for studies aiming to optimize cleaning schedules.

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

Semiconductor-insulator interfaces play an important role in the reliability of integrated devices; however, the impact of these interfaces on the physical mechanisms related to single-event effects has not been previously reported. We present experimental data that demonstrate that single-event charge collection can be impacted by changes in interface quality. The experimental data, combined with simulations, show that single-event response may depend on surface recombination at interface defects. The effect depends on strike location and increases with increasing linear energy transfer (LET). Surface recombination can affect single-event charge collection for interfaces with a surface recombination velocity (SRV) of 1000 cm/s and is a dominant charge collection mechanism with SRV > 10^{5} cm/s.

Uranyl ion, UO22+, and its aqueous complexes with organic and inorganic ligands can be the dominant species for uranium transport on the Earth surface or in a nuclear waste disposal system if an oxidizing condition is present. As an important biodegradation product, oxalate, C2O42−, is ubiquitous in natural environments and is known for its ability to complex with the uranyl ion. Oxalate can also form solid phases with uranyl ion in certain environments thus limiting uranium migration. Therefore, the determination of stability constants for aqueous and solid uranyl oxalate complexes is important not only to the understanding of uranium mobility in natural environments, but also to the performance assessment of nuclear waste disposal. Here we developed a thermodynamic model for the UO22+-Na+-H+-Cl--ClO4--C2O42--NO3--H2O system to ionic strength up to ∼11 mol•kg−1. We constrained the stability constants for UO2C2O4(aq) and UO2(C2O4)22− at infinite dilution based on our evaluation of the literature data over a wide range of ionic strengths up to ∼11 mol•kg−1. We also obtained the solubility constants at infinite dilution for solid uranyl oxalates, UO2C2O4•3H2O, based on the solubility data over a wide range of ionic strengths. The developed model will enable for the accurate stability assessment of oxalate complexes affecting uranium mobility under a wide range of conditions including those in deep geological repositories.

To increase situational awareness and support evidence-based policymaking, we formulated a mathematical model for coronavirus disease transmission within a regional population. This compartmental model accounts for quarantine, self-isolation, social distancing, a nonexponentially distributed incubation period, asymptomatic persons, and mild and severe forms of symptomatic disease. We used Bayesian inference to calibrate region-specific models for consistency with daily reports of confirmed cases in the 15 most populous metropolitan statistical areas in the United States. We also quantified uncertainty in parameter estimates and forecasts. This online learning approach enables early identification of new trends despite considerable variability in case reporting.

This work explores the effect of the ill-posed problem on uncertainty quantification for motion estimation using digital image correlation (DIC) (Sutton et al. [2009]). We develop a correction factor for standard uncertainty estimates based on the cosine of the angle between the true motion and the image gradients, in an integral sense over a subregion of the image. This correction factor accounts for variability in the DIC solution previously unaccounted for when considering only image noise, interpolation bias, contrast, and the software settings such as subset size and spacing.

Temperature- and irradiation-assisted failure mechanisms in miscible phase boundaries are clarified via atomistic calculations. We first establish the temperature-dependent brittle-to-ductile transition in U–Zr miscible phase boundaries. Our results confirm that these boundaries are mostly brittle at low temperatures and ductile at elevated temperatures. We then investigate the changes induced by irradiation on the fracture mechanisms in such phase boundaries. The irradiation-induced defect accumulation follows a multi-stage process that starts with the accumulation of isolated small dislocation loops before transitioning to the saturation and growth of larger dislocation loops and end up with a reorganization into forest dislocations. The accumulation of loops is the primary feature to participate in the delineation between brittle and ductile interfacial fracture in irradiated phase boundaries. At low damage levels, radiation defect interactions with the crack tip are limited and U–Zr miscible boundaries fail through the emission of dislocations ahead of the crack tip followed by brittle cleavage in agreement with the classical Griffith’s criterion for crack stability. At higher damage levels, the failure mode transitions from brittle crack growth to ductile void growth. In this case, the microcrack tip is blunted by the high density of pre-existing, radiation-induced defects in the vicinity of the crack. This interaction leads to the development and growth of a cavity at the interface as opposed to interfacial cleavage.

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) is critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Two-dimensional, steady-state analyses were performed on unit cells with centered, circular pores using a second-order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data. This approach can be readily generalized to study a wide variety of porous solids from ranging from structures at the nanoscale—such as nanocarbon tubes—to structures at macrolevel scales—such as geological features.

The IEC 61853 photovoltaic (PV) module energy rating standard requires measuring module power (and hence, efficiency) over a matrix of irradiance and temperature conditions. These matrix points represent nearly the full range of operating conditions encountered in the field in all but the most extreme locations and create an opportunity to develop alternative approaches for calculating system performance. In this article, a new PV module efficiency model is presented and compared with five published models using matrix data collected from four different PV module types. The results of the comparative analysis demonstrated that the new model improves on the existing ones exhibiting root-mean-square errors in normalized efficiency well below 0.01 for all cases and PV modules. The analysis also highlighted its ability to interpolate and extrapolate performance between and beyond measured matrix points of irradiance and temperature, establishing it as a robust yet relatively simple model for several applications that are detailed throughout this article.

Imbalance in the autonomic nervous system can lead to orthostatic intolerance manifested by dizziness, lightheadedness, and a sudden loss of consciousness (syncope); these are common conditions, but they are challenging to diagnose correctly. Uncertainties about the triggering mechanisms and the underlying pathophysiology have led to variations in their classification. This study uses machine learning to categorize patients with orthostatic intolerance. We use random forest classification trees to identify a small number of markers in blood pressure, and heart rate time-series data measured during head-up tilt to (a) distinguish patients with a single pathology and (b) examine data from patients with a mixed pathophysiology. Next, we use Kmeans to cluster the markers representing the time-series data. We apply the proposed method analyzing clinical data from 186 subjects identified as control or suffering from one of four conditions: postural orthostatic tachycardia (POTS), cardioinhibition, vasodepression, and mixed cardioinhibition and vasodepression. Classification results confirm the use of supervised machine learning. We were able to categorize more than 95% of patients with a single condition and were able to subgroup all patients with mixed cardioinhibitory and vasodepressor syncope. Clustering results confirm the disease groups and identify two distinct subgroups within the control and mixed groups. The proposed study demonstrates how to use machine learning to discover structure in blood pressure and heart rate time-series data. The methodology is used in classification of patients with orthostatic intolerance. Diagnosing orthostatic intolerance is challenging, and full characterization of the pathophysiological mechanisms remains a topic of ongoing research. This study provides a step toward leveraging machine learning to assist clinicians and researchers in addressing these challenges. [Figure not available: see fulltext.].

The variational multiscale (VMS) formulation is used to develop residual-based VMS large eddy simulation (LES) models for Rayleigh-Bénard convection. The resulting model is a mixed model that incorporates the VMS model and an eddy viscosity model. The Wall-Adapting Local Eddy-viscosity (WALE) model is used as the eddy viscosity model in this work. The new LES models were implemented in the finite element code Drekar. Simulations are performed using continuous, piecewise linear finite elements. The simulations ranged from Ra=106 to Ra=1014 and were conducted at Pr=1 and Pr=7. Two domains were considered: a two-dimensional domain of aspect ratio 2 with a fluid confined between two parallel plates and a three-dimensional cylinder of aspect ratio 1/4. The Nusselt number from the VMS results is compared against three dimensional direct numerical simulations and experiments. In all cases, the VMS results are in good agreement with existing literature.

The modern scientific process often involves the development of a predictive computational model. To improve its accuracy, a computational model can be calibrated to a set of experimental data. A variety of validation metrics can be used to quantify this process. Some of these metrics have direct physical interpretations and a history of use, while others, especially those for probabilistic data, are more difficult to interpret. In this work, a variety of validation metrics are used to quantify the accuracy of different calibration methods. Frequentist and Bayesian perspectives are used with both fixed effects and mixed-effects statistical models. Through a quantitative comparison of the resulting distributions, the most accurate calibration method can be selected. Two examples are included which compare the results of various validation metrics for different calibration methods. It is quantitatively shown that, in the presence of significant laboratory biases, a fixed effects calibration is significantly less accurate than a mixed-effects calibration. This is because the mixed-effects statistical model better characterizes the underlying parameter distributions than the fixed effects model. The results suggest that validation metrics can be used to select the most accurate calibration model for a particular empirical model with corresponding experimental data.

Sensitivity analysis (SA) is en route to becoming an integral part of mathematical modeling. The tremendous potential benefits of SA are, however, yet to be fully realized, both for advancing mechanistic and data-driven modeling of human and natural systems, and in support of decision making. In this perspective paper, a multidisciplinary group of researchers and practitioners revisit the current status of SA, and outline research challenges in regard to both theoretical frameworks and their applications to solve real-world problems. Six areas are discussed that warrant further attention, including (1) structuring and standardizing SA as a discipline, (2) realizing the untapped potential of SA for systems modeling, (3) addressing the computational burden of SA, (4) progressing SA in the context of machine learning, (5) clarifying the relationship and role of SA to uncertainty quantification, and (6) evolving the use of SA in support of decision making. An outlook for the future of SA is provided that underlines how SA must underpin a wide variety of activities to better serve science and society.

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

Technologies based on non-equilibrium, low-temperature plasmas are ubiquitous in today’s society. Plasma modeling plays an essential role in their understanding, development and optimization. An accurate description of electron and ion collisions with neutrals and their transport is required to correctly describe plasma properties as a function of external parameters. LXCat is an open-access, web-based platform for storing, exchangig and manipulating data needed for modeling the electron and ion components of non-equilibrium, low-temperature plasmas. The data types supported by LXCat are electron- and ion-scattering cross-sections with neutrals (total and differential), interaction potentials, oscillator strengths, and electron- and ion-swarm/transport parameters. Online tools allow users to identify and compare the data through plotting routines, and use the data to generate swarm parameters and reaction rates with the integrated electron Boltzmann solver. In this review, the historical evolution of the project and some perspectives on its future are discussed together with a tutorial review for using data from LXCat.

The development of optical systems operating in the far ultraviolet range (FUV, X= 100-200 nm) is limited by the efficiency of passivated aluminum (Al) mirrors. Although it is presently possible to obtain high-reflectivity FUV mirrors through physical vapor deposition, the process involves deposition with substrates at high temperatures, which is technically challenging for large optical elements. A novel passivation procedure for bare Al mirrors is reported. The treatment consisted of using a low-temperature electron-beam generated plasma produced in a gas mixture of Ar and SF6 to etch away the native oxide layer from the Al flm, while simultaneously promoting the generation of a thin aluminum tri-fuoride (AlF3) layer on the Al surface. In the first section we analyze the effect of varying both ion energy and SF6 concentration on the FUV reflectance, thickness, composition, and surface morphology of the resulting AlF3 protective layers. In the second section, the reflectivity of samples is optimized at selected important FUV wavelengths for astronomical observations. Notably, samples attained state-of-the-art reflectances of 75% at 108.5 nm (He Lyman y), 91% at 121.6nm(HLyman a), 90% at 130.4nm (OI), and of 95% at 155.0 nm (C IV). The stability over time of these passivated mirrors is also investigated.

Abstract: Multimodal in-situ experiments are the wave of the future, as this approach will permit multispectral data collection and analysis during real-time nanoscale observation. In contrast, the evolution of technique development in the electron microscopy field has generally trended toward specialization and subsequent bifurcation into more and more niche instruments, creating a challenge for reintegration and backward compatibility for in-situ experiments on state-of-the-art microscopes. We do not believe this to be a requirement in the field; therefore, we propose an adaptive instrument that is designed to allow nearly simultaneous collection of data from aberration-corrected transmission electron microscopy (TEM), probe-corrected scanning transmission electron microscopy, ultrafast TEM, and dynamic TEM with a flexible in-situ testing chamber, where the entire instrument can be modified as future technologies are developed. The value would be to obtain a holistic understanding of the underlying physics and chemistry of the process-structure–property relationships in materials exposed to controlled extreme environments. Such a tool would permit the ability to explore, in-situ, the active reaction mechanisms in a controlled manner emulating those of real-world applications with nanometer and nanosecond resolution. If such a powerful tool is developed, it has the potential to revolutionize our materials understanding of nanoscale mechanisms and transients. Graphical Abstract: [Figure not available: see fulltext.].

The choice of model form used to represent the anisotropic yield response of metals can depend strongly on the type and amount of data available for calibration. This two-part contribution considers the calibration (part I) of three yield functions: von Mises, Hill-48 and Yld2004-18p by Barlat and co-workers. This is followed by model verification exercises (part II). The material used was a 7079 aluminum alloy extruded tube. The calibration data were measurements of yield stress and Lankford ratio from uniaxial tension specimens cut along 12 orientations. Given that the tube was relatively thick-walled, some of the orientations included through-thickness components. This allowed the calibrations to be based exclusively on test data, without the need for parameter assumptions or supplemental crystal plasticity calculations. The Yld2004-18p function provided the best fit to the data available due to its 18 anisotropy parameters plus an unspecified exponent, compared to the quadratic Hill function with 6 anisotropy parameters and to the isotropic von Mises function. Whereas the Yld2004-18p function did not warrant further exploration due to the excellent fit it provided, the results showed that care must be taken when using Hill’s function. Finally, due to its parametrization with only 6 anisotropy parameters, it can significantly misrepresent the yield behavior depending on the calibration data used, possibly rendering it less desirable than a simple isotropic function in some applications.

Anodic stripping voltammetry (ASV) has been widely used for the detection of several heavy metal ions in neutral and acidic solution, in many cases employing electrodes and/or solutions incorporating Bi. In this work we demonstrate that Bi(OH)4− ion concentration can be measured in highly alkaline 8.5 M KOH solution using ASV. The addition of Pb in similar concentrations to the Bi(OH)4− being measured is shown to improve both the sensitivity and precision of the method. When the Pb additive is employed, a formal limit of detection of 8.5 ppb is achieved, compared to 17.3 ppb when the Pb additive is not used. Due to the use of Bi additives in alkaline battery chemistries, it follows that separators which limit Bi(OH)4− diffusion into the bulk electrolyte and away from the electrodes are of interest. For this purpose, we utilize ASV to determine Bi(OH)4− diffusion rates through Celgard 3501, cellophane 350P00, and Nafion 211. Bi(OH)4− crossover rates, as determined by ASV, are shown to be repeatable and consistent with expectations from the known separator structure.

Shortages of N95 respirators for use by medical personnel have driven consideration of novel conservation strategies, including decontamination for reuse and extended use. Decontamination methods listed as promising by the Centers for Disease Control and Prevention (CDC) (vaporous hydrogen peroxide (VHP), wet heat, ultraviolet irradiation (UVI)) and several methods considered for low resource environments (bleach, isopropyl alcohol and detergent/soap) were studied for two commonly used surgical N95 respirators (3M™ 1860 and 1870+ Aura™). Although N95 filtration performance depends on the electrostatically charged electret filtration layer, the impact of decontamination on this layer is largely unexplored. As such, respirator performance following decontamination was assessed based on the fit, filtration efficiency, and pressure drop, along with the relationship between (1) surface charge of the electret layer, and (2) elastic properties of the straps. Decontamination with VHP, wet heat, UVI, and bleach did not degrade fit and filtration performance or electret charge. Isopropyl alcohol and soap significantly degraded fit, filtration performance, and electret charge. Pressure drop across the respirators was unchanged. Modest degradation of N95 strap elasticity was observed in mechanical fatigue testing, a model for repeated donnings and doffings. CDC recommended decontamination methods including VHP, wet heat, and UV light did not degrade N95 respirator fit or filtration performance in these tests. Extended use of N95 respirators may degrade strap elasticity, but a loss of face seal integrity should be apparent during user seal checks. NIOSH recommends performing user seal checks after every donning to detect loss of appropriate fit. Decontamination methods which degrade electret charge such as alcohols or detergents should not be used on N95 respirators. The loss of N95 performance due to electret degradation would not be apparent to a respirator user or evident during a negative pressure user seal check.

Electric power is crucial to the function of other infrastructures, as well as to the stability of the economy and the social order. Disruption of commercial electric power service, even for brief periods of time, can create significant consequences to the function of other sectors, and make living in some environments untenable. This analysis, conducted in 2017 for the United States Department of Energy (DOE) as part of the Grid Modernization Laboratory Consortium (GMLC) Initiative, focuses on describing the function of each of the other infrastructure sectors and subsectors, with an eye towards those elements of these sectors that depend on primary electric power service through the commercial electric power grid. It leverages the experience of Sandia analysts in analyzing historical disruptive events, and from the development of capabilities designed to identify the physical, logical, and geographic connectivity between infrastructures. The analysis goes on to identify alternatives for the provision of primary electric power service, and the redundancy of said alternatives, to provide a picture of the sector’s ability to withstand an extended disruption.

Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for pure components in liquid, gas and super critical phases. The ANNs were tested on an experimental database of 6625 self-diffusion constants for 118 different chemical compounds. The presence of multiple phases results in a heavy skew in the distribution of diffusion constants and multiple approaches were used to address this challenge. First, an ANN was developed with the raw diffusion values to assess what the main drawbacks of this direct method were. The first approach for improving the predictions involved taking the log 10 of diffusion to provide a more uniform distribution and reduce the range of target output values used to develop the ANN. The second approach involved developing individual ANNs for each phase using the raw diffusion values. Results show that the log transformation leads to a model with the best self-diffusion constant predictions and an overall average absolute deviation (AAD) of 6.56%. The resultant ANN is a generalized model that can be used to predict diffusion across all three phases and over a diverse group of compounds. The importance of each input feature was ranked using a feature addition method revealing that the density of the compound has the largest impact on the ANN prediction of self-diffusion constants in pure compounds.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report documents a set of simplified models to predict pipeline collapse under nuclear pressure loading. After a review of pipeline design literature, a set of simple expressions have been selected to represent an approximation of the threshold pressure for failure from cross- sectional yielding, cross-section buckling, and longitudinal buckling. These expressions provide a first-order approximation on load levels needed to achieve collapse. As a demonstration, the collapse pressure for a set of representative pipelines are calculated. Estimated pressure fields are also computed for a set of nuclear detonations, providing estimates of the ground range limit for pipeline collapse.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.