Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment
The dissolution and depletion of chromium (Cr) in salt facing nickel (Ni) alloy surfaces is one of the predominant degradation mechanisms of structural components in molten salt technology. In this work, we use density functional theory to investigate the role of electronic level interactions that may underlie the depletion phenomenon of Cr in a Ni 100 surface exposed to various adsorbed salt species. Our results show that, under vacuum, Ni preferentially segregates to the surface layer. Conversely, in the presence of adsorbed anionic salt species (e.g., chlorine (Cl), fluorine (F) or the impurity oxygen (O)) Cr segregation becomes more favorable. In these cases, Cl has the weakest effect on segregation, while O has the strongest effect. Our analysis reveals the strong correlation between Cr segregation and the amount of valence charge transferred between the Cr atom and surface adsorbate: the greater the charge transfer, the lower the segregation energy. We also show that, when considered, secondary cations screen Cr-anion interactions, which in turn reduce the magnitude of the anions effect on segregation. These results shed light on the role of salt impurities likely play in the overall corrosion phenomena in molten salt environments. This work provides insights into the atomic level interactions fundamental to molten salt corrosion and on the importance of maintaining salt purity.
This SNL document contains requested radiological survey information, as part of the documentation for the MLU shipment being performed by the LANL MLU team. The survey was performed in TA-5, on October 6th, 2021. This survey was of two WIPP trailers carrying three empty TRUPACTs each.
Studying the mechanical behavior of silicon cell fractures is critical for understanding changes in PV module performance. Traditional methods of detecting cell cracks, e.g., electroluminescence (EL) imaging, utilize electrical changes and defects associated with cell fracture. Therefore, these methods reveal crack locations, but do not operate at the time or length scales required to accurately measure other physical properties of cracks, such as separation width and behavior under dynamic loads.
The adoption of digital technology into Instrumentation and Control (I&C) systems in nuclear facilities fundamentally changes the nature of these systems. Greater interconnectivity of reprogrammable, and functionally interdependent control systems has given rise to the need for computer security consideration in digital I&C Systems. The cyber security of I&C systems presents a growing risk to nuclear facilities and requires the development of educational and research tools to ensure the safety of these facilities. Currently there is a major gap in formal educational offerings on cyber security for these Operational Technology (OT) systems. To provide formal cyber security education resources, DOE’s office of International Nuclear Security (INS) partnered with the University of São Paulo (USP) to develop a training course on the cyber security of nuclear facility I&C systems using the hypothetical Nuclear Power Plant, Asherah.
We demonstrate the ability to fabricate vertically stacked Si quantum dots (QDs) within SiGe nanowires with QD diameters down to 2 nm. These QDs are formed during high-temperature dry oxidation of Si/SiGe heterostructure pillars, during which Ge diffuses along the pillars' sidewalls and encapsulates the Si layers. Continued oxidation results in QDs with sizes dependent on oxidation time. The formation of a Ge-rich shell that encapsulates the Si QDs is observed, a configuration which is confirmed to be thermodynamically favorable with molecular dynamics and density functional theory. The type-II band alignment of the Si dot/SiGe pillar suggests that charge trapping on the Si QDs is possible, and electron energy loss spectra show that a conduction band offset of at least 200 meV is maintained for even the smallest Si QDs. Our approach is compatible with current Si-based manufacturing processes, offering a new avenue for realizing Si QD devices.
The feasibility and component cost of hydrogen rail refueling infrastructure is examined. Example reference stations can inform future studies on components and systems specifically for hydrogen rail refueling facilities. All of the 5 designs considered assumed the bulk storage of liquid hydrogen on-site, from which either gaseous or liquid hydrogen would be dispensed. The first design was estimated to refuel 10 multiple unit trains per day, each train containing 260 kg of gaseous hydrogen at 350 bar on-board. The second base design targeted the refueling of 50 passenger locomotives, each with 400 kg of gaseous hydrogen on-board at 350 bar. Variations from this basic design were made to consider the effect of two different filling times, two different hydrogen compression methods, and two different station design approaches. For each design variation, components were sized, approximate costs were estimated for major components, and physical layouts were created. For both gaseous hydrogen-dispensing base designs, the design of direct-fill using a cryopump design was the lowest cost due to the high cost of the cascade storage system and gas compressor. The last three base designs all assumed that liquid hydrogen was dispensed into tender cars for freight locomotives that required 7,500 kg of liquid hydrogen, and the three different designs assumed that 5, 50, or 200 tender cars were refueled every day. The total component costs are very different for each design, because each design has a very different dispensing capacity. The total component cost for these three designs are driven by the cost of the liquid hydrogen tank; additionally, delivering that much liquid hydrogen to the refueling facility may not be practical. Many of the designs needed the use of multiple evaporators, compressors, and cryopumps operating in parallel to meet required flow rates. In the future, the components identified here can be improved and scaled-up to better fit the needs of heavy-duty refueling facilities. This study provides basic feasibility and first-order design guidance for hydrogen refueling facilities serving emerging rail applications.
Martinez-Tossas, Luis A.; Branlard, Emmanuel; Shaler, Kelsey; Vijayakumar, Ganesh; Ananthan, Shreyas; Sakievich, Philip S.; Jonkman, Jason
We study wind turbine wakes of rotors operating at high thrust coefficients (CT > 24/25) using large-eddy simulations with a rotating actuator disk model. Wind turbine wakes at high thrust coefficients are different from wakes at low thrust coefficients. Wakes behave differently at high thrust, with increased turbulence and faster recovery. Lower induction in the wake is achieved because wakes in high-thrust conditions recover much faster than in normal operating conditions. This enhanced recovery is possible thanks to the turbulence generated in the near wake. We explore the mechanism behind this behavior and propose a simple model to reproduce it. We also propose a Gaussian fit for the wakes under high-thrust conditions and use it use it to initialize an Ainslie type model within the FAST.Farm framework.
Teng, Jeffrey W.; Nergui, Delgermaa; Sepulveda-Ramos, Nelson E.; Tzintzarov, George N.; Mensah, Yaw; Cheon, Clifford D.; Rao, Sunil G.; Ringel, Brett; Gorchichko, Mariia; Li, Kan; Ying, Hanbin; Ildefonso, Adrian; Dodds, Nathaniel A.; Nowlin, Robert N.; Zhang, En X.; Fleetwood, Daniel M.; Cressler, John D.
Here, integrated silicon microwave pin diodes are exposed to 10-keV X-rays up to a dose of 2 Mrad(SiO2) and 14-MeV fast neutrons up to a fluence of 2.2×1013 cm-2. Changes in both DC leakage current and small-signal circuit components are examined. Degradation in performance due to total-ionizing dose is shown to be suppressed by non-quasi-static effects during RF operation. Tolerance to displacement damage from fast neutrons is also observed, which is explained using TCAD simulations. Overall, the characterized pin diodes are tolerant to cumulative radiation at levels consistent with space applications such as geosynchronous weather satellites.
We report a terahertz quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL) emitting around 1.9 THz with up to 10% continuous fractional frequency tuning of a single laser mode. The device shows lasing operation in pulsed mode up to 102K in a high-quality beam, with the maximum output power of 37mW and slope efficiency of 295mW/A at 77 K. Challenges for up-scaling the operating wavelength in QC metasurface VECSELs are identified.
We develop and analyze an optimization-based method for the coupling of a static peridynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.
When subjected to certain harmonic oscillations, the gas bubble in a partially liquid-filled, closed, vertical cylinder will break up. Under certain conditions, some of the gas will migrate to the bottom due to Bjerknes forces. At sufficiently large amplitudes, the bubble will break up into gas bubbles at the top and bottom ends of the cylinder. High-speed imaging captured the dynamics of bubble breakup and gas migration. Several parameters were investigated: oscillation frequency, oscillation acceleration, gas volume fraction, and liquid viscosity.
Traditional systems engineering demonstrates the importance of customer needs in scoping and defining design requirements; yet, in practice, other human stakeholders are often absent from early lifecycle phases. Human factors are often omitted in practice when evaluating and down-selecting design options due to constraints such as time, money, access to user populations, or difficulty in proving system robustness through the inclusion of human behaviors. Advances in systems engineering increasingly include non-technical influences into the design, deployment, operations, and maintenance of interacting components to achieve common performance objectives. Furthermore, such advances highlight the need to better account for the various roles of human actors to achieve desired performance outcomes in complex systems. Many of these efforts seek to infuse lessons and concepts from human factors (enhanced decision-making through Crew Resource Management), systems safety (Rasmussen's “drift toward danger”) and organization science (Giddens' recurrent human acts leading to emergent behaviors) into systems engineering to better understand how socio-technical interactions impact emergent system performance. Safety and security are examples of complex system performance outcomes that are directly impacted by varying roles of human actors. Using security performance of high consequence facilities as a representative use case, this article will outline the System Context Lenses to understand how to include various roles of human actors into systems engineering design. Several exemplar applications of this organizing lenses will be summarized and used to highlight more generalized insights for the broader systems engineering community.
Dimario, Michael; Mastin, Gary; Hodges, Ann; Lombardo, Nick; Hahn, Heidi; Professor, Nm T.
Early Systems Research and Development (ESR&D) is one of the most crucial phases in the product development process. It both blends and blurs the lines between science and engineering, and requires a risk-based, disciplined, and graded approach to effectively manage scope, cost, and complexity of the final product. Many leaders, program managers, and scientists are unwilling to involve systems engineering because of the perception that systems engineering is heavily process oriented, adds unnecessary costs, and should be applied only to mature technologies. The value of systems engineering as applied to ESR&D is unclear to these key individuals. The unfortunate result is that system engineering is not applied to ESR&D. This results in R&D efforts that may have solved the wrong problem, selected the wrong architecture, require technical rework, have difficulty transitioning later maturity levels, and result in higher R&D costs and extended development timelines. This work discusses the difficulty of introducing systems engineering to the research and early development process and their inclination perspectives of researchers, engineers, and managers. The article shall offer potential means to manage the cultural transformation of early adoption of right-sized systems engineering in ESR&D and reverse the attitudinal positions.
Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.
The U.S. Army Research Office (ARO), in partnership with IARPA, are investigating innovative, efficient, and scalable computer architectures that are capable of executing next-generation large scale data-analytic applications. These applications are increasingly sparse, unstructured, non-local, and heterogeneous. Under the Advanced Graphic Intelligence Logical computing Environment (AGILE) program, Performer teams will be asked to design computer architectures to meet the future needs of the DoD and the Intelligence Community (IC). This design effort will require flexible, scalable, and detailed simulation to assess the performance, efficiency, and validity of their designs. To support AGILE, Sandia National Labs will be providing the AGILE-enhanced Structural Simulation Toolkit (A-SST). This toolkit is a computer architecture simulation framework designed to support fast, parallel, and multi-scale simulation of novel architectures. This document describes the A-SST framework, some of its library of simulation models, and how it may be used by AGILE Performers.
LIM1TR (Lithium-Ion Modeling with 1-D Thermal Runaway) is an open-source code that uses the finite volume method to simulate heat transfer and chemical kinetics on a quasi 1-D domain. The target application of this software is to simulate thermal runaway in systems of lithium-ion batteries. The source code for LIM1TR can be found at https://github.com/ajkur/lim1tr. This user guide details the steps required to create and run simulations with LIM1TR starting with setting up the Python environment, generating an input file, and running a simulation. Additional details are provided on the output of LIM1TR as well as extending the code with custom reaction models. This user guide concludes with simple example analyses of common battery thermal runaway scenarios. The corresponding input files and processing scripts can be found in the “Examples” folder in the on-line repository, with select input files included in the appendix of this document.
In preparation for testing a lithium-helium heat exchanger at Sandia, unexpected rapid failure of the mild steel lithium preheater due to liquid metal embrittlement occurred when lithium at ~400 °C flowed into the preheater then at ~200 °C. This happened before the helium system was pressurized or heating with electron beams began. The paper presents an analysis of the preheater plus a discussion of some implications for fusion.
Understanding the capture of charge carriers by colour centres in semiconductors is important for the development of novel forms of sensing and quantum information processing, but experiments typically involve ensemble measurements, often impacted by defect proximity. Here we show that confocal fluorescence microscopy and magnetic resonance can be used to induce and probe charge transport between individual nitrogen-vacancy centres in diamond at room temperature. In our experiments, a ‘source’ nitrogen vacancy undergoes optically driven cycles of ionization and recombination to produce a stream of photogenerated carriers, one of which is subsequently captured by a ‘target’ nitrogen vacancy several micrometres away. We use a spin-to-charge conversion scheme to encode the spin state of the source colour centre into the charge state of the target, which allows us to set an upper bound to carrier injection from other background defects. We attribute our observations to the action of unscreened Coulomb potentials producing giant carrier capture cross-sections, orders of magnitude greater than those measured in ensembles.
Anwar, Ishtiaque; Hatambeigi, Mahya; Chojnicki, Kirsten; Taha, Mahmoud R.; Stormont, John C.
The stiffness of wellbore cement fracture surfaces was measured after exposing to the advective flow of nitrogen, silicone oil, and medium sweet dead crude oil for different exposure periods. The test specimens were extracted from fractured cement cylinders, where the cement fracture surfaces were exposed to the different fluids up to 15 weeks. A nanoindenter with a Berkovich indenter tip was used to measure load-indentation depth data, which was used to extract the elastic modulus (E) and nano-hardness (H) of the cement fracture surfaces. A reduction in the elastic modulus compared with an unexposed specimen were observed in all the specimens. Both elastic modulus and nano-hardness for the specimens exposed to silicone oil were lower than specimens exposed to nitrogen gas and varied with the period of exposure. The elastic modulus and nano-hardness of the specimens exposed to crude oil were the lowest with a significant decrement with the exposure period. The frequency distribution of the nanoindentation measurements shows that the volume-fraction ratio of the two types of cement hydrated nanocomposites for both the unexposed and test specimens is about 70:30%. Phase transformation beneath the indenter is observed for all of the specimens, with more obvious plastic deformation in specimens exposed to crude oil. Analytical measurements (SEM, EDS, FT-IR, and XRD) on exposed cement fracture surfaces reveal different levels of physical and chemical alteration that are consistent with the reduction in stiffness measured by nanoindentation. The study suggests that cement stiffness will decrease due to crude oil exposure, and the fracture will be sensitive to stress and pore pressure with time.
Herein, the formulation, parameter sensitivities, and usage methods for the Microstructure-Aware Plasticity (MAP) model are presented. This document is intend to serve as a reference for the underlying theory that constitutes the MAP model and as a practical guide for analysts and future developers on how aspects of this material model influence generalized mechanical behavior.
This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.
This manual describes the installation and use of the Xyce™ XDM Netlist Translator. XDM simplifies the translation of netlists generated by commercial circuit simulator tools into Xyce-compatible netlists. XDM currently supports translation from PSpice, HSPICE, and Spectre netlists into Xyce™ netlists.
This report provides a design study to produce 100% carbon-free electricity for Sandia NM and Kirtland Air Force Base (KAFB) using concentrating solar power (CSP). Annual electricity requirements for both Sandia and KAFB are presented, along with specific load centers that consume a significant and continuous amount of energy. CSP plant designs of 50 MW and 100 MW are then discussed to meet the needs of Sandia NM and the combined electrical needs of both Sandia NM and KAFB. Probabilistic modeling is performed to evaluate inherent uncertainties in performance and cost parameters on total construction costs and the levelized cost of electricity. Total overnight construction costs are expected to range between ~$300M - $400M for the 50 MW CSP plant and between ~$500M - $800M for the 100 MW plant. Annual operations and maintenance (O&M) costs are estimated together with potential offsets in electrical costs and CO2 emissions. Other considerations such as interconnections, land use and permitting, funding options, and potential agreements and partnerships with Public Service Company of New Mexico (PNM), Western Area Power Administration (WAPA), and other entities are also discussed.
We consider a class of nonlinear control synthesis problems where the underlying mathe-matical models are not explicitly known. We propose a data-driven approach to stabilize the systems when only sample trajectories of the dynamics are accessible. Our method is built on the density-function-based stability certificate that is the dual to the Lyapunov function for dynamic systems. Unlike Lyapunov-based methods, density functions lead to a convex formulation for a joint search of the control strategy and the stability certificate. This type of convex problem can be solved efficiently using the machinery of the sum of squares (SOS). For the data-driven part, we exploit the fact that the duality results in the stability theory can be understood through the lens of Perron–Frobenius and Koopman operators. This allows us to use data-driven methods to approximate these operators and combine them with the SOS techniques to establish a convex formulation of control synthesis. The efficacy of the proposed approach is demonstrated through several examples.
A common problem in developing high-reliability systems is estimating the reliability for a population of components that cannot be 100% tested. The radiation survivability of a population of components is often estimated by testing a very small sample to some multiple of the required specification level, known as an overtest. Given a successful test with a sufficient overtest margin, the population of components is assumed to have the required survivability or radiation reliability. However, no mathematical justification for such claims has been crafted without making aggressive assumptions regarding the statistics of the unknown distribution. Here we illustrate a new approach that leverages geometric bounding arguments founded on relatively modest distribution assumptions to produce conservative estimates of component reliability.
The Perovskite PV Accelerator for Commercial Technology (PACT) is an independent validation center for the evaluation of perovskite PV technologies and their bankability. The center is led by Sandia National Laboratories and the National Renewable Energy Laboratory (NREL) and includes as part of its team Los Alamos National Laboratory (LANL), CFV Labs, Black and Veatch (B&V), and the Electric Power Research Institute (EPRI). The goals of the center are to: Develop and improve indoor and outdoor performance characterization methods, Develop and validate accelerated qualification testing for early failures (5-10 years), Research degradation and failure modes, Validate outdoor performance, and Provide bankability services to US perovskite PV (PSC) industry. The importance of data and data management to the success and outcomes of the PACT center is paramount. This report describes how data will be managed and protected by PACT and identifies important data management principles that will guide our approach.
This report provides detailed documentation of the algorithms that were developed and implemented in the Plato software over the course of the Optimization-based Design for Manufacturing LDRD project.
This Content Migration Plan provides a framework and methodology for managing and executing the migration of content to the NEFC Program’s on-premises SharePoint 2016 instance, as well as guidelines regarding how to ensure that Knowledge Management Program content, both during and after the migration, is tagged properly. Analysis continues to develop a migration plan for a SharePoint Online instance in a Cloud environment.
Thermoset polymers (e.g. epoxies, vulcanizable rubbers, polyurethanes, etc.) are crosslinked materials with excellent thermal, chemical, and mechanical stability; these properties make thermoset materials attractive for use in harsh applications and environments. Unfortunately, material robustness means that these materials persist in the environment with very slow degradation over long periods of time. Balancing the benefits of material performance with sustainability is a challenge in need of novel solutions. Here, we aimed to address this challenge by incorporating boronic acid-amine complexes into epoxy thermoset chemistries, facilitating degradation of the material under pH neutral to alkaline conditions; in this scenario, water acts as an initiator to remove boron species, creating a porous structure with an enhanced surface area that makes the material more amenable to environmental degradation. Furthermore, the expulsion of the boron leaves the residual pores rich in amines which can be exploited for CO2 absorption or other functionalization. We demonstrated the formation of novel boron species from neat mixing of amine compounds with boric acid, including one complex that appears highly stable under nitrogen atmosphere up to 600 °C. While degradation of the materials under static, alkaline conditions (our “trigger”) was inconclusive at the time of this writing, dynamic conditions appeared more promising. Additionally, we showed that increasing boronic acid content created materials more resistant to thermal degradation, thus improving performance under typical high temperature use conditions.
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.
COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μr is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μcore >> μair ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.
This document summarizes the findings of a review of published literature regarding the potential impacts of electromagnetic pulse (EMP) and geomagnetic disturbance (GMD) phenomena on oil and gas pipeline systems. The impacts of telluric currents on pipelines and their associated cathodic protection systems has been well studied. The existing literature describes implications for corrosion protection system design and monitoring to mitigate these impacts. Effects of an EMP on pipelines is not a thoroughly explored subject. Most directly related articles only present theoretical models and approaches rather than specific analyses and in-field testing. Literature on SCADA components and EMP is similarly sparse and the existing articles show a variety of impacts to control system components that range from upset and damage to no effect. The limited research and the range of observed impacts for the research that has been published suggests the need for additional work on GMD and EMP and natural gas SCADA components.
In the planning for FY2020 in the U.S. DOE NE-81 Spent Fuel and Waste Science and Technology (SFWST) Campaign, the DOE requested development of a plan for activities in the Disposal Research (DR) Research and Development (R&D) over a five (5)-year period, and DOE requested periodic updates to this plan. The DR R&D 5-year plan was provided to the DOE based on the FY2020 priorities and program structure (Sassani et al., 2020) and represents a strategic guide to the work within the DR R&D technical areas (i.e., the Control Accounts), focusing on the highest priority technical thrusts. This FY2021 report is the first update to the DR R&D 5-year plan for the SFWST Campaign DR R&D activities. This 5-year plan will be a living document and is planned to be updated periodically to provide review of accomplishments and for prioritization changes based on aspects including mission progress, external technical work, and changes in SFWST Campaign objectives and/or funding levels (i.e., Program Direction). The updates to this 5-year plan will address the DR R&D that has been completed (accomplishments) and the additional knowledge gaps to be investigated, with any updates to the DR R&D priorities for the next stages of activities.
Ni-Cr alloys exhibit oscillatory segregation behaviors near low index surfaces, in which the preferred segregation species changes from Ni in the first layer to Cr in the second layer. In many dilute-alloy systems, this oscillatory pattern is attributed to the elastic release of stresses in the local lattice around the segregating solute or impurity atom. These stresses are mostly thought to originate from mismatches in the atomic size of the solute and host atoms. In Ni-Cr alloys, however, an appreciable mismatch in atomic size is not present, leading to questions about the origins of the oscillatory behavior in this alloy. Using density functional theory, we have modeled the segregation of a single Cr atom in the (100) and (111) surfaces of FCC Ni, an alloy which exhibits this oscillatory behavior. Using Bader charge analysis, we show that the negative energy correlates directly with the amount of charge on the Cr atom. As Ni atoms strip valence charge from the Cr, the Cr contracts slightly in size. The greatest contraction and highest positive charge for the Cr occurs when it is in the second layer of the surface where the system exhibits the oscillating negative segregation energy. We then find that this behavior persists in other alloy systems (Ag-Nb, Cu-Cr, Pt-Nb, and Pt-V), which exhibit similar atomic radii and electronegativity differences between host and solute to Ni-Cr. These represent alloys in which the host metal exhibits an FCC ground-state structure while the solute metal exhibits a BCC ground-state structure.
In WASH - 1400, external exposure from the finite radioactive cloud (cloudshine) is calculated by assuming that the cloud is semi-infinite, the concentration of radioactive material is uniform, and by using a correction factor to account for these approximations. This correction factor is originally based upon formulations by Healy and depends on the effective size of the plume and the distance from the plume center to the receptor. The range of the finite cloud dose correction factor table from WASH - 1400 developed using Healy formulations can be exceeded in certain situations. When the range of the table is exceeded, no extrapolation is performed; rather interpolation at the edge of the table is performed per WASH - 1400. The tabulated values of these finite cloud dose correction factors from WASH - 1400 and the interpolation at the edge of the table have been used in MACCS since its creation. An expanded table of finite cloud dose correction factors is one way to reduce the need of using interpolation at the edge of the table. The generation of an expanded finite cloud dose correction factor table for future use in MACCS is documented in this report.
A cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.
In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.
Radiographic diodes focus an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the RITS-6 Inductive Voltage Adder (IVA) driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode-cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analyzation of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below ~ few mm. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle-behavior of SMP shots is made to lend credence to the inference. There is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.
This project demonstrates that Chapel programs can interface with MPI-based libraries written in C++ without storing multiple copies of shared data. Chapel is a language for productive parallel computing using global address spaces (PGAS). We identified two approaches to interface Chapel code with the MPI-based Grafiki and Trilinos libraries. The first uses a single Chapel executable to call a C function that interacts with the C++ libraries. The second uses the mmap function to allow separate executables to read and write to the same block of memory on a node. We also encapsulated the second approach in Docker/Singularity containers to maximize ease of use. Comparisons of the two approaches using shared and distributed memory installations of Chapel show that both approaches provide similar scalability and performance.
Background: Blockchain distributed ledger technology is just starting to be adopted in genomics and healthcare applications. Despite its increased prevalence in biomedical research applications, skepticism regarding the practicality of blockchain technology for real-world problems is still strong and there are few implementations beyond proof-of-concept. We focus on benchmarking blockchain strategies applied to distributed methods for sharing records of gene-drug interactions. We expect this type of sharing will expedite personalized medicine. Basic Procedures: We generated gene-drug interaction test datasets using the Clinical Pharmacogenetics Implementation Consortium (CPIC) resource. We developed three blockchain-based methods to share patient records on gene-drug interactions: Query Index, Index Everything, and Dual-Scenario Indexing. Main Findings: We achieved a runtime of about 60 s for importing 4,000 gene-drug interaction records from four sites, and about 0.5 s for a data retrieval query. Our results demonstrated that it is feasible to leverage blockchain as a new platform to share data among institutions. Principal Conclusions: We show the benchmarking results of novel blockchain-based methods for institutions to share patient outcomes related to gene-drug interactions. Our findings support blockchain utilization in healthcare, genomic and biomedical applications. The source code is publicly available at https://github.com/tsungtingkuo/genedrug.
A mass property calculator has been developed to compute the moment of inertia properties of an assemblage of parts that make up a system. The calculator can take input from spreadsheets or Creo mass property files or it can be interfaced with Phoenix Integration Model Center. The input must include the centroidal moments of inertia of each part with respect to its local coordinates, the location of the centroid of each part in the system coordinates and the Euler angles needed to rotate from the part coordinates to the system coordinates. The output includes the system total mass, centroid and mass moment of inertia properties. The input/output capabilities allow the calculator to interface with external optimizers. In addition to describing the calculator, this document serves as its user's manual. The up-to-date version of the calculator can be found in the Git repository https://cee-gitlab.sandia.gov/cj?ete/mass-properties-calculator.
Drilling systems that use downhole rotation must react torque either through the drill-string or near the motor to achieve effective drilling performance. Problems with drill-string loading such as buckling, friction, and twist become more severe as hole diameter decreases. Therefore, for small holes, reacting torque downhole without interfering with the application of weight-on-bit, is preferred. In this paper, we present a novel mechanism that enables effective and controllable downhole weight on bit transmission and torque reaction. This scalable design achieves its unique performance through four key features: (1) mechanical advantage based on geometry, (2) direction dependent behavior using rolling and sliding contact, (3) modular scalability by combining modules in series, and (4) torque reaction and weight on bit that are proportional to applied axial force. As a result, simple mechanical devices can be used to react large torques while allowing controlled force to be transmitted to the drill bit. We outline our design, provide theoretical predictions of performance, and validate the results using full-scale testing. The experimental results include laboratory studies as well as limited field testing using a percussive hammer. These results demonstrate effective torque reaction, axial force transmission, favorable scaling with multiple modules, and predictable performance that is proportional to applied force.
Reeves, Michael J.; Tian, Dave J.; Bianchi, Antonio; Berkay Celik ZBerkay C.
Container escapes enable the adversary to execute code on the host from inside an isolated container. Notably, these high severity escape vulnerabilities originate from three sources: (1) container profile misconfigurations, (2) Linux kernel bugs, and (3) container runtime vulnerabilities. While the first two cases have been studied in the literature, no works have investigated the impact of container runtime vulnerabilities. In this paper, to fill this gap, we study 59 CVEs for 11 different container runtimes. As a result of our study, we found that five of the 11 runtimes had nine publicly available PoC container escape exploits covering 13 CVEs. Our further analysis revealed all nine exploits are the result of a host component leaked into the container. Here, we apply a user namespace container defense to prevent the adversary from leveraging leaked host components and demonstrate that the defense stops seven of the nine container escape exploits.
Echeverria, Marco J.; Galitskiy, Sergey; Mishra, Avanish; Dingreville, Remi P.; Dongare, Avinash M.
A hybrid atomic-scale and continuum-modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser-loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. The simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.
Energy utilities are evaluating emerging energy technologies to reduce reliance on carbon as an energy carrier. Hydrogen has been identified as a potential substitute for carbon-based fuels that can be blended into other gaseous energy carriers, such as natural gas. However, hydrogen blending into natural gas has important implications on safety which need to be evaluated. Designers and installers of systems that utilize hydrogen gas blending into natural gas distribution systems need to adhere to local building codes and engage with the authority having jurisdiction (AHJ) for safety and permitting approvals. These codes and standards must be considered to understand where safety gaps might be apparent when injecting hydrogen into the natural gas infrastructure. This report generates a list of relevant codes and standards for hydrogen blending on existing, upgraded, or new pipelines. Additionally, a preliminary assessment was made to identify the codes and standards that need to be modified to enable this technology as well as potential gaps due to the unique nature and safety concerns of gaseous hydrogen.
In this work, we have used the well-understood quantum Hall (QH) stripes in high quality two-dimensional electron gases to mimic charge stripes in high transition temperature (Tc) superconductors. The science question we want to address is “Can QH stripes mimic high Tc superconductor stripes and provide a controlled experimental setup to pin-down the role of stripes in high Tc superconductivity?”. We have observed anomalous superconducting transition like behavior in GaAs double quantum well systems (DQWs) when each quantum well (QW) is tuned to the charge stripe states but with different Landau level fillings. Furthermore, we have shown that the transition like behavior is sharper in the DQWs when the two QWs are more strongly coupled. Our results suggest, for the first time, experimental evidence of the paired charge stripes model, which might lead to room-temperature superconductors that have enormously wide applications in computing, energy, and transportation industries. Advancing the science of high transition temperature superconductivity will have a profound impact in advancing energy technologies, ranging from the next generation microchips, new energy transfer grid to public transportation, and thus is important to nation’s energy security and relevant across the landscape of many mission spaces. Sandia has been a leader in materials science research and development. The proposed research takes advantage of Sandia’s state-ofthe-art MBE facilities at the Center for Integrated Nanotechnologies (CINT) and utilizes Sandia’s extensive advanced materials characterization resources. We envision a significant impact on the nation’s energy research and security challenges by investing in this research.
Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. A survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.
Current wind turbine blade materials may not be damage tolerant to the extent necessary to optimize the Levelized Cost of Energy (LCOE) of wind energy plants. Traditionally, wind turbine blades have been designed using a safe-life approach, but advances in inspection techniques and structural health monitoring solutions give rise to the opportunity to design wind turbine blades using a damage tolerant approach. Materials selection is a key element of da mage tolerant design, so the extent of the damage tolerance of alternative materials has been analyzed through a literature review and discussions with industry leaders. Fabrics and resin selection significantly affect the damage tolerance of composites. Changes to fabric architecture may include through-the-thickness (TTT) fibers, stretch-broken carbon fiber (SBCF) composites, and aligned discontinuous fiber reinforced composites (ADFRCs). Previous research has demonstrated that using TTT fibers in creases damage tolerance, but additional research is necessary to demonstrate the effectiveness of SBCFs and ADFRCs in mitigating damage. Several studies have demonstrated increased damage tolerance when toughened resin systems are used. In addition to toughened resin systems, thermoplastics have been shown to be tougher than thermosets. However, thermosets have been traditionally preferred in wind turbine blade manufacturing due to ease of manufacturing. Thermoplastic resin system s have been developed that can be used with conventional manufacturing methods but have yet to be studied for its damage tolerant capabilities. Furthermore, cost and stress analyses on where to effectively implement TTT fibers, SBCF composites, ADFRCs, and toughened resin systems must be executed prior to incorporating new materials into wind turbine blade manufacturing.
Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke's law models) are implemented to describe mechanical characteristics of the system. Our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide.
We successfully demonstrated the utility of surface science techniques - namely scanning probe microscopy and thermal desorption spectroscopy - on three different material systems: incipient soot formed during fossil fuel combustion, surface oxides passivating polycrystalline nickel hydrogen uptake, and aluminum hydride cluster formation underpinning solid-state hydrogen fuel storage. For all three material systems, surface science techniques haven proven to probe intricate nanoscale phenomena that are critical to macroscale material behavior. This LDRD has gained insight into early-stage pollution formation, the impacts of common contaminants on tritium flow regulation, and the limitations of solid-state hydrogen fuel storage. Our results support the diversification of national energy technologies.
This report provides detailed documentation of the algorithms that where developed and implemented in the Plato software over the course of the Optimization-based Design for Manufacturing LDRD project.
Co-deposited, immiscible alloy systems form hierarchical microstructures under specific deposition conditions that accentuate the difference in constituent element mobility. The mechanism leading to the formation of these unique hierarchical morphologies during the deposition process is difficult to identify, since the characterization of these microstructures is typically carried out post-deposition. We employ phase-field modeling to study the evolution of microstructures during deposition combined with microscopy characterization of experimentally deposited thin films to reveal the origin of the formation mechanism of hierarchical morphologies in co-deposited, immiscible alloy thin films. Our results trace this back to the significant influence of a local compositional driving force that occurs near the surface of the growing thin film. We show that local variations in the concentration of the vapor phase near the surface, resulting in nuclei (i.e., a cluster of atoms) on the film’s surface with an inhomogeneous composition, can trigger the simultaneous evolution of multiple concentration modulations across multiple length scales, leading to hierarchical morphologies. We show that locally, the concentration must be above a certain threshold value in order to generate distinct hierarchical morphologies in a single domain.