Time resolved liquid and vapor fields of dodecane and oxymethylene ethers are measured from Spray A-3 and Spray D using high speed Rayleigh scattering and diffuse back illumination at the Engine Combustion Network (ECN) Spray A condition of 900 K and 22.8 kg/m3. Global quantities including mixture fraction, vapor and liquid penetration, as well as spreading angle are measured. The mixture fraction fields and vapor penetration profiles are well predicted by the 1-D Musculus-Kattke model. The mixture fraction field and vapor penetration from Spray A-3 are similar to those measured from Spray A in previous works. Spray D exhibits higher mixture fraction fields and vapor penetration due to its larger nozzle diameter. The quasi-steady mixture fraction fields from these injectors scale well when distance from the injector is normalized by the nozzle diameter. The turbulent dissipation structures were also analyzed based on the orientation, thickness, and magnitude of the mixing layers. The orientation and thickness are similar to other measurements in atmospheric gas jets, while the magnitude is lower. The thickness and magnitude are subject to uncertainties due to limitations in the imaging resolution of the system but still provide an order of magnitude as a useful reference for comparison against computational fluid dynamic simulations.
Mishra, Umakant; Shi, Zheng; Hoffman, Forrest M.; Xu, Min; Allison, Steven D.; Zhou, Jizhong; Randerson, James T.
Soil carbon (C) responses to environmental change represent a major source of uncertainty in the global C cycle. Feedbacks between soil C stocks and climate drivers could impact atmospheric CO2 levels, further altering the climate. Here, we assessed the reliability of Earth system model (ESM) predictions of soil C change using the Coupled Model Intercomparison Project phases 5 and 6 (CMIP5 and CMIP6). ESMs predicted global soil C gains under the high emission scenario, with soils taking up 43.9 Pg (95% CI: 9.2–78.5 Pg) C on average during the 21st century. The variation in global soil C change declined significantly from CMIP5 (with average of 48.4 Pg [95% CI: 2.0–94.9 Pg] C) to CMIP6 models (with average of 39.3 Pg [95% CI: 23.9–54.7 Pg] C). For some models, a small C increase in all biomes contributed to this convergence. For other models, offsetting responses between cold and warm biomes contributed to convergence. Although soil C predictions appeared to converge in CMIP6, the dominant processes driving soil C change at global or biome scales differed among models and in many cases between earlier and later versions of the same model. Random Forest models, for soil carbon dynamics, accounted for more than 63% variation of the global soil C change predicted by CMIP5 ESMs, but only 36% for CMIP6 models. Although most CMIP6 models apparently agree on increased soil C storage during the 21st century, this consensus obscures substantial model disagreement on the mechanisms underlying soil C response, calling into question the reliability of model predictions.
The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. We manually picked ∼17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (VP), S-wave velocity (VS), and VP = VS ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ∼700 m in areas and generally align with geologic and structural expectations. VP and VS are relatively reduced near mapped and inferred faults, with the most prominent lower VP and VS zone around the densest collection of faults. We image VP = VS ratios ranging from ∼1.5 to ∼2.4, the extremes of which occur at a depth of ∼100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ∼100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher VP = VS areas and drier rock in lower VP = VS areas. A relatively higher VP = VS area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.
In superconducting systems in which inversion and time-reversal symmetry are simultaneously broken the critical current for positive and negative current bias can be different. For superconducting systems formed by Josephson junctions (JJs) this effect is termed Josephson diode effect. In this paper, we study the Josephson diode effect for a superconducting quantum interference device (SQUID) formed by a topological JJ with a 4π-periodic current-phase relationship and a topologically trivial JJ. We show how the fractional Josephson effect manifests in the Josephson diode effect with the application of a magnetic field and how tuning properties of the trivial SQUID arm can lead to diode polarity switching. We then investigate the ac response and show that the polarity of the diode effect can be tuned by varying the ac power and discuss differences between the ac diode effect of asymmetric SQUIDs with no topological JJ and SQUIDs in which one JJ is topological.
Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. We have implemented CSP in a C++ open-source library CSPlib1 using the Kokkos2 parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS. Program summary: Program Title: CSPlib CPC Library link to program files: https://doi.org/10.17632/p9gb7z54sp.1 Developer's repository link: https://github.com/sandialabs/csplib Licensing provisions: BSD 2-clause Programming language: C++ Nature of problem: Dynamical systems can involve coupled processes with a wide range of time scales. The computational singular perturbation (CSP) method offers a reformulation of these systems which enables the use of dynamically-based diagnostic tools to better comprehend the dynamics by decoupling fast and slow processes. CSPlib is an open-source software library for analyzing general ordinary differential equation (ODE) and differential algebraic equation (DAE) systems, with specialized implementations for detailed chemical kinetic ODE/DAE systems. It relies on CSP for the analysis of these systems. CSPlib has been used in gas kinetic and heterogeneous catalytic kinetic models. Solution method: CSP analysis seeks a set of basis vectors to linearly decompose the right-hand side (RHS) of a dynamical system in a manner that decouples fast and slow processes. The CSP basis vectors are often well approximated with the right eigenvectors of the RHS Jacobian. And the left basis vectors are found by the inversion of the matrix, whose columns are the CSP basis vectors. Accordingly, the right and left CSP basis vectors are orthonormal. CSP defines mode amplitudes as the projections of the left basis vectors on the RHS; the time scales as the reciprocals of the RHS Jacobian eigenvalue magnitudes; and the CSP pointers, which are the element-wise multiplication of the transpose of the right CSP basis vectors with the left CSP basis vectors. For kinetic models that can be cast as the product of a generalized stoichiometric matrix and a rate of progress vector, CSP defines the participation index, which represents the contribution of a chemical reaction to each mode. Further, it defines the slow and fast importance indices, which describe the contribution of a chemical reaction to the slow and fast dynamics of a state variable, respectively. These indices are useful in diagnostic studies of dynamical systems and the construction of simplified models. Additional comments including restrictions and unusual features: CSPlib is a portable library that carries out many CSP analyses in parallel and can be used in modern high-performance platforms.
High Energy Arcing Faults (HEAFs) are hazardous events in which an electrical arc leads to the rapid release of energy in the form of heat, vaporized metal, and mechanical force. In Nuclear Power Plants (NPPs), these events are often accompanied by loss of essential power and complicated shutdowns. To confirm the probabilistic risk analysis (PRA) methodology in NUREG/CR-6850, which was formulated based on limited observational data, the NRC led an international experimental campaign from 2014 to 2016. The results of these experiments uncovered an unexpected hazard posed by aluminum components in or near electrical equipment and the potential for unanalyzed equipment failures. Sandia National Laboratories (SNL), in support of the NRC work, collaborated with NIST, BSI, KEMA, and NRC to support the full-scale HEAF test campaign in 2023. SNL provided high speed and real time from visible and infrared video/data of tests that collected data from copper and aluminum busses from switchgears and bus-ducts. Part of SNL work was to place cameras with high-speed data collection capability at different vantage points that provide the NRC a more complete and granular view of the test events.
Zeno Power Systems, Inc. (“Zeno”) is developing the Low-Profile Electric-Propulsion Nuclear Satellite-1 (LENS-1), which includes a radioisotope thermoelectric generator powered by the radioactive decay of strontium-90. Sandia National Laboratories was retained by Zeno to conduct a preliminary safety analysis of the LENS-1 mission. This analysis was conducted to determine whether LENS-1 is expected to meet Tier I requirements of National Security Presidential Memorandum-20. This document provides the results of an initial, deterministic safety analysis to support a Preliminary Tier Determination (PTD). This PTD is a primary element of Zeno’s revised Payload Review submitted to FAA’s Office of Commercial Space Transportation for launch approval of LENS-1.
Tropical cyclones are the leading cause of major power outages in the U.S., and their effects can be devastating for communities. However, few studies have holistically examined the degree to which socio-economic variables can explain spatial variations in disruptions and reveal potential inequities thereof. Here, we apply machine learning techniques to analyze 20 tropical cyclones and predict county-level outage duration and percentage of customers losing power using a comprehensive set of weather, environmental, and socio-economic factors. Our models are able to accurately predict these outage response variables, but after controlling for the effects of weather conditions and environmental factors in the models, we find the effects of socio-economic variables to be largely immaterial. However, county-level data could be overlooking effects of socio-economic disparities taking place at more granular spatial scales, and we must remain aware of the fact that when faced with similar outage events, socio-economically vulnerable communities will still find it more difficult to cope with disruptions compared to less vulnerable ones.
Lukin, Illya V.; Sotnikov, Andrii G.; Leamer, Jacob M.; Magann, Alicia B.; Bondar, Denys I.
We present an expression for the spectral gap, opening up new possibilities for performing and accelerating spectral calculations of quantum many-body systems. We develop and demonstrate one such possibility in the context of tensor network simulations. Our approach requires only minor modifications of the widely used simple update method and is computationally lightweight relative to other approaches. We validate it by computing spectral gaps of the 2D and 3D transverse-field Ising models and find strong agreement with previously reported perturbation theory results.
This report summarizes a gap analysis resulting from a literature review and expert interviews conducted by subject matter experts from Sandia National Laboratory, Siemens, and the Electric Power Research Institute (EPRI) in Spring 2023. The gap analysis consists of two main parts: The fault-ride through (FRT) behavior of grid-forming (GFM) inverter-based resources (IBR) and the response of state-of-the-art protection relays to the fault currents and voltages from GFM IBRs.
Due to its tunable bandgap, anisotropic behavior, and superior thermoelectric properties, device applications using layered tellurene (Te) are becoming more attractive. Here, we report a thinning technique for exfoliated tellurene nanosheets using thermal annealing in an oxygen environment. We characterize different thinning parameters, including temperature and annealing time. Based on our measurements, we show that controlled layer thinning occurs in the narrow temperature range of 325-350 °C. We also show a reliable method to form β-tellurene oxide (β-TeO2), which is an emerging wide bandgap semiconductor with promising electronic and optoelectronic properties. This wide bandgap semiconductor exhibits a broad photoluminescence (PL) spectrum with multiple peaks covering the range of 1.76-2.08 eV. This PL emission, coupled with Raman spectra, is strong evidence of the formation of 2D β-TeO2. We discuss the results obtained and the mechanisms of Te thinning and β-TeO2 formation at different temperature regimes. We also discuss the optical bandgap of β-TeO2 and show the existence of pronounced excitonic effects evident by the large exciton binding energy in this 2D β-TeO2 system that reach 1.54-1.62 eV for bulk and monolayer, respectively. Our work can be utilized to have better control over the Te nanosheet thickness. It also sheds light on the formation of well-controlled β-TeO2 layered semiconductors for electronic and optoelectronic applications.
We present a new optimization-based property-preserving algorithm for passive tracer transport. The algorithm utilizes a semi-Lagrangian approach based on incremental remapping of the mass and the total tracer. However, unlike traditional semi-Lagrangian schemes, which remap the density and the tracer mixing ratio through monotone reconstruction or flux correction, we utilize an optimization-based remapping that enforces conservation and local bounds as optimization constraints. In so doing we separate accuracy considerations from preservation of physical properties to obtain a conservative, second-order accurate transport scheme that also has a notion of optimality. Moreover, we prove that the optimization-based algorithm preserves linear relationships between tracer mixing ratios. We illustrate the properties of the new algorithm using a series of standard tracer transport test problems in a plane and on a sphere.
In the rapidly evolving field of solar energy, Photovoltaic (PV) manufacturers are constantly challenged by the degradation of PV modules due to localized overheating, commonly known as hotspots. This issue not only reduce the efficiency of solar panels but, in severe cases, can lead to irreversible damage, malfunctioning, and even fire hazards. Addressing this critical challenge, our research introduces an innovative electronic device designed to effectively mitigate PV hotspots. This pioneering solution consists of a novel combination of a current comparator and a current mirror circuit. These components are uniquely integrated with an automatic switching mechanism, notably eliminating the need for traditional bypass diodes. We rigorously tested and validated this device on PV modules exhibiting both adjacent and non-adjacent hotspots. Our findings are groundbreaking: the hotspot temperatures were significantly reduced from a dangerous 55 °C to a safer 35 °C. Moreover, this intervention remarkably enhanced the output power of the modules by up to 5.3%. This research not only contributes a practical solution to a longstanding problem in solar panel efficiency but also opens new pathways for enhancing the safety and longevity of solar PV systems.
One of the most iconic of radar waveforms is the Linear FM chirp. It is well-behaved and well-understood, and has become the gold standard against which other radar waveforms are measured. It has a number of desirable attributes, but is not without some issues. It may be processed by a number of techniques with many variations. Details of the Linear FM chirp are presented and discussed in this report.
We present a mathematical framework for constructing the most general neutrino mass matrices that yield the observed spectrum of light active neutrino masses in conjunction with arbitrarily many heavy sterile neutrinos, without the need to assume a hierarchy between Dirac and Majorana mass terms. The seesaw mechanism is a byproduct of the formalism, along with many other possibilities for generating tiny neutrino masses. We comment on phenomenological applications of this approach, in particular deriving a mechanism to address the long-standing (g-2)μ anomaly in the context of the left-right symmetric model.
Understanding temperature-dependent material decomposition and structural deformation induced by combined thermal-mechanical environments is critical for safety qualification of hardware under accident scenarios. Seeing in with X-rays elucidated the physics necessary to develop X-ray strain and thermometry diagnostics for use in optically opaque environments. Two parallel thermometry schemes were explored: X-ray fluorescence and X-ray diffraction of inorganic doped ceramics– colloquially known as thermographic phosphors. Two parallel surface strain techniques–Path-Integrated Digital Image Correlation and Frequency Multiplexed Digital Image Correlation–were demonstrated. Finally, preliminary demonstration of time-resolved digital volume correlation was performed by taking advantage of limited view reconstruction techniques. Additionally, research into blended ceramic-metal coatings was critical to generating intrinsic thermographic patterns for the future combination of X-ray strain and thermometry measurements.
The report summarizes the work and accomplishments of DOE SETO funded project 36533 “Adaptive Protection and Control for High Penetration PV and Grid Resilience”. In order to increase the amount of distributed solar power that can be integrated into the distribution system, new methods for optimal adaptive protection, artificial intelligence or machine learning based protection, and time domain traveling wave protection are developed and demonstrated in hardware-in-the-loop and a field demonstration.
Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (Sn) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.
The properties of defects in n-p-n Si bipolar junction transistors (BJTs) caused by 17-MeV Si ions are investigated via current-voltage, low-frequency (LF) noise, and deep level transient spectroscopy (DLTS) measurements. Four prominent radiation-induced defects in the base-collector junction of these transistors are identified via DLTS. At least two defect levels are observed in temperature-dependent LF 1/f noise measurements, one that is similar to a prominent defect in DLTS and another that is not. Defect microstructures are discussed. Our results show that DLTS and 1/f noise measurements can provide complementary information about defects in linear bipolar devices.
Additive manufacturing (AM) maintains a wide process window that enables complex designs otherwise unattainable via conventional production technologies. However, the lack of confidence in qualifying AM parts that leverage AM process–structure–property–performance (PSPP) relationships stymies design optimization and adoption of AM. While continuing efforts to map fundamental PSPP relationships that cover the potential design space, we first need pragmatic and then long-term solutions that overcome challenges associated with qualifying AM-designed parts. Two pragmatic solutions include: (1) AM material specifications to substantiate process reproducibility, and (2) component risk categorization to associate system risk relative to part performance and required part quality. A novel qualification paradigm under development involves efficient prediction of part performance over wide-ranging PSPP relationships through targeted testing and computational simulation. This paper describes projects at Sandia National Laboratories on PSPP relationship discovery, these pragmatic approaches, and the novel qualification approach.
Abstract: The effect of moisture on the photo-oxidative degradation of polyamide-6 (PA-6) was studied by analyzing the mechanical response after two different accelerated aging procedures. In the first aging procedure, the PA-6 was only exposed to ultra-violet (UV) radiation at 60 ∘C. In the second procedure, the same duration of UV radiation was periodically interrupted while the relative humidity was raised to 100%. Diffusion-limited and nominally homogeneous degradation conditions were investigated using bulk and film specimens, respectively. Accelerated UV aging reduced the ductility of PA-6, but the additional hygrothermal exposure had no effect on the ductility or strength, indicating that humidity did not influence the photo-oxidation of PA-6. This finding contrasts with previous studies that found thermo-oxidation of PA-6 was accelerated by moisture. Graphical abstract: (Figure presented.)
Hwang, Joonsik; Karathanassis, Ioannis K.; Koukouvinis, Phoevos; Nguyen, Tuan; Tagliante, Fabien; Pickett, Lyle M.; Sforzo, Brandon A.; Powell, Christopher F.
As modern gasoline direct injection (GDI) engines utilize sophisticated injection strategies, a detailed understanding of the air-fuel mixing process is crucial to further improvements in engine emission and fuel economy. In this study, a comprehensive evaluation of the spray process of single-component iso-octane (IC8) and multi-component gasoline surrogate E00 (36 % n-pentane, 46 % iso-octane, and 18 % n-undecane, by volume) fuels was conducted using an Engine Combustion Network (ECN) Spray G injector. High-speed extinction, schlieren, and microscopy imaging campaigns were carried out under engine-like ambient conditions in a spray vessel. Experimental results including liquid/vapor penetration, local liquid volume fraction, droplet size, and projected liquid film on the nozzle tip were compared under ECN G1 (573 K, 3.5 kg/m3), G2 (333 K, 0.5 kg/m3), and G3 (333 K, 1.01 kg/m3) conditions. In addition to the experiments, preferential evaporation process of the E00 fuel was elucidated by Large–Eddy Simulations (LES). The three-dimensional liquid volume fraction measurement enabled by the computed tomographic reconstruction showed substantial plume collapse for E00 under the G2 and G3 conditions having wider plume growth and plume-to-plume interaction due to the fuel high vapor pressure. The CFD simulation of E00 showed an inhomogeneity in the way fuel components vaporized, with more volatile components carried downstream in the spray after the end of injection. The high vapor pressure of E00 also results in ∼4 μm smaller average droplet diameter than IC8, reflecting a higher rate of initial vaporization even though the final boiling point temperature is higher. Consistent with high vapor pressure, E00 had a wider plume cone angle and enhanced interaction with the wall to cover the entire surface of the nozzle tip in a film. However, the liquid fuel underwent faster evaporation, so the final projected tip wetting area was smaller than the IC8 under the flash-boiling condition.
ORIGEN is one of the main transmutation software packages used in nuclear engineering Modeling Tungsten Boride Neutronics in ORIGEN for Z-Facilityproblems. For the case of this study, tungsten borides are studied using a coupled framework between MCNP and the ORIGEN package of scale. The input used four compositions of tungsten boride: WB with natural boron- 10 abundance, WB with 80wt% B-10 per isotope of boron, WB4 with natural boron-10 abundance, and WB4 with 80wt% B-10 per isotope of boron. Isotopic inventories were produced for WB which show the time dependent change up to 2 years after a 6-Month irradiation. This will allow for further studies of the materials to assess things material composition changes, dose contribution, and waste management requirements.
Long-duration energy storage (LDES) is critical to a stable, resilient, and decarbonized electric grid. While batteries are emerging as important LDES devices, extended, high-power discharges necessary for cost-competitive LDES present new materials challenges. Focusing on a new generation of low-temperature molten sodium batteries, we explore here unique phenomena related to long-duration discharge through a well-known solid electrolyte, NaSICON. Specifically, molten sodium symmetric cells at 110 °C were cycled at 0.1 A cm−2 for 1-23 h discharges. Longer discharges led to unstable overpotentials, reduced resistances, and decreased electrolyte strength, caused by massive sodium penetration not observed in shorter duration discharges. Scanning electron microscopy informed mechanisms of sodium penetration and even “healing” during shorter-duration cycling. Importantly, these findings show that traditional, low-capacity, shorter-duration tests may not sufficiently inform fundamental materials phenomena that will impact LDES battery performance. This case highlights the importance that candidate LDES batteries be tested under pertinent long-duration conditions.
We characterize the performance of two pixelated neutron detectors: a PMT-based array that utilizes Anger logic for pixel identification and a SiPM-based array that employs individual pixel readout. The SiPM-based array offers improved performance over the previously developed PMT-based detector both in terms of uniformity and neutron detection efficiency. Each detector array uses PSD-capable plastic scintillator as a detection medium. We describe the calibration and neutron efficiency measurement of both detectors using a 137Cs source for energy calibration and a 252Cf source for calibration of the neutron response. We find that the intrinsic neutron detection efficiency of the SiPM-based array is (30.2 ± 0.9)%, which is almost twice that of the PMT-based array, which we measure to be (16.9 ± 0.1)%.
Hydrogen continues to show promise as a viable contributor to achieving energy storage goals such as energy security and decarbonization in the United States. However, many new and expanded hydrogen use applications will require identifying methods of larger-scale storage than the solutions that currently exist for smaller storage applications. One possibility is to store large quantities of gaseous hydrogen below ground level. Underground storage of other fuels such as natural gas is already currently utilized, so much of the infrastructure and basic technologies can be used as a basis for underground hydrogen storage (UHS). A few commercial UHS facilities currently exist in the United States, including salt caverns owned and operated by Air Liquide, Linde, and Conoco Philips, but UHS is still a relatively new concept that has not been widely deployed. It is necessary to understand the safety risks and hazards associated with UHS before its use can be expanded and accepted more broadly. Many of these risks are addressed through regulations, codes, and standards (RCS) issued by governing bodies and organizations with expertise in certain hazards. This report is a review of RCS documents relevant to UHS, with a particular lens on potential technical gaps in existing guidance. These gaps may be specific to the physical properties of hydrogen or due to the different technologies relevant for hydrogen vs. natural gas storage. This is meant to be a high-level review to identify relevant documents and potential gaps. Formally addressing the individual gaps identified here within the codes and standards themselves would involve a more intensive analysis and differ based on the code or standard revision processes of the various publishing organizations. Therefore, presenting specific recommendations for revising the verbiage of the documents for UHS applications is left for future work and other publications.
This report presents analyses of the AB5 and AB6 ABCOVE sodium spray fire experiments with the MELCOR code. This code simulates the progression of accident events for analysis and auditing purposes of nuclear facilities during accident conditions. Historically, the ABCOVE experiments have contributed to the validation of aerosol physics and related phenomena. Given advancements in sodium-cooled reactor designs, characterization of the sodium spray combustion may further the review and validation of newly incorporated sodium properties and physics packages, namely, the sodium equations of state (EOS) and the sodium combustion (NAC) package within MELCOR. By analyzing the AB5 and AB6 experiments with and without the NAC package, sodium specificity for spray combustion and aerosol formation as well as speciation of the combustion products are reviewed with the new packages. This effort provides code users with a demonstration of the current code capabilities. This report provides the current best practices for the NAC package as well as a discussion of any issues observed while performing the presented analyses.
This article describes the implementation of a new numerical model of the power take-off system installed in the Monterey Bay Aquarium Research Institute wave energy converter, a device developed to provide power to various oceanic research missions. The simultaneous presence of hydraulic, pneumatic, and electrical subsystems in the power take-off system represents a significant challenge in forging an accurate model able to replicate the main dynamic characteristics of the system. The validation of the new numerical model is addressed by comparing simulations with the measurements obtained during a series of bench tests. Data from the bench tests show good agreement with the numerical model. The validated model provides deeper insights into the complex nonlinear dynamics of the power take-off system and will support further performance improvements in the future.
The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal, such as fuel degradation; cladding corrosion, embrittlement, or breaching; and the creation of a flammable environment via radiolysis of water. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents details on the quantification of residual water in the Advanced Drying Cycle Simulator (ADCS), an apparatus built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. The ADCS is outfitted with thermocouples to measure the thermal response of the ADCS to simulated decay heats and internal helium fill pressures relevant to commercial drying procedures. The ADCS is also instrumented with pressure transducers to measure the pressures and vacuum levels observed during simulated commercial drying. The most unique instrumentation used for quantifying residual water in the ADCS is a Hiden Analytical HPR-30 mass spectrometer (MS), which measures gas compositions of the ADCS internal free volume, based on partial pressures calculated from relative proportions of gas molecules detected by the MS. This report details the methodology used to implement MS measurements in quantifying residual water in the ADCS. This methodology includes the calibration of the HPR-30 MS to a Buck Research Instruments CR-4 chilled mirror hygrometer, which itself is calibrated to a NIST-traceable standard. Data collected by both the MS and the chilled mirror hygrometer from water/helium mixtures ranging from 150 to 500,000 ppmv water in helium were used to generate calibration curves, establishing a source of verification of MS measured water contents. Details regarding water content measurement uncertainties are included in this report, defining the accuracy and verifiability of the HPR-30 MS in measuring residual water content in simulated dry storage canister environments.
Motivated by increasing interest in electrochemical devices that include highly alkaline electrolytes, we investigated two force fields for potassium hydroxide (KOH) at high concentrations in water. The “FNB” model uses the SPC/E water model, while the “FHM” model uses the TIP4P/2005 water model. Here, we also developed parameters to describe zincate ions in these solutions. The density and viscosity of KOH using the FHM model are in better agreement with experiment than the values from the FNB model. Comparing the properties of the zincate solutions to the available experimental data, we find that both force fields agree reasonably well, although the FHM parameters give a better prediction of the viscosity. The developed force field parameters can be used in future simulations of zincate/KOH solutions in combination with other species of interest.
Dzara, Michael J.; Campello, Arthur C.; Breidenbach, Aeryn T.; Strange, Nicholas A.; Park, James E.; Ambrosini, Andrea A.; Coker, Eric N.; Ginley, David S.; Lee, Young S.; Bell, Robert T.; Smaha, Rebecca W.
Material design is increasingly used to realize desired functional properties, and the perovskite structure family is one of the richest and most diverse: perovskites are employed in many applications due to their structural flexibility and compositional diversity. Hexagonal, layered perovskite structures with chains of face-sharing transition metal oxide octahedra have attracted great interest as quantum materials due to their magnetic and electronic properties. Ba4MMn3O12, a member of the “12R” class of hexagonal, layered perovskites, contains trimers of face-sharing MnO6 octahedra that are linked by a corner-sharing, bridging MO6 octahedron. Here, we investigate cluster magnetism in the Mn3O12 trimers and the role of this bridging octahedron on the magnetic properties of two isostructural 12R materials by systematically changing the M4+ cation from nonmagnetic Ce4+ (f0) to magnetic Pr4+ (f1). We synthesized 12R-Ba4MMn3O12 (M= Ce, Pr) with high phase purity and characterized their low-temperature crystal structures and magnetic properties. Using substantially higher purity samples than previously reported, we confirm the frustrated antiferromagnetic ground state of 12R-Ba4PrMn3O12 below TN ≈ 7.75 K and explore the cluster magnetism of its Mn3O12 trimers. Despite being atomically isostructural with 12R-Ba4CeMn3O12, the f1 electron associated with Pr4+ causes much more complex magnetic properties in 12R-Ba4PrMn3O12. In 12R-Ba4PrMn3O12, we observe a sharp, likely antiferromagnetic transition at T2 ≈ 12.15 K and an additional transition at T1 ≈ 200 K, likely in canted antiferromagnetic order. These results suggest that careful variation of composition within the family of hexagonal, layered perovskites can be used to tune material properties using the complex role of the Pr4+ ion in magnetism.
A study of wildfire danger ratings was conducted by the Meteorological Program using the adjective ratings reported to the Weather Information Management System (WIMS). The focus of this study was to identify the shortcomings of Instruction 13-212 from Kirtland Air Force Base with regards to fire danger rating. The goal was to illustrate the benefits of using data from the Sandia National Labs Remote Automated Weather Station (RAWS) to better identify the proper fire danger rating on Kirtland Air Force Base.
Phase diagrams exhibiting extended solid-solution and lenslike melting are often reproduced using ideal solutions, where ideal mixing considers a fully random configurational entropy of mixing. In the field of irreversible thermodynamics, experimental measurements of the composition variation of high-temperature electronic transport and molten-state properties suggest, however, a strong role for short-range atomic ordering in these systems. Herein, measurements of the thermopower and resistivity are reported for Cu-Ni solid solutions as a function of temperature and composition. The electronic transport properties were interpreted with an irreversible thermodynamic framework, revealing a large electronic contribution to the entropy of mixing. By considering a cluster model for the configurational entropy that uses the electronic contribution to inform the existence of ordered associates, we rationalize such a contribution of the electronic entropy with the ideal entropy of mixing commonly used to model such systems. In conclusion, these results suggest that the short-range order of the atoms plays a significant role in both solid and liquid states, even when there are no dominant intermetallic compounds in these alloys.
Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.
Custom-form factor batteries fabricated in non-conventional shapes can maximize the overall energy density of the systems they power, particularly when used in conjunction with energy dense materials (e.g., Li metal anodes and conversion cathodes). Additive manufacturing (AM), and specifically material extrusion (ME), have been shown as effective methods for producing custom-form cell components, particularly electrodes. However, the AM of several promising energy dense materials (conversion electrodes such as iron trifluoride) have yet to be demonstrated or optimized. Furthermore, the integration of multiple AM produced cell components, such as electrodes and separators, along with a custom package remains largely unexplored. In this work, iron trifluoride (FeF3) and ionogel (IG) separators are conformally printed using ME onto non-planar surfaces to enable the fabrication of custom-form Li-FeF3 batteries. To demonstrate printing on non-planar surfaces, cathodes and separators were deposited onto cylindrical rods using a 5-axis ME printer. ME printed FeF3 was shown to have performance commensurate with FeF3 cast using conventional means, both in coin cell and cylindrical rod formats, with capacities exceeding 700 mAh/g on the first cycle and ranging between 600 and 400 mAh/g over the next 50 cycles. Additionally, a ME process for printing polyvinylidene fluoride-co-hexafluoropropylene (PVDF-HFP) based IGs directly onto FeF3 is developed and enabled using an electrolyte exchange process. In coin cells, this process is shown to produce cells with similar capacity to cells built with Celgard separators out to 50 cycles, with the exception that cycling instabilities are observed during cycles 8–20. When using printed and exchanged IGs in a custom cylindrical cell package, 6 stable high-capacity cycles are achieved. Overall, this work demonstrates approaches for producing high-energy-density Li-FeF3 cells in coin and cylindrical rod formats, which are translatable to customized, arbitrary geometries compatible with ME printing and electrolyte exchange.
The most complex challenges facing the world today comprise the work of [Department of Energy’s] (DOE’s) 17 National Laboratories: [...] From furthering U.S. energy independence and leadership in clean technologies; to promoting innovation that advances U.S. economic competitiveness; to conducting research of the highest caliber in the physical, chemical, biological, materials, computational, and information sciences to advance understanding of the world around us-the Laboratories’ purview is expansive and further their contributions are indispensable.
Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to Artificial Correlation Heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w2/3, where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.
As the field of low-dimensional materials (1D or 2D) grows and more complex and intriguing structures are continuing to be found, there is an emerging need for techniques to characterize the nanoscale mechanical properties of all kinds of 1D/2D materials, in particular in their most practical state: sitting on an underlying substrate. While traditional nanoindentation techniques cannot accurately determine the transverse Young's modulus at the necessary scale without large indentations depths and effects to and from the substrate, herein an atomic-force-microscopy-based modulated nanomechanical measurement technique with Angstrom-level resolution (MoNI/ÅI) is presented. This technique enables non-destructive measurements of the out-of-plane elasticity of ultra-thin materials with resolution sufficient to eliminate any contributions from the substrate. This method is used to elucidate the multi-layer stiffness dependence of graphene deposited via chemical vapor deposition and discover a peak transverse modulus in two-layer graphene. While MoNI/ÅI has been used toward great findings in the recent past, here all aspects of the implementation of the technique as well as the unique challenges in performing measurements at such small resolutions are encompassed.
Continued dependence on crude oil and natural gas resources for fossil fuels has caused global atmospheric carbon dioxide (CO2) emissions to increase to record-setting proportions. There is an urgent need for efficient and inexpensive carbon sequestration systems to mitigate large-scale CO2 emissions from industrial flue gas. Carbonic anhydrase (CA) has shown high potential for enhanced CO2 capture applications compared to conventional absorption-based methods currently utilized in various industrial settings. This study aims to understand structural aspects that contribute to the stability of CA enzymes critical for their applications in industrial processes, which require the ability to withstand conditions different from their native environments. Here, we evaluated the thermostability and enzyme activity of mesophilic and thermophilic CA variants at different temperature conditions and in the presence of atmospheric gas pollutants like nitrogen oxides (NOx) and sulphur oxides (SOx). Based on our enzyme activity assays and molecular dynamics simulations, we see increased conformational stability and CA activity levels in thermostable CA variants incubated week-long at different temperature conditions. The thermostable CA variants also retained high levels of CA activity despite changes in solution pH due to increasing NOx and SOx concentrations. Furthermore, a loss of CA activity was observed only at high concentrations of NOx/SOx that possibly can be minimized with appropriate buffered solutions.
Manganese dioxide is a promising cathode material for energy storage applications because of its high redox potential, large theoretical energy density, abundance, and low cost. It has been shown that the performance of MnO2 electrodes in rechargeable alkaline Zn/MnO2 batteries could be improved by nanostructuring and by increasing the concentration of defects in MnO2. However, the underlying mechanism of this improvement is not completely clear. We used an ab initio density functional computational approach to investigate the influence of nanostructuring and crystal defects on the electrochemical properties of the MnO2 cathode material. The mechanism of electrochemical discharge of MnO2 in Zn/MnO2 batteries was studied by modeling the process of H ion insertion into the structures of pyrolusite, ramsdellite, and nsutite polymorphs containing oxygen vacancies, cation vacancies, and open surfaces. Our calculations showed that the binding energies of H ions inserted into the structures of MnO2 polymorphs were strongly affected by the presence of surfaces and bulk defects. In particular, we found that the energies of H ions inserted under the surfaces and attached to the surfaces of MnO2 crystals were significantly lower than those for bulk MnO2. Furthermore, the results of our study provide an explanation for the influence of crystal defects and nanostructuring on the electrochemical reactivity of MnO2 cathodes in rechargeable alkaline Zn/MnO2 batteries.
Mode-locked vertical external cavity semiconductor lasers are a unique class of nonlinear dynamical systems driven far from equilibrium. We present a novel, to the best of our knowledge, experimental result, supported by rigorous microscopic simulations, of two coexisting mode-locked V-cavity configurations sourced by a common gain medium and operating as independent channels at angle controlled separated wavelengths. Microscopic simulations support pulses coincident on the common gain chip extracting photons from a nearby pair of coexisting kinetic holes burned in the carrier distributions.
Traditional methods of shielding fragile goods and human tissues from impact energy rely on isotropic foam materials. The mechanical properties of these foams are inferior to an emerging class of metamaterials called plate lattices, which have predominantly been fabricated in simple 2.5-dimensional geometries using conventional methods that constrain the feasible design space. In this work, additive manufacturing is used to relax these constraints and realize plate lattice metamaterials with nontrivial, locally varying geometry. The limitations of traditional computer-aided design tools are circumvented and allow the simulation of complex buckling and collapse behaviors without a manual meshing step. By validating these simulations against experimental data from tests on fabricated samples, sweeping exploration of the plate lattice design space is enabled. Numerical and experimental tests demonstrate plate lattices absorb up to six times more impact energy at equivalent densities relative to foams and shield objects from impacts ten times more energetic while transmitting equivalent peak stresses. In contrast to previous investigations of plate lattice metamaterials, designs with nonuniform geometric prebuckling in the out-of-plane direction is explored and showed that these designs exhibit 10% higher energy absorption efficiency on average and 25% higher in the highest-performing design.
The effect of doping concentration on the temperature performance of the novel split-well resonant-phonon (SWRP) terahertz quantum-cascade laser (THz QCL) scheme supporting a clean 4-level system design was analyzed using non-equilibrium Green’s functions (NEGF) calculations. Experimental research showed that increasing the doping concentration in these designs led to better results compared to the split-well direct-phonon (SWDP) design, which has a larger overlap between its active laser states and the doping profile. However, further improvement in the temperature performance was expected, which led us to assume there was an increased gain and line broadening when increasing the doping concentration despite the reduced overlap between the doped region and the active laser states. Through simulations based on NEGF calculations we were able to study the contribution of the different scattering mechanisms on the performance of these devices. We concluded that the main mechanism affecting the lasers’ temperature performance is electron-electron (e-e) scattering, which largely contributes to gain and line broadening. Interestingly, this scattering mechanism is independent of the doping location, making efforts to reduce overlap between the doped region and the active laser states less effective. Optimization of the e-e scattering thus could be reached only by fine tuning of the doping density in the devices. By uncovering the subtle relationship between doping density and e-e scattering strength, our study not only provides a comprehensive understanding of the underlying physics but also offers a strategic pathway for overcoming current limitations. This work is significant not only for its implications on specific devices but also for its potential to drive advancements in the entire THz QCL field, demonstrating the crucial role of e-e scattering in limiting temperature performance and providing essential knowledge for pushing THz QCLs to new temperature heights.
This project has 2 major thrusts and 2 minor tasks – One of the major ones is a scoping study for improving uranium isotopics by analyzing higher energy, lower yield U-235 gamma emissions. The second major one benefits past projects such as auto-enrichment, and peak-based analysis routines previously added to GADRAS-DRF [1]. The two minor tasks involve an in-person or remote GADRAS-DRF training for the IAEA, and characterizing detectors of interest.
Porous liquids (PLs), which are solvent-based systems that contain permanent porosity due to the incorporation of a solid porous host, are of significant interest for the capture of greenhouse gases, including CO2. Type 3 PLs formed by using metal-organic frameworks (MOFs) as the nanoporous host provide a high degree of chemical turnability for gas capture. However, pore aperture fluctuation, such as gate-opening in zeolitic imidazole framework (ZIF) MOFs, complicates the ability to keep the MOF pores available for gas adsorption. Therefore, an understanding of the solvent molecular size required to ensure exclusion from MOFs in ZIF-based Type 3 PLs is needed. Through a combined computational and experimental approach, the solvent-pore accessibility of exemplar MOF ZIF-8 was examined. Density functional theory (DFT) calculations identified that the lowest-energy solvent-ZIF interaction occurred at the pore aperture. Experimental density measurements of ZIF-8 dispersed in various-sized solvents showed that ZIF-8 adsorbed solvent molecules up to 2 Å larger than the crystallographic pore aperture. Density analysis of ZIF dispersions was further applied to a series of possible ZIF-based PLs, including ZIF-67, −69, −71(RHO), and −71(SOD), to examine the structure-property relationships governing solvent exclusion, which identified eight new ZIF-based Type 3 PL compositions. Solvent exclusion was driven by pore aperture expansion across all ZIFs, and the degree of expansion, as well as water exclusion, was influenced by ligand functionalization. Using these results, a design principle was formulated to guide the formation of future ZIF-based Type 3 PLs that ensures solvent-free pores and availability for gas adsorption.
Cytokines and acute-phase proteins are promising biomarkers for inflammatory disease. Despite its potential, early diagnosis based on these biomarkers remains challenging without technology enabling highly sensitive protein detection immediately after sample collection, because of the low abundance and short half-life of these proteins in bodily fluids. Enzyme-linked immunosorbent assay (ELISA) is a gold-standard method for such protein analysis, but it often requires labor-intensive and time-consuming sample handling and as well as a bulky benchtop platereader, limiting its utility in the clinical site. We developed a portable microfluidic immunoassay device capable of sensitive, quantitative, and high-throughput protein detection at point-of-need. The portable microfluidic system performs eight magnetic bead-based sandwich immunoassays from raw samples in 40 min. An innovative bead actuation strategy was incorporated into the system to automate multiple sample handling steps with minimal user intervention. The device enables quantitative protein analysis with picomolar sensitivity, as demonstrated using human samples spiked with interleukin-6 and C-reactive protein. The affinity-based assays are highly specific to the target without cross-reactivity. Therefore, we envision the reported device offering ultrasensitive and field-deployable immunoassay tests for timely and accurate clinical diagnosis.
The need for clean, renewable energy has driven the expansion of renewable energy generators, such as wind and solar. However, to achieve a robust and responsive electrical grid based on such inherently intermittent renewable energy sources, grid-scale energy storage is essential. The unmet need for this critical component has motivated extensive grid-scale battery research, especially exploring chemistries “beyond Li-ion”. Among others, molten sodium (Na) batteries, which date back to the 1960s with Na-S, have seen a strong revival, owing mostly to raw material abundance and the excellent electrochemical properties of Na metal. Recently, many groups have demonstrated important advances in battery chemistries, electrolytes, and interfaces to lower material and operating costs, enhance cyclability, and understand key mechanisms that drive failure in molten Na batteries. For widespread implementation of molten Na batteries, though, further optimization, cost reduction, and mechanistic insight is necessary. In this light, this work provides a brief history of mature molten Na technologies, a comprehensive review of recent progress, and explores possibilities for future advancements.
Metal-organic frameworks (MOFs) are a class of porous, crystalline materials that have been systematically developed for a broad range of applications. Incorporation of two or more metals into a single crystalline phase to generate heterometallic MOFs has been shown to lead to synergistic effects, in which the whole is oftentimes greater than the sum of its parts. Because geometric proximity is typically required for metals to function cooperatively, deciphering and controlling metal distributions in heterometallic MOFs is crucial to establish structure-function relationships. However, determination of short- and long-range metal distributions is nontrivial and requires the use of specialized characterization techniques. Advancements in the characterization of metal distributions and interactions at these length scales is key to rapid advancement and rational design of functional heterometallic MOFs. This perspective summarizes the state-of-the-art in the characterization of heterometallic MOFs, with a focus on techniques that allow metal distributions to be better understood. Using complementary analyses, in conjunction with computational methods, is critical as this field moves toward increasingly complex, multifunctional systems.
Rothchild, Eric; Asta, Mark D.; Chrzan, Daryl C.; Kuner, Matthew C.
The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.
We report a comparative study of three rectifying gate metals, W, Pd, and Pt/Au, on ultrawide bandgap Al0.86Ga0.14N barrier/Al0.7Ga0.3N channel high electron mobility transistors for use in extreme temperatures. The transistors were electrically characterized from 30 to 600 °C in air. Of the three gate metals, the Pt/Au stack exhibited the smallest change in threshold voltage (0.15 V, or 9% change between the 30 and 600 °C values, and a maximum change of 42%), the highest on/off current ratio (1.5 × 106) at 600 °C, and a modest forward gate leakage current (0.39 mA/mm for a 3 V gate bias) at 600 °C. These favorable results showcase AlGaN channel high electron mobility transistors' ability to operate in extreme temperature environments.
A computational model of aluminum melting is proposed which captures both the thermal fluid-solid phase transition and the mechanical effects of oxidation. The model hybridizes ideas from smoothed particle hydrodynamics and bonded particle models to simulate both hydrodynamic flows and solid elasticity. Oxidation is represented by dynamically adding and deleting spring-like bonds between surface fluid particles to represent the formation and rupture of the oxide skin. Various complex systems are simulated to demonstrate the adaptability of the method and to illustrate the significant impact of skin properties on material flow. As a result, initial comparison to experiments of a melting aluminum cantilever highlights that the computational model can reproduce key qualitative features of aluminum relocation.
Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against 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 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.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.
A generalized closed-form equation for the shaded collector fraction in solar arrays on rolling or undulating terrain is provided for single-axis tracking and fixed-tilt systems. The equation accounts for different rotation angles between the shaded and shading trackers, cross-axis slope between the two trackers, and offset between the collector plane and axis of rotation. The validity of the equation is demonstrated through comparison with numerical ray-tracing simulations and remaining minor sources of error are quantified. Additionally, a simple procedure to determine backtracking rotations for each row in an array installed on the rolling terrain (varying in the direction perpendicular to the tracker axes) is provided. The backtracking equation accounts for a desired shaded fraction (including complete shade avoidance) as well as an axis-collector offset. Test cases are provided to facilitate implementation of these equations.
Recent experimental studies suggest the use of spatially extended laser beam profiles as a strategy to control the melt pool during laser powder bed fusion (LPBF) additive manufacturing. However, linkages connecting laser beam profiles to thermal fields and resultant microstructures have not been established. Herein, we employ a coupled thermal transport-Monte Carlo model to predict the evolution of temperature fields and grain microstructures during LPBF using Gaussian, ring, and Bessel beam profiles. Simulation results reveal that the ring-shaped beam yields lower temperatures compared with the Gaussian beam. Owing to the small melt pool size when using the Bessel beam, the grains are smaller in size and more equiaxed compared to those using the Gaussian and ring beams. Our approach provides future avenues to predict the impact of laser beam shaping on microstructure development during LPBF.
Peridigm is a meshfree peridynamics code written in C++ for use on large-scale parallel computers. It was originally developed at Sandia National Laboratories and is currently managed as an open-source, community driven software project. Its primary features include bond-based, state-based, and non-ordinary state-based constitutive models, bond failure laws, contact, and support for explicit and implicit time integration. To date, Peridigm has been used primarily by methods developers focused on solid mechanics and material failure. Peridigm utilizes foundational software components from Sandia’s Trilinos project and was designed for extensibility. This paper provides an overview of the solution methods implemented in Peridigm, a discussion of its software infrastructure, and demonstrates the use of Peridigm for the solution of several example problems.
Magnetized Liner Inertial Fusion experiments have been performed at the Z facility at Sandia National Laboratories. These experiments use deuterium fuel, which produces 2.45 MeV neutrons on reaching thermonuclear conditions. To study the spatial structure of neutron production, the one-dimensional imager of neutrons diagnostic was fielded to record axial resolved neutron images. In this diagnostic, neutrons passing through a rolled edge aperture form an image on a CR-39-based solid state nuclear track detector. Here, we present a modified generalized expectation-maximization algorithm to reconstruct an axial neutron emission profile of the stagnated fusion plasma. We validate the approach by comparing the reconstructed neutron emission profile to an x-ray emission profile provided by a time-integrated pinhole camera.
Journal of Vacuum Science and Technology A: Vacuum, Surfaces and Films
Jaszewski, Samantha T.; Fields, Shelby S.; Chung, Ching C.; Jones, Jacob L.; Orson, Keithen G.; Reinke, Petra; Ihlefeld, Jon F.
The impact of the high-power impulse magnetron sputtering (HiPIMS) pulse width on the crystallization, microstructure, and ferroelectric properties of undoped HfO2 films is investigated. HfO2 films were sputtered from a hafnium metal target in an Ar/O2 atmosphere, varying the instantaneous power density by changing the HiPIMS pulse width with fixed time-averaged power and pulse frequency. The pulse width is shown to affect the ion-to-neutral ratio in the depositing species with the shortest pulse durations leading to the highest ion fraction. In situ x-ray diffraction measurements during crystallization demonstrate that the HiPIMS pulse width impacts nucleation and phase formation, with an intermediate pulse width of 110 μs stabilizing the ferroelectric phase over the widest temperature range. Although the pulse width impacts the grain size with the lowest pulse width resulting in the largest grain size, the grain size does not strongly correlate with the phase content or ferroelectric behavior in these films. These results suggest that precise control over the energetics of the depositing species may be beneficial for forming the ferroelectric phase in this material.
The complex nature of manufacturing processes stipulates electrodes to possess high variability with increased heterogeneity during production. X-ray computed tomography imaging has proved to be critical in visualizing the complicated stochastic particle distribution of as-manufactured electrodes in lithium-ion batteries. However, accurate prediction of their electrochemical performance necessitates precise evaluation of kinetic and transport properties from real electrodes. Image segmentation that characterizes voxels to particle/pore phase is often meticulous and fraught with subjectivity owing to a myriad of unconstrained choices and filter algorithms. We utilize a Bayesian convolutional neural network to tackle segmentation subjectivity and quantify its pertinent uncertainties. Otsu inter-variance and Blind/Referenceless Imaging Spatial Quality Evaluator are used to assess the relative image quality of grayscale tomograms, thus evaluating the uncertainty in the derived microstructural attributes. We analyze how image uncertainty is correlated with the uncertainties and magnitude of kinetic and transport properties of an electrode, further identifying pathways of uncertainty propagation within microstructural attributes. The coupled effect of spatial heterogeneity and microstructural anisotropy on the uncertainty quantification of transport parameters is also understood. This work demonstrates a novel methodology to extract microstructural descriptors from real electrode images through quantification of associated uncertainties and discerning the relative strength of their propagation, thus facilitating feedback to manufacturing processes from accurate image based electrochemical simulations.
The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. The method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.
Achieving commercially acceptable Zn-MnO2 rechargeable batteries depends on the reversibility of active zinc and manganese materials, and avoiding side reactions during the second electron reaction of MnO2. Typically, liquid electrolytes such as potassium hydroxide (KOH) are used for Zn-MnO2 rechargeable batteries. However, it is known that using liquid electrolytes causes the formation of electrochemically inactive materials, such as precipitation Mn3O4 or ZnMn2O4 resulting from the uncontrollable reaction of Mn3+ dissolved species with zincate ions. In this paper, hydrogel electrolytes are tested for MnO2 electrodes undergoing two-electron cycling. Improved cell safety is achieved because the hydrogel electrolyte is non-spillable, according to standards from the US Department of Transportation (DOT). The cycling of “half cells” with advanced-formulation MnO2 cathodes paired with commercial NiOOH electrodes is tested with hydrogel and a normal electrolyte, to detect changes to the zincate crossover and reaction from anode to cathode. These half cells achieved ≥700 cycles with 99% coulombic efficiency and 63% energy efficiency at C/3 rates based on the second electron capacity of MnO2. Other cycling tests with “full cells” of Zn anodes with the same MnO2 cathodes achieved ~300 cycles until reaching 50% capacity fade, a comparable performance to cells using liquid electrolyte. Electrodes dissected after cycling showed that the liquid electrolyte allowed Cu ions to migrate more than the hydrogel electrolyte. However, measurements of the Cu diffusion coefficient showed no difference between liquid and gel electrolytes; thus, it was hypothesized that the gel electrolytes reduced the occurrence of Cu short circuits by either (a) reducing electrode physical contact to the separator or (b) reducing electro-convective electrolyte transport that may be as important as diffusive transport.
Current nuclear facility emergency planning zones (EPZs) are based on outdated distance-based criteria, predating comprehensive dose and risk-informed frameworks. Recent advancements in simulation tools have permitted the development of site-specific, dose, and risk-based consequence-driven assessment frameworks. This study investigated the computation of advanced reactor (AR) EPZs using two atmospheric dispersion models: a straight-line Gaussian plume model (GPM) and a semi-Lagrangian Particle in Cell (PIC). Two case studies were conducted: (1) benchmarking the NRC SOARCA study for the Peach Bottom Nuclear Generating Station and (2) analyzing an advanced INL Heat Pipe Design A microreactor's end-of-cycle inventory. The dose criteria for both cases were 10 mSv at mean weather conditions and 50 mSv at 95th percentile weather conditions at 96 h post-release. Results demonstrated that GPM and PIC estimated similar mean peak dose levels for large boiling water reactors in the farfield case, placing EPZ limits beyond current regulations. For ARs with source terms remaining in the nearfield, PIC modeling without specific nearfield considerations could result in excessively high doses and inaccurate EPZ designations. PIC dispersion demonstrated an order of magnitude higher estimate of nearfield inhalation dose contribution when compared to GPM results. Both models significantly reduced EPZ sizing within the nearfield. Thus, reductions in the AR source term may eliminate the need for a separate EPZ.
Heating of the surficial layer of the atmosphere often generates convective vortices, known as “dust devils” when they entrain visible debris. Convective vortices are common on both Earth and Mars, where they affect the climate via dust loading, contribute to wind erosion, impact the efficiency of photovoltaic systems, and potentially result in injury and property damage. However, long-duration terrestrial convective vortex activity records are rare. We have developed a high-precision and high-recall method to extract convective vortex signatures from infrasound microbarometer data streams. The techniques utilizes a wavelet-based detector to capture potential events and then a template matching system to extract the duration of the vortex. Since permanent and temporary infrasound sensors networks are present throughout the globe (many with open data), our method unlocks a vast new convective vortex dataset without requiring the deployment of specialized instrumentation. SIGNIFICANCE STATEMENT: Convective vortices, or “dust devils,” contribute to regional dust loading in Earth’s atmosphere. However, long-duration convective vortex activity records are rare. We came up with a way to autonomously detect the pressure signatures left by convective vortices striking low-frequency sound, or “infrasound,” sensors. Since permanent infrasound stations have been active for decades, our method has the potential to add ordersof-magnitude more events than previously catalogued.
Seleson, Pablo; Pasetto, Marco; John, Yohan; Trageser, Jeremy; Reeve, Samuel T.
PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. Due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.
Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. Though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.
Laccases from white-rot fungi catalyze lignin depolymerization, a critical first step to upgrading lignin to valuable biodiesel fuels and chemicals. In this study, a wildtype laccase from the basidiomycete Fomitiporia mediterranea (Fom_lac) and a variant engineered to have a carbohydrate-binding module (Fom_CBM) were studied for their ability to catalyze cleavage of β-O-4′ ether and C–C bonds in phenolic and non-phenolic lignin dimers using a nanostructure-initiator mass spectrometry-based assay. Fom_lac and Fom_CBM catalyze β-O-4′ ether and C–C bond breaking, with higher activity under acidic conditions (pH < 6). The potential of Fom_lac and Fom_CBM to enhance saccharification yields from untreated and ionic liquid pretreated pine was also investigated. Adding Fom_CBM to mixtures of cellulases and hemicellulases improved sugar yields by 140% on untreated pine and 32% on cholinium lysinate pretreated pine when compared to the inclusion of Fom_lac to the same mixtures. Adding either Fom_lac or Fom_CBM to mixtures of cellulases and hemicellulases effectively accelerates enzymatic hydrolysis, demonstrating its potential applications for lignocellulose valorization. We postulate that additional increases in sugar yields for the Fom_CBM enzyme mixtures were due to Fom_CBM being brought more proximal to lignin through binding to either cellulose or lignin itself.
During hypersonic flight, compressional and viscous heating of the air can form a plasma layer which encases the aircraft. If the boundary layer becomes turbulent, then the electron density fluctuations can effect a parasitic modulation in microwave signals transmitted through the plasma. We developed an approach for studying the interaction of microwave signals with a turbulent, hypersonic plasma layer. The approach affords a great deal of flexibility in both the plasma layer model and the antenna configuration. We then analyzed a situation in which microwaves, transmitted from a rectangular aperture antenna, propagate through a turbulent plasma layer to a distant receiver. We characterized the first- and second-order statistics of the computed parasitic modulation and quantified the depolarization of the signal. The amplitude fluctuations are lognormally distributed at low frequencies and Rice-distributed at high frequencies. Fluctuations in the copolarized phase and amplitude of the far-field signal are strongly anticorrelated. We used a multioutput Gaussian process (MOGP) to model these quantities. The efficacy of the MOGP model is demonstrated by recovering the time evolution of the copolarized phase given the copolarized amplitude and occasional measurements of the phase.
The size of a pressure transducer is known to affect the accuracy of measurements of wall-pressure fluctuations beneath a turbulent boundary layer because of spatial averaging over the sensing area of the transducer. In this paper, the effect of finite transducer size is investigated by applying spatial averaging or wavenumber filters to a database of hypersonic wall pressure generated from a direct numerical simulation (DNS) that simulates the turbulent portion of the boundary layer over a sharp 7° half-angle cone at nominally Mach 8. A good comparison between the DNS and the experiment in the Sandia Hypersonic Wind Tunnel at Mach 8 is achieved after spatial averaging is applied to the DNS data over an area similar to the sensing area of the transducer. The study shows that a finite sensor size similar to that of the PCB132 transducer can cause significant attenuation in the root-mean-square and power spectral density (PSD) of wall-pressure fluctuations, and the attenuation effect is identical between cone and flat plate configurations at the same friction Reynolds number. The Corcos theory is found to successfully compensate for the attenuated highfrequency components of the wall-pressure PSD.
During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na2O-SiO2(NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na2O-SiO2(NS15) composition, the 25% Na2O-SiO2(NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na2O-SiO2 compositions. Coordination analysis of the Na+ ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na+ ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.
Ground heat flux (G0) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G0 data (measured or modeled) for snow-free period as a reference. When observed G0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.
This document is a guide to setting the system and processing configuration for the Geophysical Monitoring System (GMS) Interactive Analysis (IAN) application.
This document lays out a set of near-future investigations in salt, the third phase of BATS (BATS 3). This phase is planned to answer the few remaining issues from the first two phases of BATS (BATS 1 and BATS 2), and to prepare for a subsequent large-scale demonstration phase. The BATS experiments are the first part of a larger plan to conduct field experiments to answer specific technical questions, improve the technical basis for disposal of heat-generating radioactive waste in salt (Stauffer et al., 2015; SNL et al., 2020), and demonstrate readiness for disposal of radioactive waste in salt, including large, hot waste packages.