Discharge of sodium coolant into containment from a sodium-cooled fast reactor vessel can occur in the event of a pipe leak or break. In this situation, some of the liquid sodium droplets discharged from the coolant system will react with oxygen in the air before reaching the containment. This phase of the event is normally termed the sodium spray fire phase. Unreacted sodium droplets pool on the containment floor where continued reaction with containment atmospheric oxygen occurs. This phase of the event is normally termed the sodium pool fire phase. Both phases of these sodium-oxygen reactions (or fires) are important to model because of the heat addition and aerosol generation that occur. Any fission products trapped in the sodium coolant may also be released during this progression of events, which if released from containment could pose a health risk to workers and the public. The paper describes progress of an international collaborative research in the area of the sodium fire modeling in the sodium-cooled fast reactors between the United States and Japan under the framework of the Civil Nuclear Energy Research and Development Working Group. In this collaboration between Sandia National Laboratories and Japan Atomic Energy Agency, the validation basis for and modeling capabilities of sodium spray and pool fires in MELCOR of Sandia National Laboratories and SPHINCS of Japan Atomic Energy Agency are being enhanced. This study documents MELCOR and SPHINCS sodium pool fire model validation exercises against the JAEA's sodium pool fire experiments, F7-1 and F7-2. The proposed enhancement of the sodium pool fire models in MELCOR through addition of thermal hydraulic and sodium spreading models that enable a better representation of experimental results is also described.
The domain wall-magnetic tunnel junction (DW-MTJ) is a versatile device that can simultaneously store data and perform computations. These three-terminal devices are promising for digital logic due to their nonvolatility, low-energy operation, and radiation hardness. Here, we augment the DW-MTJ logic gate with voltage-controlled magnetic anisotropy (VCMA) to improve the reliability of logical concatenation in the presence of realistic process variations. VCMA creates potential wells that allow for reliable and repeatable localization of domain walls (DWs). The DW-MTJ logic gate supports different fanouts, allowing for multiple inputs and outputs for a single device without affecting the area. We simulate a systolic array of DW-MTJ multiply-accumulate (MAC) units with 4-bit and 8-bit precision, which uses the nonvolatility of DW-MTJ logic gates to enable fine-grained pipelining and high parallelism. The DW-MTJ systolic array provides comparable throughput and efficiency to state-of-the-art CMOS systolic arrays while being radiation-hard. These results improve the feasibility of using DW-based processors, especially for extreme-environment applications such as space.
Iodide redox reactions in molten NaI/AlCl3 are shown to generate surface-blocking films, which may limit the useful cycling rates and energy densities of molten sodium batteries below 150 °C. An experimental investigation of electrode interfacial stability at 110 °C reveals the source of the reaction rate limitations. Electrochemical experiments in a 3-electrode configuration confirm an increase of resistance on the electrode surface after oxidation or reduction current is passed. Using chronopotentiometry, chronoamperometry, cyclic voltammetry, and electrochemical impedance spectroscopy, the film formation is shown to depend on the electrode material (W, Mo, Ta, or glassy carbon), as well as the Lewis acidity and molar ratio of I−/I3− in the molten salt electrolytes. These factors impact the amount of charge that can be passed at a given current density prior to developing excessive overpotential due to film formation that blocks the electrode surface. The results presented here guide the design and use of iodide-based molten salt electrolytes and electrode materials for grid scale battery applications.
New technology for electric vehicles (EVs) must meet the requirements of higher energy usage, lower costs, and more sustainable source materials. One promising material for EV power system components is iron nitride (IN) soft magnetic composites (SMCs) because of their competitive magnetic properties and high abundance of the source materials. As part of an ongoing program at Sandia National Laboratories, this project focused on using computer modeling to optimize the prototyping process for an iron nitride SMC toroidal inductor to reach a target inductance of 600 μH. Four inductors with different combinations of wiring (26 AWG and 20 AWG) and vol% loading of iron nitride (65 vol% and 50 vol%) were fabricated at Cal Poly and characterized using an LCR meter. These inductors were also modeled using COMSOL Multiphysics™ with the Magnetic Fields module. The inductance data from the experiment and the model show that the 65 vol% IN prototypes and models agree with about 8% difference, while the 50 vol% IN samples show about a 9% difference between the prototype and the model. These results suggest that the model can predict inductance with both accuracy and precision with low confidence for the given sample size of four. An additional parameter of AC resistance is studied but the AC resistance results from the inductors and from the model generally do not agree closely, suggesting that the current model used in the project does not fully capture the mechanisms behind AC resistance of the inductor. With the focus of the project on inductance, the percent difference results of less than 9% across the four inductors that were tested increases confidence in the model’s predictive capabilities for inductance only. Using the inductance results from both the model and experiment, the final suggested inductor design is a 65 vol% core with 150 windings of 20 AWG wire that is 8 cm across and 1.5 cm tall to reach the inductance goal of 600 μH based on analysis using the optimized COMSOLTM model.
Yttrium iron garnet (Y3Fe5O12; YIG) has a unique combination of low magnetic damping, high spin-wave conductivity, and insulating properties that make it a highly attractive material for a variety of applications in the fields of magnetics and spintronics. While the room-temperature magnetization dynamics of YIG have been extensively studied, there are limited reports correlating the low-temperature magnetization dynamics to the material structure or growth method. Here, in this study, we investigate liquid phase epitaxy grown YIG films and their magnetization dynamics at temperatures down to 10 K. We show there is a negligible increase in the ferromagnetic resonance linewidth down to 10 K, which is unique when compared with YIG films grown by other deposition methods. From the broadband ferromagnetic resonance measurements, polarized neutron reflectivity, and scanning transmission electron microscopy, we conclude that these liquid phase epitaxy grown films have negligible rare-earth impurities present, specifically the suppression of Gd diffusion from the Gd3Ga5O12 (GGG) substrate into the Y3Fe5O12 film, and therefore negligible magnetic losses attributed to the slow-relaxation mechanism. Overall, liquid phase epitaxy YIG films have a YIG/GGG interface that is five times sharper and have ten times lower ferromagnetic resonance linewidths below 50 K than comparable YIG films by other deposition methods. Thus, liquid phase epitaxy grown YIG films are ideal for low-temperature experiments/applications that require low magnetic losses, such as quantum transduction and manipulation via magnon coupling.
IEEE Transactions on Terahertz Science and Technology
Curwen, Christopher A.; Kawamura, Jonathan H.; Hayton, Darren J.; Addamane, Sadhvikas J.; Reno, John L.; Williams, Benjamin S.; Karasik, Boris S.
We report high-resolution frequency study and phase locking have been performed on a terahertz (THz) quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL) operating around 2.5 THz. A subharmonic diode mixer is used to down convert the THz signal to a 100 MHz intermediate frequency that is phase locked to a stable 100 MHz microwave reference. Between 90% and 95% of the QC-VECSEL signal is locked within 2 Hz of the multiplied RF reference, and amplitude fluctuations on the order of 1%–10% are observed, depending on the bias point of the QC-VECSEL. The bandwidth of the locking loop is ~1 MHz. Many noise peaks in the IF signal are observed, likely corresponding to mechanical resonances in the 10 Hz–10 kHz. These peaks are generally -30 to -60 dB below the main tone and are below the phase noise level of the multiplied RF reference that ultimately limits the phase noise of the locked QC-VECSEL.
Zheng, Jeffrey X.; Fiagbenu, Merrilyn M.A.; Esteves, Giovanni; Musavigharavi, Pariasadat; Jariwala, Deep; Stach, Eric A.; Olsson, Roy H.
Ferroelectric Al1−xScxN has raised much interest in recent years due to its unique ferroelectric properties and complementary metal oxide semiconductor back-end-of-line compatible processing temperatures. Potential applications in embedded nonvolatile memory, however, require ferroelectric materials to switch at relatively low voltages. One approach to achieving a lower switching voltage is to significantly reduce the Al1−xScxN thickness. In this work, ferroelectric behavior in 5-27 nm films of sputter deposited Al0.72Sc0.28N has been studied. We find that the 10 kHz normalized coercive field increases from 4.4 to 7.3 MV/cm when reducing the film thickness from 27.1 to 5.4 nm, while over the same thickness range, the characteristic breakdown field of a 12.5 μm radius capacitor increases from 8.3 to 12.1 MV/cm. The 5.4 nm film demonstrates ferroelectric switching at 5.5 V when excited with a 500 ns pulse and a switching speed of 60 ns.
Sulfur-deficient polycrystalline two-dimensional (2D) molybdenum disulfide (MoS2) memtransistors exhibit gate-tunable memristive switching to implement emerging memory operations and neuromorphic computing paradigms. Grain boundaries and sulfur vacancies are critical for memristive switching; however, the underlying physical mechanisms are not fully understood. Furthermore, the adsorption of water and gaseous species strongly perturbs electronic transport in monolayer MoS2, and little work has been done to explore the influence of surface interactions on defect-related kinetics that produces memristive switching. Here, we study the switching kinetics of back-gated MoS2 memtransistors using current transient measurements in a controlled atmosphere chamber. We observe that adsorbed water molecules lead to suppression of the electronic trap-filling processes concomitant with the resistive switching process, resulting in altered kinetics of the resistive switching. Additionally, using the transient response from “bunched” drain voltage pulse trains performed as a function of temperature, we extract the energy of the affected trap state and find that it places the trap roughly midgap [ E T = E C - 0.7 ( ± 0.4 ) eV]. Our results highlight the importance of controlling for surface interactions that may affect switching kinetics in 2D memtransistors, synaptic transistors, and related memory devices.
Cerjan, Alexander; Loring, Terry A.; Cheng, Wenting; Ying Chen, Ssu; Prodan, Camelia; Prodan, Emil
Topological metals are conducting materials with gapless band structures and nontrivial edge-localized resonances. Their discovery has proven elusive because traditional topological classification methods require band gaps to define topological robustness. Inspired by recent theoretical developments that leverage techniques from the field of C*-algebras to identify topological metals, here, we directly observe topological phenomena in gapless acoustic crystals and realize a general experimental technique to demonstrate their topology. Specifically, we not only observe robust boundary-localized states in a topological acoustic metal, but also re-interpret a composite operator—mathematically derived from the K-theory of the problem—as a new Hamiltonian whose physical implementation allows us to directly observe a topological spectral flow and measure the topological invariants. Our observations and experimental protocols may offer insights for discovering topological behaviour across a wide array of artificial and natural materials that lack bulk band gaps.
This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSMs) from Orano (formerly Transnuclear Inc.) and provides a summary of the surface sampling activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. One goal of this testing is to collect defensible and detailed dust deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion cracking (CISCC) research. Another goal is to provide data for the validation of computational fluid dynamics (CFD) based deposition modeling. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the decay heat and thermal hydraulic environment. Test and heater design is supported by detailed CFD modeling. Instrumentation throughout the canister, storage module, and environment will provide extensive information about the thermal-hydraulic behavior of horizontal dry cask storage systems. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer detailed quantification and composition of the deposited particulates from a realistic storage environment. Discussions of a potential host site for the Canister Deposition Field Demonstration (CDFD) are ongoing. Until a host site is chosen, testing of key CDFD hardware components including the heater assemblies, power skid, and remote data acquisition system will continue. Functional testing of the finalized heater assemblies and test apparatus started this fiscal year. These initial heater tests have shown the assemblies are performing within design specifications. Staged surface sampling of a mockup of a canister outside the AHSM on a transfer skid was also performed. Refinements to the sampling procedures and techniques were captured from observation of these activities and lessons-learned debriefs. These updated sampling procedures and techniques are planned to be tested again in the field using the mockup in order to assure personnel are using the most accurate and repeatable methods possible prior to deployment for actual CDFD testing.
This report describes the results of preliminary testing of aerosol monitoring equipment that will be used to continuously monitor the aerosol source term for the multi-year Canister Deposition Field Demonstration (CDFD). These data are required inputs for the development and validation of models for the deposition of dust and potentially corrosive salts on the surface of spent nuclear fuel (SNF) dry storage canisters. Surface salt loads correlate with the extent of corrosion damage on a metal surface, and potentially to the likelihood and timing of initiation of stress corrosion cracks. Aerosols will be monitored at the CDFD site using three instruments. A Dekati® ELPI+ cascade impactor will be used for real-time monitoring of aerosol particle sizes. It will also collect dust in 14 size bins on impactor targets that can be chemically analyzed to determine the soluble salts present as a function of particle size. However, this instrument can only measure dried aerosols, with a diameter of <10 µm. The second instrument is a Topas laser particle size spectrometer, which provides real-time monitoring of aerosol particle sizes up to ~40 µm in size. It monitors both the ambient (potentially deliquesced) aerosol particle size distributions required for the dust deposition models and the distributions of the equivalent dried particles, allowing correlation with the Dekati® data. However, it does not discriminate between inert dust particles and salt aerosols, and it does not retain samples of the different particle sizes for later analysis. The third instrument that will monitor aerosols at the CDFD site is a Clean Air Status and Trends Network (CASTNET) tower, which uses a multiple canister system to collect weekly samples for analysis to total suspended aerosol particle compositions and atmospheric gas concentrations. This status report describes work in FY23 to develop the capabilities for using these tools. In two training exercises, the cascade impactor and laser particle sizer were deployed in two different testing environments, one indoor and one outdoor. For the cascade impactor, the tests provided opportunities for the operators to familiarize themselves with impactor substrate preparation, and post-test sample removal and analysis. For the laser particle sizer, the tests were used to evaluate different instrument parameters, to determine the most appropriate settings for capturing transient events. Data and samples were collected for weeks to months for each test, and the results are presented here. In addition to the preliminary testing, contracts were developed with WSP Analytical Labs for sample preparation and analysis of the cascade impactor samples. The impactor tower from outdoor test was delivered to WSP and used to train the staff there in disassembly, sample extraction, sample analysis, and tower reassembly with new target substrates. These are tasks that WSP will be performing routinely for the CDFD project. The CASTNET system cannot be purchased or tested until an actual site has been selected for the CDFD test. Work for this FY has been restricted to preparation of contracts for purchasing the CASTNET tower, and for sample analysis, once the tower is in operation.
Solar-powered photochemical water splitting using suspensions of photocatalyst nanoparticles is an attractive route for economical production of green hydrogen. SrTiO3-based photocatalysts have been intensely investigated due to their stability and recently demonstrated near-100% external quantum yield (EQY) for water splitting using wavelengths below 360 nm. To extend the optical absorption into the visible, SrTiO3nanoparticles have been doped with various transition metals. Here we demonstrate that doping SrTiO3nanoparticles with 1% Rh introduces midgap acceptor states which reduce the free electron concentration by 5 orders of magnitude, dramatically reducing built-in potentials which could otherwise separate electron-hole (e-h) pairs. Rhodium states also function as recombination centers, reducing the photocarrier lifetime by nearly 2 orders of magnitude and the maximum achievable EQY to 10%. Furthermore, the absence of built-in electric fields within Rh-doped SrTiO3nanoparticles suggests that modest e-h separation can be achieved by exploiting a difference in mobility between electrons and holes.
Here this paper introduces a publicly available PyTorch-ABAQUS deep-learning framework of a family of plasticity models where the yield surface is implicitly represented by a scalar-valued function. In particular, our focus is to introduce a practical framework that can be deployed for engineering analysis that employs a user-defined material subroutine (UMAT/VUMAT) for ABAQUS, which is written in FORTRAN. To accomplish this task while leveraging the back-propagation learning algorithm to speed up the neural-network training, we introduce an interface code where the weights and biases of the trained neural networks obtained via the PyTorch library can be automatically converted into a generic FORTRAN code that can be a part of the UMAT/VUMAT algorithm. To enable third-party validation, we purposely make all the data sets, source code used to train the neural-network-based constitutive models, and the trained models available in a public repository. Furthermore, the practicality of the workflow is then further tested on a dataset for anisotropic yield function to showcase the extensibility of the proposed framework. A number of representative numerical experiments are used to examine the accuracy, robustness and reproducibility of the results generated by the neural network models.
A significant amount of uncertainty exists regarding potential human exposure to laboratory biomaterials and organisms in Biosafety Level 2 (BSL-2) research laboratories. Computational fluid dynamics (CFD) modeling is proposed as a way to better understand potential impacts of different combinations of biomaterials, laboratory manipulations, and exposure routes on risks to laboratory workers. Here, in this study, we use CFD models to simulate airborne concentrations of contaminants in an actual BSL-2 laboratory under different configurations. Results show that ventilation configuration, sampling location, and contaminant source location can significantly impact airborne concentrations and exposures. Depending on the source location and airflow patterns, the transient and time-integrated concentrations varied by several orders of magnitude. Contaminant plumes from sources located near a return vent (or exhaust like a fume hood or ventilated biosafety cabinet) are likely to be more contained than sources that are further from the exhaust. Having a direct flow between the source and the exhaust (through-flow condition) may reduce potential exposures to individuals outside the air flow path. Designing a BSL-2 room with ventilation and airflow patterns that maximize through-flow conditions to the return/exhaust vents and minimize dispersion and mixing throughout the room is, therefore, recommended. CFD simulations can also be used to assist in characterizing the impacts of supply and return vent locations, room layout, and source locations on spatial and temporal contaminant concentrations. In addition, proper placement of particle sensors can also be informed by CFD simulations to provide additional characterization and monitoring of potential exposures in BSL-2 facilities.
This report is the revised (Revision 10) Task F specification for DECOVALEX-2023. Task F is a comparison of the models and methods used in deep geologic repository performance assessment. The task proposes to develop a reference case for a mined repository in a fractured crystalline host rock (Task F1) and a reference case for a mined repository in a salt formation (Task F2). Teams may choose to participate in the comparison for either or both reference cases. For each reference case, a common set of conceptual models and parameters describing features, events, and processes that impact performance will be given, and teams will be responsible for determining how best to implement and couple the models. The comparison will be conducted in stages, beginning with a comparison of key outputs of individual process models, followed by a comparison of a single deterministic simulation of the full reference case, and moving on to uncertainty propagation and uncertainty and sensitivity analysis. This report provides background information, a summary of the proposed reference cases, and a staged plan for the analysis.
Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.
Pb-Zr-Ti-O (PZT) perovskites span a large solid-solution range and have found widespread use due to their piezoelectric and ferroelectric properties that also span a large range. Crystal structure analysis via Rietveld refinement facilitates materials analysis via the extraction of the structural parameters. These parameters, often obtained as a function of an additional dimension (e.g., pressure), can help to diagnose materials response within a use environment. Often referred to as in-situ studies, these experiments provide an abundance of data. Viewing structural changes due to applied pressure conditions can give much-needed insight into materials performance. However, challenges exist for viewing/presenting results when the details are inherently three-dimensional (3D) in nature. For PZT perovskites, the use of polyhedra (e.g., Zr/Ti-O6 octahedra) to view bonding/connectivity is beneficial; however, the visualization of the octahedra behavior with pressure dependence is less easily demonstrated due to the complexity of the added pressure dimension. We present a more intuitive visualization by projecting structural data into virtual reality (VR). We employ previously published structural data for Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3 as an exemplar for VR visualization of the PZT R3c crystal structure between ambient and 0.62 GPa pressure. This is accomplished via our in-house CAD2VR™ software platform and the new CrystalVR plugin. The use of the VR environment enables a more intuitive viewing experience, while enabling on-the-fly evaluation of crystal data, to form a detailed and comprehensive understanding of in-situ datasets. Discussion of methodology and tools for viewing are given, along with how recording results in video form can enable the viewing experience.
Chong, Lina; Gao, Guoping; Wen, Jianguo; Li, Haixia; Xu, Haiping; Green, Zach; Sugar, Joshua D.; Kropf, A.J.; Xu, Wenqian; Lin, Xiao M.; Xu, Hui; Wang, Lin W.; Di Liu, Jia
Discovery of earth-abundant electrocatalysts to replace iridium for the oxygen evolution reaction (OER) in a proton exchange membrane water electrolyzer (PEMWE) represents a critical step in reducing the cost for green hydrogen production. We report a nanofibrous cobalt spinel catalyst codoped with lanthanum (La) and manganese (Mn) prepared from a zeolitic imidazolate framework embedded in electrospun polymer fiber. The catalyst demonstrated a low overpotential of 353 millivolts at 10 milliamperes per square centimeter and a low degradation for OER over 360 hours in acidic electrolyte. A PEMWE containing this catalyst at the anode demonstrated a current density of 2000 milliamperes per square centimeter at 2.47 volts (Nafion 115 membrane) or 4000 milliamperes per square centimeter at 3.00 volts (Nafion 212 membrane) and low degradation in an accelerated stress test.
Glusa, Christian; D'Elia, Marta; Alla, Alessandro; Oliveira, Hugo
In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.
The efficient utilization of lignin, the direct source of renewable aromatics, into value-added renewable chemicals is an important step towards sustainable biorefinery practices. Nevertheless, owing to the random heterogeneous structure and limited solubility, lignin utilization has been primarily limited to burning for energy. The catalytic depolymerization of lignin has been proposed and demonstrated as a viable route to sustainable biorefinery, however, low yields and poor selectivity of products, high char formation, and limited to no recycling of transition-metal-based catalyst involved in lignin depolymerization demands attention to enable practical-scale lignocellulosic biorefineries. In this study, we demonstrate the catalytic depolymerization of ionic liquid-based biorefinery poplar lignin into guaiacols over a reusable zirconium phosphate supported palladium catalyst. The essence of the study lies in the high conversion (>80 %), minimum char formation (7–16 %), high yields of guaiacols (up to 200 mg / g of lignin), and catalyst reusability. Both solid residue, liquid stream, and gaseous products were thoroughly characterized using ICP-OES, PXRD, CHN analysis, GC-MS, GPC, and 2D NMR to understand the hydrogenolysis pathway.
Oriented attachment (OA) of nanoparticles is an important crystal growth pathway in the synthesis of hierarchical structures. Although a significant understanding of OA has been made, the effect of atomistic misalignment and the roles of solvent/particle and particle/particle interactions on the structure-energy relationship during an OA remain elusive. In this study, we perform molecular dynamics simulations to calculate the potential of mean force (PMF) profile for gibbsite particle translation on a gibbsite slab with 1 or 2 intervening water layers (1W or 2W). The structures of the gibbsite surfaces and the confined water are analyzed to determine how the number and type of hydrogen bonds (H-bonds) influence the free energy profile during the translation. The PMF profile exhibits a periodicity of length 5.078 Å, consistent with the gibbsite unit cell size along the translation direction. The changes in the surface-water and water-water hydrogen bond network and water and surface OH groups’ orientations during the translation are strongly coupled with the changes in the PMF profile in the 1W case. However, when increasing the number of intervening water layers from 1W to 2W, the particle/slab misalignment becomes a dominant factor controlling the behavior of the PMF profile. We also establish a method to quantify misalignment between the particle and the slab, which exhibits a strong correlation with the free energy for the 2W case. These results shed more light into the roles of particle/slab misalignment and hydrogen bond network in the OA of mineral particles in aqueous solution.
Using first-principles density functional theory (DFT) methods and size-converged supercell models, we analyze the electronic and atomic structure of magnetic $3d$ transition metal dopants in cubic gallium nitride (c-GaN). All stable defect charge states for Fermi levels across the full experimental gap are computed using a method that correctly resolves the boundary condition problem (without a jellium approximation) and eliminates finite-size errors. The resulting computed defect levels are not impacted by the DFT band-gap problem, they span a width consistent with the experimental gap rather than being limited to the single-particle DFT gap. All defects with electronically degenerate (half-metal) $T$d ground states are found to have significant distortions, relaxing to $D$2d structures driven by the Jahn-Teller instability. This leads to insulating ground states for all substitutional $3d$ dopants, refuting claims in the literature that +$U$ or hybrid functional methods are required to avoid artificial half-metal results. Interpreting the $d$n atomic occupations within a crystal-field model and exchange splittings, we identify a systematic trend across the $3d$ transition metal series. Approaches to estimate excited-state energies as observed in photoluminescence from defect centers are assessed, ranging from a Koopmans-type single-particle energy interpretation to relaxed total energy differences in fully self-consistent DFT. The single-particle interpretations are found to be qualitatively predictive and the calculations are consistent with the limited available experimental data across the $3$d dopant series. These results provide a baseline understanding to guide future studies and a conceptual framework within which to interpret new results.
Polymorphism and phase transitions in sodium diuranate, Na2U2O7, are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, β- and γ-Na2U2O7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na2U2O7 phase (R3̄m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are for P21/a α-Na2U2O7 and C2/m β-Na2U2O7, respectively. Gibbs free energy calculations reveal that α-Na2U2O7 (P21/a) and β-Na2U2O7 (C2/m) are almost energetically degenerate at low temperature, with β-Na2U2O7 becoming slightly more stable than α-Na2U2O7 as temperature increases. These findings are consistent with XRD data showing a mixture of α and β phases after cooling of γ-Na2U2O7 to room temperature and the observation of a sluggish α → β phase transition above ca. 600 K. A recently observed α-Na2U2O7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct β → γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial β-phase decomposition, followed by recrystallization into γ-phase as temperature increases.
Balogun, Shuaib A.; Lively, Ryan P.; Losego, Mark D.; Ren, Yi; Steiner, Adam M.
Vapor phase infiltration (VPI) is a post-polymerization modification technique that infuses inorganics into polymers to create organic–inorganic hybrid materials with new properties. Much is yet to be understood about the chemical kinetics underlying the VPI process. The aim of this study is to create a greater understanding of the process kinetics that govern the infiltration of trimethyl aluminum (TMA) and TiCl4 into PMMA to form inorganic-PMMA hybrid materials. To gain insight, this paper initially examines the predicted results for the spatiotemporal concentrations of inorganics computed from a recently posited reaction–diffusion model for VPI. This model provides insight on how the Damköhler number (reaction versus diffusion rates) and non-Fickian diffusional processes (hindering) that result from the material transforming from a polymer to a hybrid can affect the evolution of inorganic concentration depth profiles with time. Subsequently, experimental XPS depth profiles are collected for TMA and TiCl4 infiltrated PMMA films at 90 °C and 135 °C. The functional behavior of these depth profiles at varying infiltration times are qualitatively compared to various computed predictions and conclusions are drawn about the mechanisms of each of these processes. TMA infiltration into PMMA appears to transition from a diffusion-limited process at low temperatures (90 °C) to a reaction-limited process at high temperatures (135 °C) for the film thicknesses investigated here (200 nm). While TMA appears to fully infiltrate these 200 nm PMMA films within a few hours, TiCl4 infiltration into PMMA is considerably slower, with full saturation not occurring even after 2 days of precursor exposure. Infiltration at 90 °C is so slow that no clear conclusions about mechanism can be drawn; however, at 135 °C, the TiCl4 infiltration into PMMA is clearly a reaction-limited process, with TiCl4 permeating the entire thickness (at low concentrations) within only a few minutes, but inorganic loading continuously increasing in a uniform manner over a course of 2 days. Near-surface deviations from the uniform-loading expected for a reaction-limited process also suggest that diffusional hindering is high for TiCl4 infiltration into PMMA. In conclusion, these results demonstrate a new, ex situ analysis approach for investigating the rate-limiting process mechanisms for vapor phase infiltration.
The propagation of self-sustained formation reactions in sputter-deposited Co/Al multilayers is known to exhibit a design-dependent instability. Multilayers having thin bilayers (<55 nm period) exhibit stable propagating waves, whereas those with a larger period react unstably. The specific two-dimensional (2D) instability observed involves the transverse propagation of a band in front of a stalled front commonly referred to as a “spin band.” Previous finite-element studies have shown that these instabilities are thermodynamically driven by the forward conduction of heat away from the flame front. However, the magnitude of that loss is inherently tied to the bilayer design in traditional bimetallic multilayers, which couples any proposed stability criteria to a varying critical diffusion distance. This work utilizes a recently developed class of materials known as “inert-mediated reactive multilayers” to decouple the thermodynamic and kinetic contributions to propagating wave stability by reducing the stored chemical energy density in normally stable bilayer designs. By depositing an inert product phase (B2-CoAl) within the mid-plane of Co and Al reactant layers, spin instabilities arise as a function of both diluted volume and critical diffusion distance. From there, a stability criterion is determined for Co/Al multilayers based on enthalpy loss from the reaction zone, and its physical significance is explored.
Wexler, Robert B.; Sai Gautam, Gopalakrishnan; Bell, Robert T.; Shulda, Sarah; Strange, Nicholas A.; Trindell, Jamie T.; Sugar, Joshua D.; Nygren, Eli; Sainio, Sami; Mcdaniel, Anthony H.; Ginley, David; Carter, Emily A.; Stechel, Ellen B.
Modeling-driven design of redox-active off-stoichiometric oxides for solar thermochemical H2 production (STCH) seldom has resulted in empirical demonstration of competitive materials. We report the theoretical prediction and experimental evidence that the perovskite Ca2/3Ce1/3Ti1/3Mn2/3O3 is synthesizable with high phase purity, stable, and has desirable redox thermodynamics for STCH, with a predicted average neutral oxygen vacancy (VO) formation energy, Ev = 3.30 eV. Flow reactor experiments suggest potentially comparable or greater H2 production capacity than recent promising Sr-La-Mn-Al and Ba-Ce-Mn metal oxide perovskites. Utilizing quantum-based modeling of a solid solution on both A and B sub-lattices, we predict the impact of nearest-neighbor composition on Ev and determine that A-site Ce4+ reduction dominates the redox-activity of Ca2/3Ce1/3Ti1/3Mn2/3O3. X-ray absorption spectroscopy measurements provide evidence that supports these predictions and reversible Ce4+-to-Ce3+ reduction. Our models predict that Ce4+ reduces even when it is not nearest-neighbor to the VO, suggesting that refinement of Ce stoichiometry has the possibility of further enhancing performance.
Measured salt compositions in dust collected over roughly the last decade from surfaces of in-service stainless-steel alloys at four locations around the United States are presented, along with the predicted brine compositions that would result from deliquescence of these salts. The salt compositions vary greatly from ASTM seawater and from laboratory salts (i.e., NaCl or MgCl2) commonly used on corrosion testing. The salts contained relatively high amounts of sulfates and nitrates, evolved to basic pH values, and exhibited deliquescence relative humidity values (RH) higher than seawater. Additionally, inert dust in components were quantified and considerations for laboratory testing are presented. The observed dust compositions are discussed in terms of the potential corrosion behavior and are compared to commonly used accelerated testing protocols. Finally, ambient weather conditions and their influence on diurnal fluctuations in temperature (T) and RH on heated metal surfaces are evaluated and a relevant diurnal cycle for laboratory testing a heated surface has been developed. Suggestions for future accelerated tests are proposed that include exploration of the effects of inert dust particles on atmospheric corrosion, chemistry considerations, and realistic diurnal fluctuations in T and RH. Understanding mechanisms in both realistic and accelerated environments will allow development of a corrosion factor (i.e., scaling factor) for the extrapolation of laboratory-scale test results to real world applications.
Machine learning is on a bit of a tear right now, with advances that are infiltrating nearly every aspect of our lives. In the domain of materials science, this wave seems to be growing into a tsunami. Yet, there are still real hurdles that we face to maximize its benefit. This Matter of Opinion, crafted as a result of a workshop hosted by researchers at Sandia National Laboratories and attended by a cadre of luminaries, briefly summarizes our perspective on these barriers. By recognizing these problems in a community forum, we can share the burden of their resolution together with a common purpose and coordinated effort.
The burgeoning field of materials informatics necessitates a focus on educating the next generation of materials scientists in the concepts of data science, artificial intelligence (AI), and machine learning (ML). In addition to incorporating these topics in undergraduate and graduate curricula, regular hands-on workshops present the most effective medium to initiate researchers to informatics and have them start applying the best AI/ML tools to their own research. With the help of the Materials Research Society (MRS), members of the MRS AI Staging Committee, and a dedicated team of instructors, we successfully conducted workshops covering the essential concepts of AI/ML as applied to materials data, at both the Spring and Fall Meetings in 2022, with plans to make this a regular feature in future meetings. Here, in this article, we discuss the importance of materials informatics education via the lens of these workshops, including details such as learning and implementing specific algorithms, the crucial nuts and bolts of ML, and using competitions to increase interest and participation.
In magneto-inertial fusion, the ratio of the characteristic fuel length perpendicular to the applied magnetic field R to the α-particle Larmor radius Q α is a critical parameter setting the scale of electron thermal-conduction loss and charged burn-product confinement. Using a previously developed deep-learning-based Bayesian inference tool, we obtain the magnetic-field fuel-radius product B R ∝ R / Q α from an ensemble of 16 magnetized liner inertial fusion (MagLIF) experiments. Observations of the trends in BR are consistent with relative trade-offs between compression and flux loss as well as the impact of mix from 1D resistive radiation magneto-hydrodynamics simulations in all but two experiments, for which 3D effects are hypothesized to play a significant role. Finally, we explain the relationship between BR and the generalized Lawson parameter χ. Our results indicate the ability to improve performance in MagLIF through careful tuning of experimental inputs, while also highlighting key risks from mix and 3D effects that must be mitigated in scaling MagLIF to higher currents with a next-generation driver.
The accurate quantification of the performance loss rate of photovoltaic systems is critical for project economics. Following the current research activities in the photovoltaic performance and reliability field, this work presents a comparative assessment between common change point methods for performance loss rate estimation of fielded photovoltaic installations. An extensive testing campaign was thus performed to evaluate time series analysis approaches for performance loss rate evaluation of photovoltaic systems. Historical electrical data from eleven photovoltaic systems installed in Nicosia, Cyprus, and the locations’ meteorological measurements over a period of 8 years were used for this investigation. The application of change point detection algorithms on the constructed monthly photovoltaic performance ratio series revealed that the obtained trend might not always be linear. Specifically, thin film photovoltaic systems showed nonlinear behavior, while nonlinearities were also detected for some crystalline silicon photovoltaic systems. When applying several change point techniques, different numbers and locations of changes were detected, resulting in different performance loss rate values (varying by up to 0.85%/year even for the same number of change points). The results highlighted the importance of the application of nonlinear techniques and the need to extract a robust nonlinear model for detecting significant changes in time series data and estimating accurately the performance loss rate of photovoltaic installations.
The Sandia National Laboratories, in California (Sandia/CA) is a research and development facility, owned by the U.S. Department of Energy’s National Nuclear Security Administration agency (DOE/NNSA). The laboratory is located in the City of Livermore (the City) and is comprised of approximately 410 acres. The Sandia/CA facility is operated by National Technology and Engineering Solutions of Sandia, LLC (NTESS) under a contract with the DOE/NNSA. The DOE/ NNSA’s Sandia Field Office (SFO) oversees the operations of the site. North of the Sandia/CA facility is the Lawrence Livermore National Laboratory (LLNL), in which Sandia/CA’s sewer system combines with before discharging to the City’s Publicly Owned Treatment Works (POTW) for final treatment and processing. The City’s POTW authorizes the wastewater discharge from Sandia/CA via the assigned Wastewater Discharge Permit #1251 (the Permit), which is issued to the DOE/NNSA’s main office for Sandia National Laboratories, located in New Mexico (Sandia/NM). The Permit requires the submittal of this Monthly Sewer Monitoring Report to the City by the twenty-fifth day of each month.
High speed analog-to-digital converters (ADC), switched-capacitor delay elements, and pulsed radio frequency (RF) systems all require switches in the signal path operating at high switching speeds, providing low resistance when enabled, and providing high signal isolation when disabled. In semiconductor technologies such as CMOS, the enabled state resistance directly scales with the sizing of the switch device, where a larger width switch provides a lower enabled state resistance. As the device width is increased, so is the capacitance formed between the gate, drain, and source of the device.
The Z Machine at Sandia National Laboratories is a pulsed power facility for high-energy density physics experiments that can shock materials to extreme temperatures and pressures through a focused energy release of up to ∼ 25 MJ in < 100 nanoseconds. It has been in operation for more than two decades and conducts up to ∼ 100 experiments, or “shots,” per year. Based on a set of 74 known shot times from 2018, we determined that Z Machine shots produce detectable ∼ 3–17 Hz ground motion 12 km away at the Albuquerque Seismological Laboratory, New Mexico (ANMO), borehole seismograph, with peak signal at ∼ 7 Hz. The known shot waveforms were used to create a three-component template, leading to the detection of 2339 Z Machine shots since 1998 through single-station cross-correlation. Local seismic magnitude estimates range from local magnitude (ML) -2 to -1.3 and indicate that only a small fraction of the shot energy is transmitted by seismic phases observable at 12 km distance. The most recent major facility renovation, which was intended to decrease mechanical dissipation, is associated with an abrupt decrease in observed seismic amplitudes at ANMO despite stable maximum shot energy. The highly repetitive impulsive sources are well suited to coda-wave interferometry to investigate time-dependent velocity structures. Relative velocity variations (dv/v) show an annual cycle with amplitude of ∼ 0.2%. Local minima are observed in the late spring, and dv/v increases through the summer monsoon rainfall, possibly reflecting patchy saturation as rainfall infiltrates near the eastern edge of the Albuquerque basin. The cumulative results demonstrate that forensic seismology can provide insight into long-term operation of facilities such as pulsed-power laboratories, and that their recurring signals may be valuable for studies of time-dependent structure.
X-ray diffraction (XRD) is a necessary technique for understanding states of materials under static and dynamic loading conditions. The higher-pressure Equation of State (EOS) of many materials can only be explored via shock or ramp compression at temperatures and pressures of interest. While static XRD work has yielded EOS measurements in the 100 - 200 GPa regime, dynamic X-ray diffraction (DXRD) can explore EOS phases in the TPa regime, which closely resembles inner-core planetary conditions. DXRD hinges on the ability to measure the exact phase or phase change of a material while under dynamic loading conditions. Macroscopic diagnostic systems (e.g. velocimetry and pyrometry) can infer a phase change but not identify the specific phase entered by a material. While microscopic (atomic-level) diagnostic systems (e.g. DXRD) have been designed and implemented in Department of Energy’s (DOE) National Laboratories complex, the unique nature of Sandia National Laboratories’ Pulsed Power Facility (Z Machine) prohibits the use of such devices. The destructive nature of Z experiments presents a challenge to data capture and retrieval. Furthermore there are electromagnetic interference, X-ray background, and mechanical constraints to consider. Thus, a multi-part X-ray diagnostic for use on the Z Machine and Z-Beamlet Laser system has been designed and analyzed. Portions of this new DYnamic SCintillator Optic (DYSCO) have been built, tested and fielded. A data analysis software has been written. Finally, the radiance profile of the DYSCO’s scintillator has been characterized through experiments performed at the University of Arizona.
Induced seismicity is an inherent risk associated with geologic carbon storage (GCS) in deep rock formations that could contain undetected faults prone to failure. Modeling-based risk assessment has been implemented to quantify the potential of injection-induced seismicity, but typically simplified multiscale geologic features or neglected multiphysics coupled mechanisms because of the uncertainty in field data and computational cost of field-scale simulations, which may limit the reliable prediction of seismic hazard caused by industrial-scale CO2 storage. The degree of lateral continuity of the stratigraphic interbedding below the reservoir and depth-dependent fault permeability can enhance or inhibit pore-pressure diffusion and corresponding poroelastic stressing along a basement fault. This study presents a rigorous modeling scheme with optimal geological and operational parameters needed to be considered in seismic monitoring and mitigation strategies for safe GCS.
Interval Assignment (IA) is the problem of selecting the number of mesh edges (intervals) for each curve for conforming quad and hex meshing. The intervals x is fundamentally integer-valued. Many other approaches perform numerical optimization then convert a floating-point solution into an integer solution, which is slow and error prone. We avoid such steps: we start integer, and stay integer. Incremental Interval Assignment (IIA) uses integer linear algebra (Hermite normal form) to find an initial solution to the meshing constraints, satisfying the integer matrix equation Ax=b. Solving for reduced row echelon form provides integer vectors spanning the nullspace of A. We add vectors from the nullspace to improve the initial solution, maintaining Ax=b. Heuristics find good integer linear combinations of nullspace vectors that provide strict improvement towards variable bounds or goals. IIA always produces an integer solution if one exists. In practice we usually achieve solutions close to the user goals, but there is no guarantee that the solution is optimal, nor even satisfies variable bounds, e.g. has positive intervals. We describe several algorithmic changes since first publication that tend to improve the final solution. The software is freely available.
Deploying Tidal Energy Converters for electricity generation requires prior-knowledge of the potential Annual Energy Production (AEP) at the site, Ideally using a year-long tidal current record at the proposed site to minimize uncertainty. However, such records are often unavailable. Fortunately, using the periodic nature of tidal variability, the International Electrotechnical Commission Technical Specification for tidal energy resource assessment requires AEP calculation using at least 90 days of tidal current records at each turbine location. The sensitivity of AEP to different record durations has not been fully assessed. This is the goal of our study. The study utilized the U.S. tidal energy geodatabase to simulate tidal currents with various lengths, during 100 years of the 21st century. We then consider two frameworks for evaluating AEP: (a) The long-term (months) fixed instrument (FI) measurement at each proposed tidal turbine location, and (b) one FI measurement and short-term (hours) boat-based moving vessel measurements. Under the two scenarios, we examine the AEP assessed from short tidal current records, including how the AEP uncertainties vary spatially and temporally, and how they are associated with various astronomical factors. This helps provide guidance on choosing the appropriate assessment methodologies to reduce the AEP uncertainties and project cost.
Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.
Villa, Daniel L.; Schostek, Tyler; Govertsen, Krissy; Macmillan, Madeline
Applying extreme temperature events for future conditions is not straightforward for infrastructure resilience analyses. This work introduces a stochastic model that fills this gap. The model uses at least 50 years of daily extreme temperature records, climate normals with 10%–90% confidence intervals, and shifts/offsets for increased frequency and intensity of heat wave events. Intensity and frequency are shifted based on surface temperature anomaly from 1850–1900 for 32 models from CMIP6. A case study for Worcester, Massachusetts passed 85% of cases using the two-sided Kolmogorov–Smirnov p-value test with 95% confidence for both temperature and duration. Future shifts for several climate scenarios to 2020, 2040, 2060, and 2080 had acceptable errors between the shifted model and 10- and 50-year extreme temperature event thresholds with the largest error being 2.67°C. The model is likely to be flexible enough for other patterns of extreme weather such as extreme precipitation and hurricanes.
Solving large number of small linear systems is increasingly becoming a bottleneck in computational science applications. While dense linear solvers for such systems have been studied before, batched sparse linear solvers are just starting to emerge. In this paper, we discuss algorithms for solving batched sparse linear systems and their implementation in the Kokkos Kernels library. The new algorithms are performance portable and map well to the hierarchical parallelism available in modern accelerator architectures. The sparse matrix vector product (SPMV) kernel is the main performance bottleneck of the Krylov solvers we implement in this work. The implementation of the batched SPMV and its performance are therefore discussed thoroughly in this paper. The implemented kernels are tested on different Central Processing Unit (CPU) and Graphic Processing Unit (GPU) architectures. We also develop batched Conjugate Gradient (CG) and batched Generalized Minimum Residual (GMRES) solvers using the batched SPMV. Our proposed solver was able to solve 20,000 sparse linear systems on V100 GPUs with a mean speedup of 76x and 924x compared to using a parallel sparse solver with a block diagonal system with all the small linear systems, and compared to solving the small systems one at a time, respectively. We see mean speedup of 0.51 compared to dense batched solver of cuSOLVER on V100, while using lot less memory. Thorough performance evaluation on three different architectures and analysis of the performance are presented.
The purpose of this protocol is to define procedures and practices to be used by the PACT center for field testing of metal halide perovskite (MHP) photovoltaic (PV) modules. The protocol defines the physical, electrical, and analytical configuration of the tests and applies equally to mounting systems at a fixed orientation or sun tracking systems. While standards exist for outdoor testing of conventional PV modules, these do not anticipate the unique electrical behavior of perovskite cells. Further, the existing standards are oriented toward mature, relatively stable products with lifetimes that can be measured on the scale of years to decades. The state of the art for MHP modules is still immature with considerable sample to sample variation among nominally identical modules. Version 0.0 of this protocol does not define a minimum test duration, although the intent is for modules to be fielded for periods ranging for weeks to months. This protocol draws from relevant parts of existing standards, and where necessary includes modifications specific to the behavior of perovskites.
A new non-neutral generalized Ohm's law (GOL) model for atomic plasmas is presented. This model differs from previous models of this type in that quasi-neutrality is not assumed at any point. Collisional effects due to ionization, recombination, and elastic scattering are included, and an expression for the associated plasma conductivity is derived. An initial set of numerical simulations are considered that compare the GOL model to a two-fluid model in the ideal (collisionless) case. The results demonstrate that solutions obtained from the two models are essentially indistinguishable in most cases when the ion-electron mass ratio is within the range of physical values for atomic plasmas. Additionally, some limitations of the model are discussed.
Protecting against multi-step attacks of uncertain start times and duration forces the defenders into indefinite, always ongoing, resource-intensive response. To allocate resources effectively, the defender must analyze and respond to an uncertain stream of potentially undetected multiple multi-step attacks and take measures of attack and response intensity over time into account. Such response requires estimation of overall attack success metrics and evaluating effect of defender strategies and actions associated with specific attack steps on overall attack metrics. We present a novel game-theoretic approach GPLADD to attack metrics estimation and demonstrate it on attack data derived from MITRE's ATT&CK Framework and other sources. In GPLADD, the time to complete attack steps is explicit; the attack dynamics emerges from attack graph and attacker-defender capabilities and strategies and therefore reflects 'physics' of attacks. The time the attacker takes to complete an attack step is drawn from a probability distribution determined by attacker and defender strategies and capabilities. This makes time a physical constraint on attack success parameters and enables comparing different defender resource allocation strategies across different attacks. We solve for attack success metrics by approximating attacker-defender games as discrete-time Markov chains and show evaluation of return on detection investments associated with different attack steps. We apply GPLADD to MITRE's APT3 data from ATT&CK Framework and show that there are substantial and un-intuitive differences in estimated real-world vendor performance against a simplified APT3 attack. We focus on metrics that reflect attack difficulty versus attacker ability to remain hidden in the system after gaining control. This enables practical defender optimization and resource allocation against multi-step attacks.
This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of 192 µm. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier-Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant C2 = 2 compared to the accepted value C2 = 2/3. Results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.
As part of the project “Designing Resilient Communities (DRC): A Consequence-Based Approach for Grid Investment,” funded by the United States (US) Department of Energy’s (DOE) Grid Modernization Laboratory Consortium (GMLC), Sandia National Laboratories (Sandia) partnered with a variety of government, industry, and university participants to develop and test a framework for community resilience planning focused on modernization of the electric grid. This report provides a summary of the development, description, and demonstration of the resulting Resilient Community Design Framework.
Azizur-Rahman, Khalifa M.; Mah, Jasmine J.; Liang, Baolai; Huffaker, Diana L.; Nolde, Jill; Aifer, Edward; Hun Park, Jeung; Kim, Richard S.; Dahl, Russel
Diesel piston-bowl shape is a key design parameter that affects spray-wall interactions and turbulent flow development, and in turn affects the engine’s thermal efficiency and emissions. It is hypothesized that thermal efficiency can be improved by enhancing squish-region vortices as they are hypothesized to promote fuel-air mixing, leading to faster heat-release rates. However, the strength and longevity of these vortices decrease with advanced injection timings for typical stepped-lip (SL) piston geometries. Dimple stepped-lip (DSL) pistons enhance vortex formation at early injection timings. Previous engine experiments with such a bowl show 1.4% thermal efficiency gains over an SL piston. However, soot was increased dramatically [SAE 2022-01-0400]. In a previous study, a new DSL bowl was designed using non-combusting computational fluid dynamic simulations. This improved DSL bowl is predicted to promote stronger, more rotationally energetic vortices than the baseline DSL piston: it employs shallower, narrower, and steeper-curved dimples that are placed further out into the squish region. In the current experimental study, this improved bowl is tested in a medium-duty diesel engine and compared against the SL piston over an injection timing sweep at low-load and part-load operating conditions. No substantial thermal efficiency gains are achieved at the early injection timing with the improved DSL design, but soot emissions are lowered by 45% relative to the production SL piston, likely due to improved air utilization and soot oxidation. However, these benefits are lost at late injection timings, where the DSL piston renders a lower thermal efficiency than that of the SL piston. Energy balance analyses show higher wall heat transfer with the DSL piston than with the SL piston despite a 1.3% reduction in the piston surface area. Vortex enhancement may not necessarily lead to improved efficiency as more energetic squish-region vortices can lead to higher convective heat transfer losses.
In this work, we use the Brillouin flow analytic framework to examine the physics of Magnetically Insulated Transmission Lines (MITL). We derive a model applicable to any particle species, including both positive and negative ions, in planar and cylindrical configurations. We then show how to self-consistently solve for two-species simultaneously, using magnetically insulated electrons and positive ions as an example. We require both layers to be spatially separated and magnetically insulated (mutually magnetically insulated); for a 7.5 cm gap with a 2 MV bias voltage, this condition requires magnetic fields in excess of 2.73 T. We see a close match between mutually insulated MITL performance and “superinsulated” (high degree of magnetic insulation) electron-only theory, as may be expected for these high magnetic fields. However, the presence of ions leads to several novel effects: (1) Opposite to electron-only theory, total electron currents increase rather than decrease as the degree of magnetic insulation becomes stronger. The common assumption of neglecting electrons for superinsulated MITL operation must be revisited when ions are present—we calculate up to 20× current enhancement. (2) The electron flow layer thickness increases up to double, due to ion space-charge enhancement. (3) The contributions from both ions and electrons to the MITL flow impedance are calculated. The flow impedance drops by over 50% when ions fill the gap, which can cause significant reflections at the load if not anticipated and degrade performance. Additional effects and results from the inclusion of the ion layer are discussed.
This document provides the instructions for participating in the 2021 blind photovoltaic (PV) modeling intercomparison organized by the PV Performance Modeling Collaborative (PVPMC). It describes the system configurations, metadata, and other information necessary for the modeling exercise. The practical details of the validation datasets are also described. The datasets were published online in open access in April 2023, after completing the analysis of the results.
We propose primal–dual mesh optimization algorithms that overcome shortcomings of the standard algorithm while retaining some of its desirable features. “Hodge-Optimized Triangulations” defines the “HOT energy” as a bound on the discretization error of the diagonalized Delaunay Hodge star operator. HOT energy is a natural choice for an objective function, but unstable for both mathematical and algorithmic reasons: it has minima for collapsed edges, and its extrapolation to non-regular triangulations is inaccurate and has unbounded minima. We propose a different extrapolation with a stronger theoretical foundation, and avoid extrapolation by recalculating the objective just beyond the flip threshold. We propose new objectives, based on normalizations of the HOT energy, with barriers to edge collapses and other undesirable configurations. We propose mesh improvement algorithms coupling these. When HOT optimization nearly collapses an edge, we actually collapse the edge. Otherwise, we use the barrier objective to update positions and weights and remove vertices. By combining discrete connectivity changes with continuous optimization, we more fully explore the space of possible meshes and obtain higher quality solutions.
CuCr2O4 spinel is a candidate coating material for central receivers in concentrating solar power to protect structural alloys against high temperature oxidation and related degradation. Coating performance and microstructure of dip-coated and sintered coatings is dictated by the initial particle size of the CuCr2O4 and sintering temperature, but can be compromised by particle agglomeration. Here in this study, sub-micron particles were synthesised through the Pechini and modified Pechini sol–gel methods. Phase composition was confirmed via X-ray diffraction. Particle growth during calcination of the nanoparticles at different temperatures (650°C, 750°C, 850°C) and times (between 1 and 24 h) was measured via laser diffraction and scanning electron microscopy. The modified Pechini method displayed evidence of smaller particle sizes and greater agglomeration. The kinetics of particle growth observed are consistent with a diffusion limited inhibited grain growth model.
Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.
HyRAM+ is a toolkit that includes fast-running models for the unconstrained (i.e., no wall interactions) dispersion and flames for non-premixed fuels. The models were developed for use with hydrogen, but the toolkit was expanded to include propane and methane in a recent release. In this work we validate the dispersion and flame models for these additional fuels, based on reported literature data. The validation efforts spanned a range of release conditions, from subsonic to underexpanded jets and flames for a range of mass flow rates. In general, the dispersion model works well for both propane and methane although the width of the jet/plume is predicted to be wider than observed in some cases. The flame model tends to over-predict the induced buoyancy for low-momentum flames, while the radiative heat flux agrees with the experimental data reasonably well, for both fuels. The models could be improved but give acceptable predictions for propane and methane behavior for the purposes of risk assessment.