The National Nuclear Security Agency (NNSA) created a Minority Serving Institution Partnership Plan (MSIPP) to 1) align investments in a university capacity and workforce development with the NNSA mission to develop the needed skills and talent for NNSA's enduring technical workforce at the laboratories and production plants, and 2) to enhance research and education at underrepresented colleges and universities. Out of this effort, MSIPP launched a new program in early FY17 focused on Tribal Colleges and Universities (TCUs). The following report summarizes the project focus and status updates during this reporting period.
The catastrophic nuclear reactor accident at Fukushima damaged public confidence in nuclear energy and a demand for new engineered safety features that could mitigate or prevent radiation releases to the environment in the future. We have developed a novel use of sacrificial material (SM) to prevent the molten corium from breaching containment during accidents as well as a validated, novel, high-fidelity modeling capability to design and optimize the proposed concept. Some new reactor designs employ a core catcher and a SM, such as ceramic or concrete, to slow the molten corium and avoid the breach of the containment. However, existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials (current designs are limited to water). The SM proposed in this Laboratory Development Research and Development (LDRD) project is based on granular carbonate minerals that can be used in existing light water reactor plants. This new SM will induce an endothermic reaction to quickly freeze the corium in place, with minimal hydrogen explosion and maximum radionuclide retention. Because corium spreading is a complex process strongly influenced by coupled chemical reactions (with underlying containment material and especially with the proposed SM), decay heat and phase change. No existing tool is available for modeling such a complex process. This LDRD project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. We have demonstrated small-scale to large-scaled experiments using lead oxide (Pb0) as surrogate for molten corium, which showed that the reaction of the SM with molten Pb0 results in a fast solidification of the melt and the formation of open pore structures in the solidified Pb0 because of CO2 released from the carbonate decomposition.
This report describes how random pileup calculations are performed by the Gamma Detector Response and Analysis Software (GADRAS) Version 19.1. The computational approach and examples are presented for gamma-ray detectors with and without pileup rejectors. This pileup algorithm executes more quickly and the results are more accurate than previous versions of GADRAS. The detector response function can be refined to characterize distortions in peak shapes that occur at high-count rates. The empirical refinement can also be applied to describe the response of partially-effective pileup rejectors. Implications are discussed for the analysis of both static measurements and dynamic collections of the type acquired with radiation portals. ACKNOWLEDGEMENTS This work was funded by the Defense Threat Reduction Agency (DTRA) and the Department of Homeland Security (DHS) Counter Weapons of Mass Destruction (CWMD) office.
The Terry Turbine Expanded Operating Band Project is currently conducting testing at Texas A&M University, and the resulting data has been incorporated into MELCOR models of the Terry turbines used in nuclear power plants. These improved models have produced improvements in the Fukushima Daiichi Unit 2 simulations while providing new insights into the behavior of the plant. The development of future experimental test efforts is ongoing. Development of and refinements to the plans for full-scale steam and steam-water turbine ingestion testing has been performed. These full-scale steam-based tests will complement the testing occurring at Texas A&M University, and will resolve the remaining questions regarding scale or working fluid. Planning work has also begun for future testing intended to explore the uncontrolled RCIC self-regulation theorized to have occurred in Fukushima Daiichi Unit 2.
This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 11932 and 14959. Revision 11932 represents the last official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 14959. Along with the newly updated MELCOR Users' Guide and Reference Manual, users will be aware and able to assess the new capabilities for their modeling and analysis applications. Following the official release an addendum section has been added to this report detailing modifications made to the official release which support the accompanying patch release. The addendums address user reported issues and previously known issues within the official code release which extends the original Quick look document to also support the patch release. Furthermore, the addendums section documents the recent changes to input records in the Users' Guide applicable to the patch release and corrects a few issues in the revision 14959 release as well.
Crystal plasticity-finite element method (CP-FEM) is now widely used to understand the mechanical response of polycrystalline materials. However, quantitative mesh convergence tests and verification of the necessary size of polycrystalline representative volume elements (RVE) are often overlooked in CP-FEM simulations. Mesh convergence studies in CP-FEM models are more challenging compared to conventional finite element analysis (FEA) as they are not only computationally expensive but also require explicit discretization of individual grains using many finite elements. Resolving each grains within a polycrystalline domain complicates mesh convergence study since mesh convergence is strongly affected by the initial crystal orientations of grains and local loading conditions. In this work, large-scale CP-FEM simulations of single crystals and polycrystals are conducted to study mesh sensitivity in CP-FEM models. Various factors that may affect the mesh convergence in CP-FEM simulations, such as initial textures, hardening models and boundary conditions are investigated. In addition, the total number of grains required to obtain adequate RVE is investigated. This work provides a list of guidelines for mesh convergence and RVE generation in CP-FEM modeling.
Sandia National Laboratories (SNL) National Solar Thermal Test Facility (NSTTF) and Tech Library have been collaborating over the course of the FY19 period to establish and maintain the first and only digital collection in the world of Concentrating Solar Power (CSP) related historical documents, dating back to the CSP program inception here at Sandia in the 1970's thru to the present. The unclassified, unrestricted (UUR) collection, comprised of internally generated Sandia documents as well as a significant number of external reports will be searchable via both the Sandia website and OSTI, DOE's document repository. DOE is currently championing efforts to get the collection launched, where international partners, which include Australia and Germany, plan to forward related documents to be included in the CSP archive. Advancing this transformative project will make the CSP collection accessible to the Sandia and global communities.
Chen, Yizhang; Cogulu, Egecan; Roy, Debangsu; Ding, Jinjun; Mohammadi, Jamileh B.; Kotula, Paul G.; Missert, Nancy A.; Wu, Mingzhong; Kent, Andrew D.
We present a study of the transport properties of thermally generated spin currents in an insulating ferrimagnetic-antiferromagnetic-ferrimagnetic trilayer over a wide range of temperature. Spin currents generated by the spin Seebeck effect (SSE) in a yttrium iron garnet (YIG) YIG/NiO/YIG trilayer on a gadolinium gallium garnet (GGG) substrate were detected using the inverse spin Hall effect (ISHE) in Pt. By studying samples with different NiO thicknesses, the spin diffusion length of NiO was determined to be ∼3.8 nm at room temperature. Surprisingly, a large increase of the SSE signal was observed below 30 K, and the field dependence of the signal closely follows a Brillouin function for an S=7/2 spin. The increase of the SSE signal at low temperatures could thus be associated with the paramagnetic SSE from the GGG substrate. Besides, a broad peak in the SSE response was observed around 100 K. These observations are important in understanding the generation and transport properties of spin currents through magnetic insulators and the role of a paramagnetic substrate in spin current generation.
This research objective of this EELDRD study was to learn to electrodeposit Pt Au alloys with independently controlled composition and grain size. What was accomplished was the capability to electrodeposit PtAu alloys with controlled composition and a nanocrysolline grain size. Nanocrystalline metals as a class and, specifically, the Pt0.9Au0.1 alloy developed in 2015-17 via sputtering at Sandia National Labs have clear advantages in strength, wear resistance, and fatigue tolerance over commercially-available structural alloys. With this capability befitting coating of complex components and implementable at existing vendors, we can upgrade the electrical contact component reliability of selected Labs systems.
Cylindrical dog-bone (or dumbbell) shaped samples have become a common design for dynamic tensile tests of ductile materials with a Kolsky tension bar. When a direct measurement of displacement between the bar ends is used to calculate the specimen strain, the actual strain in the specimen gage section is overestimated due to strain in the specimen shoulder and needs to be corrected. The currently available correction method works well for elastic-perfectly plastic materials but may not be applicable to materials that exhibit significant work-hardening behavior. In this study, we developed a new specimen strain correction method for materials possessing an elastic-plastic with linear work-hardening stress–strain response. A Kolsky tension bar test of a Fe-49Co-2V alloy (known by trade names Hiperco and Permendur) was used to demonstrate the new specimen strain correction method. This new correction method was also used to correct specimen strains in Kolsky tension bar experiments on two other materials: 4140 alloy, and 304L-VAR stainless steel, which had different work-hardening behavior.
Slabs of cast calcium sulfate dihydrate have been used for decades as a heat transfer barrier in commercial and residential construction and are assessed in accordance with the parameters defined by ASTM C1396. The study described herein hypothesizes that minor impurities have significant effects on the high-temperature performance of these casts. Five thermocouples are cast into a 20-mm-thick slab at approximately 4-mm intervals are measured between ambient and approximately 1000°C in a furnace. Thermal diffusivity and inertia are estimated from this temperature profile. This study compares natural, flue gas desulfurized (FGD) and reagent grade hemihydrates that have been converted into the dihydrate form with distilled water; the common additives kaolin and borax are used to modify the cast. The thermocouple data allow an effective thermal diffusivity (α′) and an effective thermal inertia (I′) to be calculated. The hemihydrate source and kaolin content are found to affect the high-temperature performance; the FGD source increases the thermal inertia, and kaolin inhibits the formation of other borate compounds through intercalation.
Hansen, Nils H.; Leon, Larisa; Ruwe, Lena; Moshammer, Kai; Seidel, Lars; Shrestha, Krishna P.; Wang, Xiaoxiao; Mauss, Fabian
A comprehensive, chemically detailed mechanism for the combustion of 2-methyl-2-butene and n-pentane is presented to provide insights into the different sooting tendencies of these two structurally different C5 hydrocarbons. A hierarchically assembled mechanism has been developed to specifically target speciation data from low-pressure premixed flames of 2-methyl-2-butene [Ruwe et al., Combust. Flame, 175, 34-46, 2017] and newly measured mole fraction data for a fuel-rich (ɸ=1.8) n-pentane flame, in which species profiles up to phenol were quantified. The partially isomer-resolved chemical composition of this flame was determined using flame-sampling molecular-beam mass spectrometry with single-photon ionization by tunable, synchrotron-generated vacuum-ultraviolet radiation. The presented model, which includes a newly determined, consistent set of the thermochemistry data for the C5 species, presents overall satisfactory capabilities to predict the mole fraction profiles of common combustion intermediates. The analysis of the model predictions revealed the fuel-structure dependencies (i.e. saturated vs. unsaturated and linear vs. branched) of the formation of small aromatic species that are considered as soot precursors. The propensity of the 2-methyl-2-butene flame to form larger concentrations of aromatic species was traced back to the readily available formation routes of several small precursor molecules and the efficient formation of “first aromatic rings” beyond benzene.
Recent event-related brain potential (ERP) experiments have demonstrated parafoveal N400 expectancy and congruity effects, showing that semantic information can be accessed from words in parafoveal vision (a conclusion also supported by some eye-tracking work). At the same time, it is unclear how higher-order integrative aspects of language comprehension unfold across the visual field during reading. In the current study, we recorded ERPs in a parafoveal flanker paradigm, while readers were instructed to read passively for comprehension or to judge the plausibility of sentences in which target words varied in their semantic expectancy and congruity. We directly replicated prior work showing graded N400 effects for parafoveal viewing, which are then not duplicated when the target words are processed foveally. Critically, although N400 effects were not modulated by task goals, a posteriorly distributed late positive component thought to reflect semantic integration processes was observed to semantic incongruities only in the plausibility judgment task. However, this effect was observed at a considerable delay, appearing only after words had moved into foveal vision. Our findings thus suggest that semantic access can be initiated in parafoveal vision, whereas central foveal vision may be necessary to enact higher-order (and task-dependent) integrative processing.
The Hydrogen Materials—Advanced Research Consortium (HyMARC) is the core storage material research team of the DOE/EERE Fuel Cell Technologies Office (FCTO) and is comprised of Sandia National Laboratories (Livermore, CA; SNL), Lawrence Livermore National Laboratory (LLNL), and Lawrence Berkeley Laboratory (LBNL). Its objective is to overcome critical scientific barriers limiting the use of solid-state materials for vehicular hydrogen storage, thereby enabling design and discovery of breakthrough storage materials. Over the three-year lifetime of the project, HyMARC "moved the bar" relative to compressed gas storage by identifying the most promising material improvement strategies, obtaining thermodynamic data that was either missing or inaccurate in the literature, and filling major gaps in the toolkit of computational models. The HyMARC team also developed many new capabilities in the areas of material synthesis and characterization that address specific roadblocks to discovery of successful storage materials.
High-harmonic generation (HHG) is a signature optical phenomenon of strongly driven, nonlinear optical systems. Specifically, the understanding of the HHG process in rare gases has played a key role in the development of attosecond science1. Recently, HHG has also been reported in solids, providing novel opportunities such as controlling strong-field and attosecond processes in dense optical media down to the nanoscale2. Here, we report HHG from a low-loss, indium-doped cadmium oxide thin film by leveraging the epsilon-near-zero (ENZ) effect3–8, whereby the real part of the material’s permittivity in certain spectral ranges vanishes, as well as the associated large resonant enhancement of the driving laser field. We find that ENZ-assisted harmonics exhibit a pronounced spectral redshift as well as linewidth broadening, resulting from the photo induced electron heating and the consequent time-dependent ENZ wavelength of the material. Our results provide a new platform to study strong-field and ultrafast electron dynamics in ENZ materials, reveal new degrees of freedom for spectral and temporal control of HHG, and open up the possibilities of compact solid-state attosecond light sources.
The creation of microelectromechanical systems (MEMS) that can operate through elevated temperatures would enable systems diagnostics and controls that are not possible with conventional-off-the-shelf components. The integration of silicon carbide (SiC) with aluminum nitride (AlN) has led to the fabrication of devices that can withstand elevated temperature anneals >935 °C. The results from a piezoelectric micromachined ultrasonic transducer (PMUT) and a microresonator are reported as demonstrations of the fabrication process. Testing the PMUT response before and after annealing at 935 °C led to a change in resonant frequency of less than 1%, which is attributable to a shift in film stress. The response of the microresonator was RF tested in situ up to 500 °C and showed no degradation in its electromechanical coupling coefficient. The resonant frequency decreased with temperature due to the temperature coefficient of Young's modulus, and the quality factor decreased with temperature and remained unrecoverable upon cooling. The degradation in the quality factor is suspected to be a result of oxidation of the titatium nitride (TiN) top electrode, which increases the resistivity and leads to an unrecoverable reduction in the quality factor. The robust piezoelectric response of AlN at these temperatures show that AlN is a very promising candidate for elevated temperature applications.
Fluid inclusions are found within mineral crystals or along grain boundaries in many sedimentary rocks, notably in evaporite formations, and can migrate along a thermal or hydro-mechanical gradient. Shale and salt rocks have been considered potential host rocks for radioactive waste disposal, due to their low permeability. Previously stagnant inclusions may become mobilised by a perturbation of the in situ state by a geotechnical installation or the emplacement of heat-generating waste. The migration of fluid inclusions can thus have important impacts on the long-term performance of a geologic repository for high-level radioactive waste disposal. As a part of the international research project DECOVALEX-2019, two aspects of fluid inclusion migration in rock salt are currently investigated under different boundary conditions: a) altered hydro-mechanical conditions as a consequence of tunnel excavation or borehole drilling and b) coupled thermo-hydro-mechanical-chemical conditions during the heating period of the post-closure phase of a repository. To obtain a mechanistic understanding of underlying physical processes for fluid inclusion migration, a multi-scale modelling strategy has been developed. Microscale hydraulic and time-dependent mechanical conditions related to the creep behaviour of rock salt are constrained by considering the macroscale stress evolution of an underground excavation. An analysis using a coupled two-phase flow and elasto-plastic model with a consideration of permeability variation indicates that a pathway dilation along the halite grain boundary may increase the permeability by two orders of magnitude. The calculated high flow velocity may explain the fast pressure build-up observed in the field. In addition, a mathematical model for the migration and morphological evolution of a single fluid inclusion under a thermal gradient has been formulated. A first-order analysis of the model leads to a simple mathematical expression that is able to explain the key observations of thermally driven inclusion migration in salt. Finally, numerical methods such as a phase field method for solving a moving boundary problem of fluid inclusion migration have also been explored.
Gold-plated copper alloys are used extensively in electrical contacts where diffusional processes are known to cause contact degradation. An in situ transmission electron microscopy (TEM) heating study was carried out to provide fundamental understanding of the aging phenomena in reasonable timescales. Samples to visualize the interface in TEM were prepared by focused ion beam (FIB) microscopy and heated in situ up to 350°C while holding at intermediate temperatures to enable imaging. The grain boundaries in Au coatings, specifically the columnar boundaries, provided rapid pathways for diffusion of Cu all the way to the Au surface. This unequal diffusion created vacancies in Cu which coalesced into Kirkendall voids. This in situ technique has been applied to visualize the diffusion pathways in electroplated and sputtered Au films deposited directly on Cu, as well the role of Ni and NiP as barrier layers for mitigating Cu diffusion.
Microstructure reconstruction problems are usually limited to the representation with finitely many number of phases, e.g. binary and ternary. However, images of microstructure obtained through experimental, for example, using microscope, are often represented as a RGB or grayscale image. Because the phase-based representation is discrete, more rigid, and provides less flexibility in modeling the microstructure, as compared to RGB or grayscale image, there is a loss of information in the conversion. In this paper, a microstructure reconstruction method, which produces images at the fidelity of experimental microscopy, i.e. RGB or grayscale image, is proposed without introducing any physics-based microstructure descriptor. Furthermore, the image texture is preserved and the microstructure image is represented with continuous variables (as in RGB or grayscale images), instead of binary or categorical variables, which results in a high-fidelity image of microstructure reconstruction. The advantage of the proposed method is its quality of reconstruction, which can be applied to any other binary or multiphase 2D microstructure. The proposed method can be thought of as a subsampling approach to expand the microstructure dataset, while preserving its image texture. Moreover, the size of the reconstructed image is more flexible, compared to other machine learning microstructure reconstruction method, where the size must be fixed beforehand. In addition, the proposed method is capable of joining the microstructure images taken at different locations to reconstruct a larger microstructure image. A significant advantage of the proposed method is to remedy the data scarcity problem in materials science, where experimental data is scare and hard to obtain. The proposed method can also be applied to generate statistically equivalent microstructures, which has a strong implication in microstructure-related uncertainty quantification applications. The proposed microstructure reconstruction method is demonstrated with the UltraHigh Carbon Steel micrograph DataBase (UHCSDB).
Harden and optimize the ROCm based AMD GPU backend, develop a prototype backend for the Intel ECP Path Forward architecture, and improve the existing prototype Remote Memory Space capabilities.
The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (FCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling. These priorities are directly addressed in the SFWST Geologic Disposal Safety Assessment (GDSA) control account, which is charged with developing a geologic repository system modeling and analysis capability, and the associated software, GDSA Framework, for evaluating disposal system performance for nuclear waste in geologic media. GDSA Framework is supported by SFWST Campaign and its predecessor the Used Fuel Disposition (UFD) campaign.
We demonstrate a platform for phase and amplitude modulation in silicon nitride photonic integrated circuits via piezo-optomechanical coupling using tightly mechanically coupled aluminum nitride actuators. The platform, fabricated in a CMOS foundry, enables scalable active photonic integrated circuits for visible wavelengths, and the piezoelectric actuation functions without performance degradation down to cryogenic temperatures. As an example of the potential of the platform, we demonstrate a compact (∼40 µm diameter) silicon nitride ring resonator modulator operating at 780 nm with intrinsic quality factors in excess of 1.5 million, >10 dB change in extinction ratio with 2 V applied, a switching time less than 4 ns, and a switching energy of 0.5 pJ/bit. We characterize the exemplary device at room temperature and 7 K. At 7 K, the device obtains a resistance of approximately 20 teraohms, allowing it to operate with sub-picowatt electrical power dissipation. We further demonstrate a Mach-Zehnder modulator constructed in the same platform with piezoelectrically tunable phase shifting arms, with 750 ns switching time constant and 20 nW steady-state power dissipation at room temperature.
A mechanism to electrically tune the frequency of terahertz quantum cascade lasers (QCLs) is developed that allows for tuning, while the QCL is operated close to its peak bias and temperature. Two optically coupled but electrically isolated cavities are used in which the bias of a control cavity tunes the resonant-mode of the coupled QCL cavity independent of the QCL's operating bias. Approximately 4 GHz electrical tuning is realized for a 3.6 THz distributed-feedback QCL operating in pulsed mode at 58 K in a Stirling cooler. The single-mode QCL emits near-constant peak-power in the range of 5 - 5.3 mW through the tuning range and radiates in a narrow single-lobed beam with a far-field divergence of ∼ 4 ° × 11 °. The superlattice structure of the QCL is designed to implement a low-voltage intersubband absorption transition that is detuned from that of its gain transition, the strength of which could be controlled sensitively with applied voltage utilizing resonant-tunneling injection of electrons in the absorption subband. The tuning is realized by the application of small bias voltages (∼ 6 - 7 V) and requires a narrow bias range (∼ 1 V, ∼ 40 A / cm 2) to traverse across the entire tuning range, and the method should be generally applicable to all intersubband lasers including mid-infrared QCLs.
As demands for memory-intensive applications continue to grow, the memory capacity of each computing node is expected to grow at a similar pace. In high-performance computing (HPC) systems, the memory capacity per compute node is decided upon the most demanding application that would likely run on such system, and hence the average capacity per node in future HPC systems is expected to grow significantly. However, since HPC systems run many applications with different capacity demands, a large percentage of the overall memory capacity will likely be underutilized; memory modules can be thought of as private memory for its corresponding computing node. Thus, as HPC systems are moving towards the exascale era, a better utilization of memory is strongly desired. Moreover, upgrading memory system requires significant efforts. Fortunately, disaggregated memory systems promise better utilization by defining regions of global memory, typically referred to as memory blades, which can be accessed by all computing nodes in the system, thus achieving much better utilization. Disaggregated memory systems are expected to be built using dense, power-efficient memory technologies. Thus, emerging nonvolatile memories (NVMs) are placing themselves as the main building blocks for such systems. However, NVMs are slower than DRAM. Therefore, it is expected that each computing node would have a small local memory that is based on either HBM or DRAM, whereas a large shared NVM memory would be accessible by all nodes. Managing such system with global and local memory requires a novel hardware/software co-design to initiate page migration between global and local memory to maximize performance while enabling access to huge shared memory. In this paper we provide support to migrate pages, investigate such memory management aspects and the major system-level aspects that can affect design decisions in disaggregated NVM systems
This paper takes a critical look at the materials aspects of thermal runaway of lithium-ion batteries and correlates contributions from individual cell components to thermal runaway trends. An accelerating rate calorimeter (ARC) was used to evaluate commercial lithium-ion cells based on LiCoO2 (LCO), LiFePO4 (LFP), and LiNixCoyAl1-x-yO2 (NCA) at various states of charge (SOC). Cells were disassembled and the component properties were evaluated by thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and temperature-resolved X-ray diffraction (TR-XRD). The whole cell thermal runaway onset temperature decreases and peak heating rate increases with SOC due to cathode destabilization. LCO and NCA cathodes are metastable, with NCA cells exhibiting the highest thermal runaway rates. By contrast, the LFP cathode is stable to >500 °C, even when charged. For anodes, the decomposition and whole cell self-heating onset temperature is generally independent of SOC. DSC exotherm onset temperatures of the anodes were generally within 10 °C of the onset of self-heating in whole cell ARC. However, onset temperatures of the cathodes were typically observed above the ARC onset of whole cell runaway. This systematic evaluation of component to whole cell degradation provides a scientific basis for future thermal modeling and design of safer cells.
Industry seems to have the isolated-single-turbine-in-a-flat-field problem "solved". Wide-scale deployment and cost competitiveness require large wind plants — tens to hundreds of turbines in complex terrain. Predictive wind farm simulation is an exascale problem.
We demonstrate all-optical switching of high quality factor quasibound states in the continuum resonances in broken symmetry GaAs metasurfaces. By slightly breaking the symmetry of the GaAs nanoresonators, we enable leakage of symmetry protected bound states in the continuum (BICs) to free space that results in sharp spectral resonances with high quality factors of ∼500. We tune the resulting quasi-BIC resonances with ultrafast optical pumping at 800 nm and observe a 10 nm spectral blue shift of the resonance with pump fluences of less than 100 μJ cm-2. The spectral shift is achieved in an ultrafast time scale (<2.5 ps) and is caused by a shift in the refractive index mediated by the injection of free carriers into the GaAs resonators. An absolute reflectance change of 0.31 is measured with 150 μJ cm-2. Our results demonstrate a proof-of-concept that these broken symmetry metasurfaces can be modulated or switched at ultrafast switching speeds with higher contrast at low optical fluences (<100 μJ cm-2) than conventional Mie-metasurfaces.
This milestone shows a demonstration of the surface map model to map the surface temperature for single-phase computational fluid dynamics simulation data from STAR-CCM+ to the subchannel code CTF using a two-step process that captures global and local rod surface temperature behavior based on simulation boundary conditions and position. This model can be used to improve results from CTF by transforming surface temperatures obtained from CTF to model data that resembles data from STAR-CCM+. A summary of the current two-step model process and results of the initial investigations with a single subchannel and a 5 x 5 set of fuel rods are shown and discussed within, with suggestions for further improvements.
This report documents the completion of milestone STPM12-19 Documented Kokkos application usecases. The goal of this milestone was to develop use case examples for common patterns users implement with Kokkos. This work was performed in the fourth quarter of FY19 and resulted in use case descriptions available in the Kokkos Wiki, with code examples.
This report documents the completion of milestone STPRO4-25 Harden and optimize the ROCm based AMD GPU backend, develop a prototype backend for the Intel ECP Path Forward architecture, and improve the existing prototype Remote Memory Space capabilities. The ROCM code was hardened up to the point of passing all Kokkos unit tests - then AMD deprecated the programming model, forcing us to start over in FY20 with HIP. The Intel ECP Path Forward architecture prototype was developed with some initial capabilities on simulators - but plans changed, so that work will not continue. Instead SYCL will be developed as a backend for Aurora. Remote Spaces was improved. Development is ongoing part of a collaboration with NVIDIA.
This report documents the completion of milestone STPRO4-24 "Provide high quality (production) Kokkos support and consultation for ASC applications and libraries.". The Kokkos team resolved 344 issues reported to github (no explicit tracking of ASC vs non-ASC was performed). We engaged actively on the Kokkos slack channel, which now averages about 60 unique users per week. A survey of ASC customers was conducted with regards to the experienced support. The feedback indicates that the ASC customers are satisfied by Kokkos' support efforts.
This report documents the completion of milestone STPRO4-26 Engaging the C++ Committee. The Kokkos team attended the three C++ Committee meetings in San Diego, Hawaii, and Cologne with multiple members, updated multiple in-flight proposals (e.g. MDSpan, atomic ref), contributed to numerous proposals central for future capabilities in C++ (e.g. executors, affinity) and organized a new effort to introduce a Basic Linear Algebra library into the C++ standard. We also implemented a production quality version of mdspan as the basis for replacing the vast majority of the implementation of Kokkos::View, and thus start the transitioning of one of the core features in Kokkos to its future replacement.
Neuromorphic computers based on analogue neural networks aim to substantially lower computing power by reducing the need to shuttle data between memory and logic units. Artificial synapses containing nonvolatile analogue conductance states enable direct computation using memory elements; however, most nonvolatile analogue memories require high write voltages and large current densities and are accompanied by nonlinear and unpredictable weight updates. Here, we develop an inorganic redox transistor based on electrochemical lithium-ion insertion into LiXTiO2 that displays linear weight updates at both low current densities and low write voltages. The write voltage, as low as 200 mV at room temperature, is achieved by minimizing the open-circuit voltage and using a low-voltage diffusive memristor selector. We further show that the LiXTiO2 redox transistor can achieve an extremely sharp transistor subthreshold slope of just 40 mV/decade when operating in an electrochemically driven phase transformation regime.
How does a structure fail in an explosion? Is there an engineering design opportunity to mitigate the effects of a blast? These are questions that may be addressed with simulation; however, a computational simulation capability to answer these questions must include developments in material physics modeling and numerical algorithms that we believe are possible, but do not yet exist. This LDRD project proposes to develop the physics models and computational approaches to create a computational simulation approach to answer these physics and engineering questions.
An overview of the Uintah code and benchmark case is presented. Unitah provides a parallel, adaptive, multi-physics framework and solves time-dependent PDEs in parallel.
Sevougian, Stephen D.; Stein, Emily R.; Laforce, Tara; Perry, Frank V.; Nole, Michael; Chang, Kyung W.
The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy, Office of Spent Fuel and Waste Disposition (SFWD), has been conducting research and development on generic deep geologic disposal systems (i.e., geologic repositories). This report describes specific activities in the second half of Fiscal Year (FY) 2019 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package within the SFWST Campaign. The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and system performance assessment (PA) models in several host-rock environments, and to simulate and analyze these generic repository concepts and models using the GDSA Framework toolkit, and other tools as needed.
In this paper, we study, analyze, and validate some important zero-dimensional physics-based models for vanadium redox batch cell (VRBC) systems and formulate an adequate physics-based model that can predict the battery performance accurately. In the model formulation process, a systems approach to multiple parameters estimation has been conducted using VRBC systems at low C-rates (~C/30). In this batch cell system, the effect of ions' crossover through the membrane is dominant, and therefore, the capacity loss phenomena can be explicitly observed. Paradoxically, this means that using the batch system might be a better approach for identifying a more suitable model describing the effect of ions transport. Next, we propose an efficient systems approach, which enables to help understand the battery performance quickly by estimating all parameters of the battery system. Finally, open source codes, executable files, and experimental data are provided to enable people's access to robust and accurate models and optimizers. In battery simulations, different models and optimizers describing the same systems produce different values of the estimated parameters. Providing an open access platform can accelerate the process to arrive at robust models and optimizers by continuous modification from the users' side.
The fillers R&D program, mostly experimental, is part of a broader R&D program that includes new process modeling and performance assessment of criticality effects and the overall importance of criticality to repository performance (consequence screening). A literature research and consultation effort with experts by Hardin and Brady (2018) identified several potentially effective and workable filler materials including cements (primarily phosphate based), molten-metal alloys, and low-temperature glasses. Filler attributes were defined and the preliminary lists were compared qualitatively. Further comparative analysis will be done (e.g., cost estimates) after experimental screening has narrowed the list of alternatives. The following cement filler compositions were selected for experimental development work and accelerated testing in FY19: Aluminum phosphate cements (APCs); more specifically aluminum oxide / aluminum phosphate (Al2O3/ AlPO4) cements in which Al2O3 serves as the filler material bound by an AlPO4 binder formed by the reaction of Al2O3 with H3PO4; Calcium phosphate cements (CPCs); more specifically composed of pure or nearly pure hydroxyapatite (Ca5(PO4)3(OH)); Magnesium potassium phosphate cements (MKPs) composed of magnesium oxide / magnesium potassium phosphate (MgO / MgKPO4) cements in which MgO serves as the filler and MgKPO4 serves as the binder formed by the reaction of MgO with monopotassium phosphate (KH2PO4) and tricalcium phosphate ((Ca3(PO4)2); Two additional potential cement materials were explored preliminarily as the result of: (1) continued literature investigations into other filler candidates (wollastonite-based phosphate ceramic) and (2) the experimental discovery of a well-consolidated fly ash phosphate cement during the evaluation of fly ash as a potential filler material with Al2O3in APCs. Fly ash phosphate cements, more specifically in which a fly ash material composed primarily of mullite and quartz serves as the filler and is reacted with H3PO4 to form amorphous phosphate phase(s) as the binder; Wollastonite aluminum phosphate cements (WAPC), specifically wollastonite / aluminum phosphate (CaSiO3/ AlPO4) in which CaSiO3 serves as the filler material and AlPO4 serves as the binder formed by Al(OH)3 or metakaolin as Al sources and H3PO4 or ammonium dihydrogen phosphate (ADP) (NH4H2PO4) as phosphate sources. The FY19 effort focused on the optimization of compositions and subsequent processing of these five materials to achieve dense and well-consolidated monolithic samples with relatively low porosity. Once these goals were met basic material properties screening evaluations were performed including an assessment of dissolution resistance in water at elevated temperature (200 °C) and mechanical testing including unconfined compressive strength (UCS) testing. To date, the aluminum phosphate cements (APCs) appear to show the most promise for continued development. They are easily prepared and form smooth pourable slurries that remain stable for days with relatively low viscosities of several thousand centipoise (cP). They are then set at elevated temperatures (e.g., 170 °C) under ambient (0.1 MPa) or elevated pressure (~1MPa). Overall, they demonstrate the best dissolution resistance in water at elevated temperature (200 °C) and good compressive strengths. However, additional effort is required to optimize the APC slurry formulations and the process used for thermal curing these materials. The calcium phosphate cements (CPCs) can be formed at room temperature to produce a well-consolidated body. However, their slurry viscosities are very high (and difficult to measure) and they exhibit relatively short cure times of 2 to 3 hours. Also, dissolution resistance is very poor, the poorest of all the cements examined The same is the case for the small number of MKP cements fabricated; they cure very quickly (10 minutes or less) and disintegrate within a few hours upon immersion in distilled water. Surprisingly, fly ash reacts with phosphoric acid to form dense and well-consolidated cements but the mixture rapidly sets at room temperature (less than 30 minutes) and the subsequent conversion of the binder to an amorphous phosphate phase(s) as a function of temperature is complicated. Finally, the wollastonite aluminum phosphate cements (WAPC) are easily prepared and form smooth pourable slurries that remain stable for several hours. They are then set at 130 °C. A WAPC sample exhibited the highest compressive strengths of all the materials we evaluated but in general their dissolution resistance to water is poor.
Keto-hydroperoxides (KHPs) are reactive, partially oxidized intermediates that play a central role in chain-branching reactions during the gas-phase low-temperature oxidation of hydrocarbons and oxygenated species. Although multiple isomeric forms of the KHP intermediate are possible in complex oxidation environments when multiple reactant radicals exist that contain nonequivalent O2 addition sites, isomer-resolved data of KHPs have not been reported. In this work, we provide partially isomer-resolved detection and quantification of the KHPs that form during the low-temperature oxidation of tetrahydrofuran (THF, cycl.-O-CH2CH2CH2CH2-). We describe how these short-lived KHPs were detected, identified, and quantified using integrated experimental and theoretical approaches. The experimental approaches were based on direct molecular-beam sampling from a jet-stirred reactor operated at near-atmospheric pressure and at temperatures between 500 and 700 K, followed by mass spectrometry with single-photon ionization via tunable synchrotron-generated vacuum-ultraviolet radiation, and the identification of fragmentation patterns. The interpretation of the experiments was guided by theoretical calculations of ionization thresholds, fragment appearance energies, and photoionization cross sections. On the basis of the experimentally observed and theoretically calculated ionization and fragment appearance energies, KHP isomers could be distinguished as originating from H-abstraction reactions from either the α-C adjacent to the O atom or the β-C atoms. Temperature-dependent concentration profiles of the partially resolved isomeric KHP intermediates were determined in the range of 500-700 K, and the results indicate that the observed KHP isomers are formed overwhelmingly (∼99%) from the α-C THF radical. Comparisons of the partially isomer-resolved quantification of the KHPs to up-to-date kinetic modeling results reveal new opportunities for the development of a next-generation THF oxidation mechanism.
To ensure reliable and predictable service in the electrical grid between distributed renewable distributed energy resources (DERs) it is important to gauge the level of trust present within critical components and DER aggregators (DERAs). Although trust throughout a smart grid is temporal and dynamically varies according to measured states, it is possible to accurately formulate communications and service level strategies based on such trust measurements. Utilizing an effective set of machine learning and statistical methods, it is shown that establishment of trust levels between DERAs using behavioral pattern analysis is possible. Further, it is also shown that the establishment of such trust can facilitate simple secure communications routing between DERAs. Providing secure routing between DERAs enables a grid operator to maintain service level agreements to its customers, reduce the attack surface and increase operational resiliency.
King, Michael P.; Ryder, Kaitlyn L.; Ryder, Landen D.; Sternberg, Andrew L.; Kozub, John A.; Zhang, En X.; Khachatrian, Ani; Buchner, Steven P.; Mcmorrow, Dale P.; Hales, Joel M.; Zhao, Yuanfu; Wang, Liang; Wang, Chuanmin; Weller, Robert A.; Schrimpf, Ronald D.; Weiss, Sharon M.; Reed, Robert A.; Black, Dolores A.
A sensitive volume is developed using pulsed laser-induced collected charge for two bias conditions in an epitaxial silicon diode. These sensitive volumes show good agreement with experimental two photon absorption laser-induced collected charge at a variety of focal positions and pulse energies. When compared to ion-induced collected charge, the laser-based sensitive volume over predicts the experimental collected charge at low bias and agrees at high bias. Here, a sensitive volume based on ion-induced collected charge adequately describes the ion experimental results at both biases. Differences in the amount of potential modulation explain the differences between the ion-and laser-based sensitive volumes at the lower bias. Truncation of potential modulation by the highly doped substrate at the higher bias results in similar sensitive volumes.
The relationship between the structure and thermodynamic properties of schoepite, an important uranyl phase with formula [(UO2)8O2(OH)12]·12H2O formed upon corrosion of UO2, has been investigated within the framework of density functional perturbation theory (DFPT). Experimental crystallographic lattice parameters are well reproduced in this study using standard DFT. Phonon calculations within the quasi-harmonic approximation predict standard molar entropy and isobaric heat capacity of S0 = 179.60 J mol-1 K-1 and C0P = 157.4 J mol-1 K-1 at 298.15 K, i.e., ~6% and ~4% larger than existing DFPT-D2 calculations. The computed variation of the standard molar isobaric heat capacity with water content from schoepite (UO3·xH2O, x = 2.25) to dehydrated schoepite (x = 1) is predicted to be essentially linear along isotherms ranging from 100 to 500 K. Finally, these findings have important implications for the dehydration of layered uranyl corrosion phases and hygroscopic materials.
There have been recent efforts towards the development of biologically-inspired neuromorphic devices and architecture. Here, we show a synapse circuit that is designed to perform spike-timing-dependent plasticity which works with the leaky, integrate, and fire neuron in a neuromorphic computing architecture. The circuit consists of a three-terminal magnetic tunnel junction with a mobile domain wall between two low-pass filters and has been modeled in SPICE. The results show that the current flowing through the synapse is highly correlated to the timing delay between the pre-synaptic and post-synaptic neurons. Using micromagnetic simulations, we show that introducing notches along the length of the domain wall track pins the domain wall at each successive notch to properly respond to the timing between the input and output current pulses of the circuit, producing a multi-state resistance representing synaptic weights. We show in SPICE that a notch-free ideal magnetic device also shows spike-timing dependent plasticity in response to the circuit current. This work is key progress towards making more bio-realistic artificial synapses with multiple weights, which can be trained online with a promise of CMOS compatibility and energy efficiency.
Tensile properties, fatigue crack initiation, fatigue crack growth rate, and fatigue life are evaluated in 304L austenitic stainless steel fabricated by directed energy deposition (DED). Large lack of fusion (LoF) defects (often >1 mm in length) significantly reduce ultimate tensile strength and ductility, as well as accelerate fatigue crack initiation and reduce fatigue life. In comparison, small spherical defects (<100 μm in diameter) have less effect on tensile and fatigue properties. Fatigue crack growth rate is less severely affected by defects than other properties, showing only local acceleration in the proximity of LoF defects. Therefore, shorter fatigue life is attributed to the role of LoF defects on facilitating fatigue crack initiation and to a lesser extent fatigue crack propagation. Additionally, the fatigue life can be normalized for defects by considering their effect on ultimate tensile strength, suggesting that in the limit of low defect population, the fatigue strength of additively manufactured stainless steel is similar to conventional wrought materials.
SNAP potentials are inter-atomic potentials for molecular dynamics that enable simulations at accuracy levels comparable to density functional theory(DFT) at a fraction of the cost. As such, SNAP scales to on the order of 104 — 106 atoms. In this work, we explore CPU optimization of potentials computation using SIMD. We note that efficient use of SIMD is non-obvious as the application features an irregular iteration space for various potential terms, necessitating use of SIMD across atoms in a cross matrix, batched fashion. We present a preliminary analytic model to determine the correct batch size for several CPU architectures across several vendors, and show end-to-end speedups between 1.66x and 3.22x compared to the original.
The aim and scope of this project was the development of a capability to prepare high-quality, epitaxial beta gallium oxide films by oxide reactive molecular-beam epitaxy. The purpose was to demonstrate that beta gallium oxide could be grown by such a method using Sandia’s existing oxide molecular-beam epitaxy instrument. The key activity in this project was the installation of a gallium oxide capability on the Sandia instrument. This required the acquisition of several custom items for the instrument, including: a gallium effusion cell, appropriate cell power supplies and temperature controllers, a shutter to block beam flux, installation of an existing ozone generator with a directed gas nozzle and controlled leak valve, and re-routing the chilled water system to accommodate the cell. In addition, beta gallium oxide single crystals were acquired and their surfaces characterized by reflection high energy electron diffraction.
Diffraction series data have been acquired and analyzed via multivariate statistical analysis. For two different data series analyzed, the data analysis was able to reduce the raw diffraction data series into a much smaller easier-to-interpret solution consisting mainly of crystallographic phase and orientation information.
EmrE is a small, homodimeric membrane transporter that exploits the established pH gradient across the E. coli inner membrane to export polyaromatic cations that might otherwise inhibit cellular growth. While herculean efforts through experimental studies have established many fundamental facts about the specificity and rate of substrate transport in EmrE, the low resolution of the available structures have hampered efforts to tie those findings to the EmrE coupling mechanism between proton and small molecule substrates. Here we present a full three-dimensional structure of EmrE optimized against available cyro-EM data to delineate the critical interactions by which EmrE regulates its conformation. We use the generated structural model to conduct equilibrium and nonequilibrium molecular dynamics simulations to probe EmrE dynamics under different substrate loading states, representing different states in the transport cycle. The model is stable under extended simulation, and reveals that water dynamics within the EmrE lumen change substantially with the loading state. The water dynamics cause hydrogen bonding networks to shift radically when the protonation states change for a pair of solvent-exposed glutamate residues (E14) within the lumen of the transporter, which are proposed to act as proton binding sites during the transport cycle. One specific hydrogen bond from a tyrosine (Y60) of one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. Furthermore, the hydrogen bond provided by Y60 lowers the pKa of the interacting glutamate relative to its partner on the opposite monomer such that it will protonate second, establishing the need for both glutamates to be protonated for the hydrogen bond to break and a substrate-free transition to take place.
Thin films are materials systems that are widely used in applications ranging from electronics and optical devices to industrial and biomedical ones. However, these systems are unstable against various homogenization processes and aging mechanisms (grain growth, coarsening, surface evolution, and diffusion of species) even at low service temperatures. In this work, we examine the role of various aspects of microstructure (grain boundary types and characters, free surfaces, surface diffusion, and thermal grooves) on the thermal aging of such systems. Existing experimental tools will be leveraged to characterize thin films, i.e., in-situ quantitative thermal annealing, and provide direct comparisons to predications emerging from a recently developed meso-scale model via Precession Electron Diffraction (PED). Parametric studies will be conducted to gain insights on the role of each of the aforementioned aspects of microstructure on the dynamics. Herein, we will focus on “hard” gold thin films (with alloying elements like Co, Ni, or Fe) as they are materials of choice in a wide range of applications at Sandia. Initially, pure gold systems are examined and future studies will focus on microstructure evolution in the presence of alloying elements and second-phase particles.
Density functional theory (DFT) is undergoing a shift from a descriptive to a predictive tool in the field of solid state physics, heralded by a spike in “high-throughput” studies. However, methods to rigorously evaluate the validity and accuracy of these studies is lacking, raising serious questions when simulation and experiment disagree. In response, we have developed the V-DM/16 test set, designed to evaluate the experimental accuracy of DFT’s various implementations for pe riodic transition metal solids. Our test set evaluates 26 transition metal elements and 80 transition metal alloys across three physical observables: lattice constants, elastic coefficients, and formation energy of alloys. Whether or not a functional can accurately evaluate the formation energy offers key insights into whether the relevant physics are being captured in a simulation, an especially impor tant question in transition metals where active d-electrons can thwart the accuracy of an otherwise well-performing functional. Our test set captures a wide variety of cases where the unique physics present in transition metal binaries can undermine the effectiveness of “traditional” functionals. By application of the V/DM-16 test set, we aim to better characterize the performance of existing functionals on transition metals, and to offer a new tool to rigorously evaluate the performance of new functionals in the future.
Symmetric four-point bending (S4PB) and anti-symmetric four-point bending (AS4PB) were applied to assess the effect of hydrogen on crack initiation and propagation in a stable, nitrogen-strengthened austenitic stainless steel (21Cr-6Ni-9Mn). Specimens of a high strength aluminum alloy (AA2219-T851), which has also been identified as hydrogen-compatible, were used for test method development prior to completing the stainless steel test matrix. Single edge notched bend (SEN(B)) specimens were extracted from forged 21Cr-6Ni-9Mn bar and AA2219 plate. The 21Cr-6Ni-9Mn specimens were then hydrogen charged with ~200 wt. ppm hydrogen; AA2219 specimens were not charged. Aluminum specimens were tested in S4PB to induced mode I (pure bending), and AS4PB to attain varying levels of mode I/II mixity and mode II (pure shear), where the ratio of mode I to mode II varies with the position of the crack plane relative to the load line. After test method troubleshooting and validation, hydrogen charged stainless steel specimens were then subjected to mode I and two ratios of mixed mode I/II. Mode II loading was not achieved due to high load limitations. Analyses of fracture profiles for both materials reveal a marked effect of loading mode mixity on initial crack propagation orientation, however a specific contribution of hydrogen was not readily identifiable. Fracture initiation toughness in the presence of hydrogen, JIH, was calculated following the J-integral approach; the mode I JIH calculated for stainless steel samples fractured in S4PB were consistent with published values determined from compact tension specimens. The peak loads required to initiate fracture during AS4PB far exceeded those required during S4PB. This is attributed to the increasing shear force applied to the specimens as the degree of mode mixity increases. Ultimately, understanding the fracture response of hydrogen exposed stainless steels subjected to mixed mode I/II loading is critical for designing hydrogen containment vessels or gas transfer systems (GTS).
The purpose of this report is to review technical issues relevant to the performance evaluation of dry storage systems during vacuum drying and long-term storage operations. It also provides updates on experimental components under development that are vital for pursuing advanced studies. Validation of the extent of water removal in a multi-assembly dry storage system using an industrial vacuum drying procedure is needed, as operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system. Water remaining in canisters/casks upon completion of vacuum drying can lead to cladding corrosion, embrittlement, and breaching, as well as fuel degradation. Therefore, additional information is needed to evaluate the potential impacts of water retention on extended long-term dry storage. A general lack of data and experience modeling the drying process necessitates the testing of advanced concepts focused on the simulation of industrial vacuum drying. Smaller-scale tests that incorporate relevant physics and well-controlled boundary conditions are necessary to provide insight and guidance to the modeling of prototypic systems undergoing drying processes. This report describes the development and testing of waterproof, electrically-heated spent fuel rod simulators as a proof of concept to enable experimental simulation of the entire dewatering and drying process. This report also describes the preliminary development of specially-designed, unheated mock fuel rods for monitoring internal rod pressures and studying water removal from simulated failed fuel rods. A variety of moisture monitoring instrumentation is also being considered and will be downselected for the tracking of dewpoints of gas samples. The effects of cladding oxidation and crud on water retention in dry storage systems can be explored via separate effects tests (SETs) that would measure chemisorbed and physisorbed water content on cladding samples. The concepts listed above will be incorporated into an advanced dry cask simulator with multiple fuel assemblies in order to account for important inter-assembly heat-transfer physics. Plans are described for harvesting up to five full-length 5x5 laterally truncated assemblies from commercial 17x17 PWR skeleton components with the goal of constructing this simulator.
Computer Methods in Applied Mechanics and Engineering
Carlberg, Kevin T.; Guzzetti, Sofia; Alvarez, L.A.; Blanco, P.J.; Veneziani, A.
Numerical simulations of the cardiovascular system are affected by uncertainties arising from a substantial lack of data related to the boundary conditions and the physical parameters of the mathematical models. Quantifying the impact of this uncertainty on the numerical results along the circulatory network is challenged by the complexity of both the morphology of the domain and the local dynamics. Here, we propose to integrate (i) the Transverse Enriched Pipe Element Methods (TEPEM) as a reduced-order model for effectively computing the 3D local hemodynamics; and (ii) a combination of uncertainty quantification via Polynomial Chaos Expansion and classical relaxation methods – called network uncertainty quantification (NetUQ) – for effectively propagating random variables that encode uncertainties throughout the networks. The findings demonstrate the computational effectiveness of computing the propagation of uncertainties in networks with nontrivial topology, including portions of the cerebral and the coronary systems.
Magnetically driven experiments supporting pulsed-power utilize a wide range of configurations, including wire-arrays, gas-puffs, flyer plates, and cylindrical liners. This experimental flexibility is critical to supporting radiation effects, dynamic materials, magneto-inertial-fusion (MIF), and basic high energy density laboratory physics (HEDP) efforts. Ultimately, the rate at which these efforts progress is limited by our understanding of the complex plasma physics of these systems. Our effort has been to begin to develop an advanced algorithmic structure and a R&D code implementation for a plasma physics simulation capability based on the five-moment multi-fluid / full-Maxwell plasma model. This model can be used for inclusion of multiple fluid species (e.g., electrons, multiple charge state ions, and neutrals) and allows for generalized collisional interactions between species, models for ionization/recombination, magnetized Braginskii collisional transport, dissipative effects, and can be readily extended to incorporate radiation transport physics. In the context of pulsed-power simulations this advanced model will help to allow SNL to computationally simulate the dense continuum regions of the physical load (e.g. liner implosions, flyer plates) as well as partial power-flow losses in the final gap region of the inner MITL. In this report we briefly summarize results of applying a preliminary version of this model in the context of verification type problems, and some initial magnetic implosion relevant prototype problems. The MIF relevant prototype problems include results from fully-implicit / implicit-explicit (IMEX) resistive MHD as well as full multifluid EM plasma formulations.
This report has been accepted for publication in the journal Fusion Science and Technology, in the special issue associated with the Tritium 2016 conference, where the work was presented. Scanning calorimetry of a confined, reversible hydrogen sorbent material has been previously proposed as a method to determine compositions of unknown mixtures of diatomic hydrogen isotopologues and helium. Application of this concept could result in greater process knowledge during the handling of these gases. Previously published studies have focused on mixtures that do not include tritium. This paper focuses on modeling to predict the effect of tritium in mixtures of the isotopologues on a calorimetry scan. The model predicts that tritium can be measured with a sensitivity comparable to that observed for hydrogen-deuterium mixtures, and that under some conditions, it may be possible to determine the atomic fractions of all three isotopes in a gas mixture.
The linear energy-momentum dispersion which arises from graphene’s underlying honeycomb lattice gives graphene its unique electronic properties unfound in conventional semiconductors. Theoretically speaking, when an electrostatic potential with hexagonal or honeycomb symmetry is imposed onto a two-dimensional electron/hole system, the band structure is modified in a way that the same linear energy-momentum dispersion could exist. Experimentally, there has not been any evidence from transport demonstrating the so-called “artificial graphene”. In this project, we attempt to create an artificial superlattice potential with hexagonal symmetry for two dimensional carriers in an undoped SiGe heterostructure by patterning a nanoscale hole array in a metallic gate. Using undoped heterostructures allows us to access a very wide density range, which covers the magic densities at which the Dirac points are expected. A process flow for fabricating such field-effect-transistor devices with a lattice constant as small as 90 nm is reported. Magneto-transport measurements performed at 0.3 K show that the superlattice potential in the quantum well in which the two-dimensional system resides is indeed modulated by the gates. However, no signature of the sought-after linear dispersion is observed in the transport data.
At Sandia National Laboratories, I spent the summer of 2016 interning at the Combustion Research Facility. As the title suggests, the main goal of the facility is to investigate combustion with the goal of optimizing engine efficiency and reducing pollutants and emissions. At the lab, I worked in the field of spectroscopy, lasers and optics. Specifically, I was involved in a research group that investigated innovative methods of probing matter and chemical reactions with laser light.
Combined picosecond (10-12 s) and femtosecond (10-15 s) laser pulses can give sensitive, low-noise measurements of important quantities in reacting flows, such as species concentrations and temperature. Emmanuel’s work focused on the development of an instrument for tailoring the time profile of picosecond laser pulses for use in nonlinear optical spectroscopic methods created from broad bandwidth femtosecond pulses. In addition to constructing the device, Emmanuel produced a LabView-based automation code, building off skills he developed in a previous CCI internship at Sandia.
Radiation-imaging devices are important to the nuclear security and safeguards communities due to their ability to detect and localize radioactive sources. Because of their low natural background, difficulty to shield, and unique association with SNM, fast-neutron imaging provides a promising means for the detection of SNM. To make these R&D systems useful for end-user applications, robust reconstruction and analysis algorithms must be developed that provide detailed information on the location, energy spectrum, and intervening material. To date, deconvolution, match filtering, and MLEM have been demonstrated. However, the information provided is limited to an approximate location of the source.
In this report Chenhall presents a new Monte Carlo (MC) transport method developed as part of his PhD thesis work funded via fellowship by Sandia National Laboratories. The model is used to simulate the electron thermal transport within inertial confinement fusion (ICF) type problems. The new model aims to improve upon the currently used implicit Schurtz, Nicolai, and Busquet (iSNB) method, in particular by using finite particle ranges in comparison to the exponential solution of a diffusion method and by improved higher order angular modeling. The new method is being developed using the 1D LILAC and 2D DRACO codes developed by the Rochester Laboratory for Laser Energetics. The model is compared to iSNB for several ICF type problems: Omega shot 60303 a shock timing experiment, Omega shot 68951 a cryo target implosion and a NIF phase plate polar direct drive design. Overall, the MC method performs at least as well as the iSNB method and appears to indicate a lower predicted preheat ahead of the shock fronts. However, the difference in overall results between iSNB and MC transport is small for the choice of test problems and further simulations will need to be run to fully understand the difference between the models.
An initial foray into the design of specimens that can be used to provide data about the quasistatic ductile failure of metals when subjected to shear-dominated (low triaxiality) states of stress was undertaken. Four specimen geometries made from two materials with different ductility (Al 7075, lower ductility and steel A286, higher ductility) were considered as candidates. Based on results from analysis and experimentation, it seems that two show promise for further consideration. Whereas preliminary results indicate that the Johnson-Cook model fit the failure data for Al 7075 well, it did not fit the data for steel A286. Further work is needed to consolidate the results and evaluate other failure models that may fit the steel data better, as well as to extend the results of this work to the dynamic loading regime.
Concern over Arctic methane (CH4) emissions has increased following recent discoveries of poorly understood sources and predictions that methane emissions from known sources will grow as Arctic temperatures increase. New efforts are required to detect increases and explain sources without being confounded by the multiple sources. Methods for distinguishing different sources are critical. We conducted measurements of atmospheric methane and source tracers and performed baseline global atmospheric modeling to begin assessing the climate impact of changes in atmospheric methane. The goal of this project was to address uncertainties in Arctic methane sources and their potential impact on climate by (1) deploying newly developed trace-gas analyzers for measurements of methane, methane isotopologues, ethane, and other tracers of methane sources in the Barrow, AK, (2) characterizing methane sources using high-resolution atmospheric chemical transport models and tracer measurements, and (3) modeling Arctic climate using the state-of-the-art high- resolution Spectral Element Community Atmosphere Model (CAM-SE).
In experiments conducted on the Z-machine at Sandia National Laboratories, dynamic material properties cannot be analyzed using traditional analytic methods, necessitating solving an inverse problem. Bayesian model calibration is a statistical framework for solving an inverse problem to estimate parameters input into a computational model in the presence of multiple uncertainties. Disentangling input parameter uncertainty and model misspecification is often poorly identified problem. When using computational models for physical parameter estimation, the issue of parameter identifiability must be carefully considered to obtain accurate and precise estimates of physical parameters. Additionally, in dynamic material properties applications, the experimental output is a function, velocity over time. While we can sample an arbitrarily large number of points from the measured velocity, these curves only contain a finite amount of information about the calibration parameters. In this report, we propose modifications to the Bayesian model calibration framework to simplify and improve the estimation of physical parameters with functional outputs. Specifically, we propose scaling the likelihood function by an effective sample size rather than modeling the discrepancy function; and modularizing input nuisance parameters with weakly identified parameters. We evaluate the performance of these proposed methods using a statistical simulation study and then apply these methods to estimate parameters of the tantalum equation of state. We conclude that these proposed methods can provide simple, fast, and statistically valid alternatives to the full Bayesian model calibration procedure; and that these methods can be used to estimate parameters of the equation of state for tantalum.
A 150 lbf thrust class, modular, bi-propellant, rocket engine/gas-generator and supporting test infrastructure has been developed in a cooperative effort between Sandia National Laboratories and the New Mexico Institute of Mining and Technology’s (NMIMT’s) Energetic Materials Research and Testing Center (EMRTC). This modular test engine design consists of a head end fuel-oxidizer injector, a spark ignition gaseous H2/O2 torch igniter, combustion chamber and nozzle module. This robust design allows for rapid configuration changes as well as economical repair should hardware become damaged in testing. The engine interfaces with a permanently installed pressurizing system capable of delivering liquid nitrous oxide and a variety of liquid fuels for both rocket engine development and propellant performance evaluation. The regulated high pressure systems allow for delivery of liquefied gases above their saturation pressure as well as allowing for high pressure rocket engine/gas-generator operation. The facility test cell houses a 1 ton thrust capacity test stand leaving room for larger scale engine development.
Color centers (defect complexes such as SiV) in diamond have shown potential in fields ranging from metrology, cybersecurity to quantum computation. Demonstrations in these fields have pushed the envelope of state-of-the-art operations – for example, single photon sources (SPS) making use of SiV centers in diamond for quantum key distribution have demonstrated all the requirements for SPS operation including: (1) stable operation with second correlation function <<1, (2) electrically driven single photon emission and (3) compatibility with frequenc y conversion to telecommunication frequencies. To-date, however, all these demonstrations have been on lab-scale one-off devices. The key question behind how to deterministically fabricate these devices, namely activation yield has been overlooked. For context, Si based semiconductor devices are hugely successful because we have a high activation yield for implanted dopants. This is not yet true for diamond color centers. As currently understood, the color center yield is dominated by a lack of vacancies in the immediate area of the implantation. We propose to optimize the activation yield of color center using a combination of (1) focused single ion implantation with in-situ detection to count the number of implanted Si ions and (2) localized point defect (vacancy) creation using a focused Li ion beam to improve the yield. These experiments build on the unique capabilities of the SNL nanoImplanter (nI) to produce focused ion beam with spatial resolution of < 10 nm of both Si and Li ions. This work will also leverage our world-leading single ion implantation and detection capabilities.
The effort to develop larger-scale computing systems introduces a set of related challenges: Large machines are more difficult to synchronize. The sheer quantity of hardware introduces more opportunities for errors. New approaches to hardware, such as low-energy or neuromorphic devices are not directly programmable by traditional methods.
Schunk, Randy; Tjiptowidjojo, Kristianto; Cochrane, Andrew; Bonnecaze, Roger T.
The inextensible cylindrical shell theory and lubrication theory combine into a model for the elastohydrodynamics of a rolling-imprint modality of nanoimprint lithography (NIL). Foil-bearing theory describes the formation of the lubrication gap due to relative motion between a tensioned substrate and a rigid, cylindrical surface. Reproduction of the results of foil-bearing theory for both stiff and perfectly flexible substrates validates this coupled model and reveals a highly predictable region of uniformity that provides low shear stress conditions ideal for UV-cure. These results show theoretical limitations that are used to construct an operating window for predicting rolling-mode NIL process feasibility.
Efficiency bottlenecks inherent to conventional computing in executing neural algorithms have spurred the development of novel devices capable of “in-memory” computing. Commonly known as “memristors,” a variety of device concepts including conducting bridge, vacancy filament, phase change, and other types have been proposed as promising elements in artificial neural networks for executing inference and learning algorithms. In this article, we review the recent advances in memristor technology for neuromorphic computing and discuss strategies for addressing the most significant performance challenges, including nonlinearity, high read/write currents, and endurance. As an alternative to two-terminal memristors, we introduce the three-terminal electrochemical memory based on the redox transistor (RT), which uses a gate to tune the redox state of the channel. Decoupling the “read” and “write” operations using a third terminal and storage of information as a charge-compensated redox reaction in the bulk of the transistor enables high-density information storage. These properties enable low-energy operation without compromising analog performance and nonvolatility. Finally, we discuss the RT operating mechanisms using organic and inorganic materials, approaches for array integration, and prospects for achieving the device density and switching speeds necessary to make electrochemical memory competitive with established digital technology.
Engineered solid-liquid interfaces will play an important role in the development of future energy storage and conversion (ESC) devices. In the present study, defective graphene oxide (GO) and reduced graphene oxide (rGO) structures were used as engineered interfaces to tune the selectivity and activity of Pt disk electrodes. GO was deposited on Pt electrodes via the Langmuir-Blodgett technique, which provided compact and uniform GO films, and these films were subsequently converted to rGO by thermal reduction. Electrochemical measurements revealed that both GO and rGO interfaces on Pt electrodes exhibit selectivity toward the oxygen reduction reaction (ORR), but they do not have an impact on the activity of the hydrogen oxidation reaction in acidic environments. Scanning transmission electron microscopy at atomic resolution, along with Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), and scanning electron microscopy (SEM), revealed possible diffusion sites for H2 and O2 gas molecules and functional groups relevant to the selectivity and activity of these surfaces. Based on these insights, rGO interfaces are further demonstrated to exhibit enhanced activity for the ORR in nonaqueous environments and demonstrate the power of our ex situ engineering approach for the development of next-generation ESC devices.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.
In this work, a novel computational procedure, based on the principles of flat-histogram Monte Carlo, is developed for facile prediction of the adsorption thermodynamics of intrinsically flexible adsorbents. We then demonstrate how an accurate prediction of methane deliverable capacity in a metal–organic framework (MOF) with significant intrinsic flexibility requires use of such a method. Dynamic side chains in the framework respond to methane adsorbates and reorganize to exhibit a more conducive pore space at high adsorbate densities while simultaneously providing a less conducive pore space at low adsorbate densities. This “responsive pore” MOF achieves ~20% higher deliverable capacity than if the framework were rigid and elucidates a strategy for designing high deliverable capacity MOFs in the future.