Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present a combination of machine-learned models that predicts the surface elastic properties of general free surfaces in face-centered cubic (FCC) metals. These models are built by combining a semi-analytical method based on atomistic simulations to calculate surface properties with the artificial neural network (ANN) method or the boosted regression tree (BRT) method. The latter is also used to link bulk properties and surface orientation to surface properties. The surface elastic properties are represented by their invariants considering plane elasticity within a polar method. The resulting models are shown to accurately predict the surface elastic properties of seven pure FCC metals (Cu, Ni, Ag, Au, Al, Pd, Pt). The BRT model reveals the correlations between bulk and corresponding surface properties in terms of invariants, which can be used to guide the design of complex nano-sized particles, wires and films. Finally, by expressing the surface excess energy density as a function of surface elastic invariants, fast predictions of surface energy as a function of in-plane deformations can be made from these model constructs.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

For decades, it has been observed that the commonly used Borgnakke-Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke-Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke-Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke-Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). The consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The optically segmented single volume scatter camera (OS-SVSC) aims to image neutron sources for nuclear nonproliferation applications using the kinematic reconstruction of elastic double-scatter events. We report on the design, construction, and calibration of one module of a new prototype. The module includes 16 EJ-204 organic plastic scintillating bars individually wrapped in Teflon tape, each measuring 0.5 {\mathrm {cm}}\times 0.5 {\mathrm {cm}}\times 20 {\mathrm {cm}}. The scintillator array is coupled to two custom silicon photomultiplier (SiPM) boards consisting of a 2\times 8 array of SensL J-Series-60035 SiPMs, which are read out by a custom 16 channel DRS4 based digitizer board. The electrical crosstalk between SiPMs within the electronics chain is measured as 0.76\% \,\pm \,0.11\% among all 16 channels. We report the detector response of one module including interaction position, time, and energy, using two different optical coupling materials: EJ-560 silicone rubber optical coupling pads and EJ-550 optical coupling grease. We present results in terms of the overall mean and standard deviation of the z -position reconstruction and interaction time resolutions for all 16 bars in the module. We observed the 1\sigma z -position resolution for gamma interactions in the 0.3-0.4 MeVee range to be 2.24 cm ± 1.10 cm and 1.45 cm ± 0.19 cm for silicone optical coupling pad and optical grease, respectively. The observed 1\sigma interaction time resolution is 265 ps ± 29 ps and 235 ps ± 10 ps for silicone optical coupling pad and optical grease, respectively.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Errors in quantum logic gates are usually modeled by quantum process matrices (CPTP maps). But process matrices can be opaque and unwieldy. We show how to transform the process matrix of a gate into an error generator that represents the same information more usefully. We construct a basis of simple and physically intuitive elementary error generators, classify them, and show how to represent the error generator of any gate as a mixture of elementary error generators with various rates. Finally, we show how to build a large variety of reduced models for gate errors by combining elementary error generators and/or entire subsectors of generator space. We conclude with a few examples of reduced models, including one with just 9N2 parameters that describes almost all commonly predicted errors on an N-qubit processor.

Abstract not provided.

Electrochemical characteristics and semiconducting behavior of additively manufactured electron beam melted (EBM) and wrought (WR) Ti–6Al–4V (Ti-G5) are compared in Ringer’s physiological solution. X-ray diffraction (XRD) and field emission scanning electron microscopy (FE-SEM) confirmed the α + β structure of the tested materials, with two different microstructure types of “bimodal” and “basket-weave” for WR and EBM, respectively. Potentiodynamic polarization (PDP) revealed that the corrosion current density for EBM (icorr = 0.27 ± 0.06 μA cm−2) is less than the WR (icorr = 0.70 ± 0.05 μA cm−2). Moreover, potentiostatic polarization (PS) that was employed to form the passive layers at three different potentials of 300, 500, and 700 mVAg/AgCl, showed that the passive films on the EBM sample are thinner. This finding was confirmed by electrochemical impedance spectroscopy (EIS). Furthermore, through Mott–Schottky (M–S) analysis, donor densities on WR passive films were found to be ~ 1.5 times larger than EBM. Although PS and EIS confirmed that the passive layer on EBM is thinner, it provides higher corrosion resistance than WR. The passive layer on both samples were found to have n-type characteristics with a duplex structure. Graphical abstract: [Figure not available: see fulltext.]

Abstract not provided.

Abstract not provided.

Abstract not provided.

Advances in machine learning algorithms and increased computational efficiencies give engineers new capabilities and tools to apply to engineering design. Machine learning models can approximate complex functions and, therefore, can be useful for various tasks in the engineering design workflow. This paper investigates using reinforcement learning (RL), a subset of machine learning that teaches an agent to complete a task through accumulating experiences in an interactive environment, to automate the designing of 2D discretized topologies. RL agents use past experiences to learn sequential sets of actions to best achieve some objective. In the proposed environment, an RL agent can make sequential decisions to design a topology by removing elements to best satisfy compliance minimization objectives. After each action, the agent receives feedback by evaluating how well the current topology satisfies the design objectives. After training, the agent was tasked with designing optimal topologies under various load cases. The agent's proposed designs had similar or better compliance minimization performance to those produced by traditional gradient-based topology optimization methods. These results show that a deep RL agent can learn generalized design strategies to satisfy multi-objective design tasks and, therefore, shows promise as a tool for arbitrarily complex design problems across many domains.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. The novel SM proposed here is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (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. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Metals subjected to irradiation environments undergo microstructural evolution and concomitant degradation, yet the nanoscale mechanisms for such evolution remain elusive. Here, we combine in situ heavy ion irradiation, atomic resolution microscopy, and atomistic simulation to elucidate how radiation damage and interfacial defects interplay to control grain boundary (GB) motion. While classical notions of boundary evolution under irradiation rest on simple ideas of curvature-driven motion, the reality is far more complex. Focusing on an ion-irradiated Pt Σ3 GB, we show how this boundary evolves by the motion of 120° facet junctions separating nanoscale {112} facets. Our analysis considers the short- and mid-range ion interactions, which roughen the facets and induce local motion, and longer-range interactions associated with interfacial disconnections, which accommodate the intergranular misorientation. We suggest how climb of these disconnections could drive coordinated facet junction motion. These findings emphasize that both local and longer-range, collective interactions are important to understanding irradiation-induced interfacial evolution.

Abstract not provided.

This report presents a specification for the Portals 4 network programming interface. Portals 4 is intended to allow scalable, high-performance network communication between nodes of a parallel computing system. Portals 4 is well suited to massively parallel processing and embedded systems. Portals 4 represents an adaption of the data movement layer developed for massively parallel processing platforms, such as the 4500-node Intel TeraFLOPS machine. Sandia's Cplant cluster project motivated the development of Version 3.0, which was later extended to Version 3.3 as part of the Cray Red Storm machine and XT line. Version 4 is targeted to the next generation of machines employing advanced network interface architectures that support enhanced offload capabilities.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Although mass spectrometry is a widely used analytical tool, age-appropriate, interactive outreach activities for laboratory visitors, especially children, are lacking. The presented interactive demonstration, "Did you eat a MOLEcule today?", introduces all ages to molecular weight concepts and mass spectrometry in a research laboratory, while connecting the concepts to real-world applications. Through real-time breath analysis, participants explore the concepts of molecular weight, electrostatic field manipulation of charged molecules, and analyte identification by mass analysis. This module is rapid and highly adaptable for outreach activities but also includes age- or classroom-appropriate variations to decrease or increase difficulty levels. The presented interactive demonstration has repeatedly been implemented, with over 2300 participants during six annual "Take Our Daughters & Sons to Work Day" and two corporate "Family Day" outreach activities, successfully engaging, exciting, and educating both kids and parents.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

High gain in hotspot-ignition inertial confinement fusion (ICF) implosions requires the propagation of thermonuclear burn from a central hotspot to the surrounding cold dense fuel. As ICF experiments enter the burning plasma regime, diagnostic signatures of burn propagation must be identified. In previous work [A. J. Crilly et al., Phys. Plasmas 27(1), 012701 (2020)], it has been shown that the spectral shape of the neutron backscatter edges is sensitive to the dense fuel hydrodynamic conditions. The backscatter edges are prominent features in the ICF neutron spectrum produced by the 180° scattering of primary deuterium–tritium fusion neutrons from ions. In this work, synthetic neutron spectra from radiation-hydrodynamics simulations of burning ICF implosions are used to assess the backscatter edge analysis in a propagating burn regime. Significant changes to the edge's spectral shape are observed as the degree of burn increases, and a simplified analysis is developed to infer scatter-averaged fluid velocity and temperature. The backscatter analysis offers direct measurement of the increased dense fuel temperatures that result from burn propagation.

Abstract not provided.

We present a computational study and framework that allows us to study and understand the crack nucleation process from forging flaws. Forging flaws may be present in large steel rotor components commonly used for rotating power generation equipment including gas turbines, electrical generators, and steam turbines. The service life of these components is often limited by crack nucleation and subsequent growth from such forging flaws, which frequently exhibit themselves as non-metallic oxide inclusions. The fatigue crack growth process can be described by established engineering fracture mechanics methods. However, the initial crack nucleation process from a forging flaw is challenging for traditional engineering methods to quantify as it depends on the details of the flaw, including flaw morphology. We adopt the peridynamics method to describe and study this crack nucleation process. For a specific industrial gas turbine rotor steel, we present how we integrate and fit commonly known base material property data such as elastic properties, yield strength, and S-N curves, as well as fatigue crack growth data into a peridynamic model. The obtained model is then utilized in a series of high-performance two-dimensional peridynamic simulations to study the crack nucleation process from forging flaws for ambient and elevated temperatures in a rectangular simulation cell specimen. The simulations reveal an initial local nucleation at multiple small oxide inclusions followed by micro-crack propagation, arrest, coalescence, and eventual emergence of a dominant micro-crack that governs the crack nucleation process. The dependence on temperature and density of oxide inclusions of both the details of the microscopic processes and cycles to crack nucleation is also observed. The results are compared with fatigue experiments performed with specimens containing forging flaws of the same rotor steel.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Glacier forelands provide an excellent opportunity to investigate vegetation succession and soil development along the chronosequence; however, there are few studies on soil biogeochemical changes from environmental factors, aside from time. This study aimed to investigate soil development and biogeochemical changes in the glacier foreland of Midtre Lovénbreen, Svalbard, by considering various factors, including time. Eighteen vegetation and soil variables were measured at 38 different sampling sites of varying soil age, depth, and glacio-fluvial activity. Soil organic matter (SOM) was quantitatively measured, and the compositional changes in SOM were determined following size-density fractionation. In the topsoil, the soil organic carbon (SOC) and total nitrogen (N) content was found to increase along the soil chronosequence and were highly correlated with vegetation-associated variables. These findings suggest that plant-derived material was the main driver of the light fraction of SOM accumulation in the topsoil. The heavy fractions of SOM were composed of microbially transformed organic compounds, eventually contributing to SOM stabilization within short 90-yr deglaciation under harsh climatic conditions. In addition to time, the soil vertical profiles showed that other environmental parameters, also affected the soil biogeochemical properties. The high total phosphorous (P) content and electrical conductivity in the topsoil were attributed to unweathered subglacial materials and a considerable amount of inorganic ions from subglacial meltwater. The high P and magnesium content in the subsoil were attributed to parent materials, while the high sodium and potassium content in the surface soil were a result of sea-salt deposition. Glacio-fluvial runoff hampered ecosystem development by inhibiting vegetation development and SOM accumulation. This study emphasizes the importance of considering various soil-forming factors, including parent/subglacial materials, aeolian deposition, and glacio-fluvial runoff, as well as soil age, to obtain a comprehensive understanding of the ecosystem development in glacier forelands.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, nonsooting flames.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Response to elongational flow is fundamental to soft matter and directly impacts new developments in a broad range of technologies form polymer processing and microfluidics to controlled flow in biosystems. Of particular significance are the effects of elongational flow on self-assembled systems where the interactions between the fundamental building blocks control their adaptation. Here we probe the effects of associating groups on the structure and dynamics of linear polymer melts in uniaxial elongation using molecular dynamics simulations. We study model polymers with randomly incorporated backbone associations with interaction strengths varying from 1kBT to 10kBT. These associating groups drive the formation of clusters in equilibrium with an average size that increases with interaction strength. Flow drives these clusters to continuously break and reform as chains stretch. These flow-driven cluster dynamics drive a qualitative transition in polymer elongation dynamics from homogeneous to nanoscale localized yield and cavitation as the association strength increases.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This article describes a calculation of the spontaneous emission limited linewidth of a semiconductor laser consisting of hybrid or heterogeneously integrated, silicon and III–V intracavity components. Central to the approach are a) description of the multi-element laser cavity in terms of composite laser/free-space eigenmodes, b) use of multimode laser theory to treat mode competition and multiwave mixing, and c) incorporation of quantum-optical contributions to account for spontaneous emission effects. Application of the model is illustrated for the case of linewidth narrowing in an InAs quantum-dot laser coupled to a high- (Formula presented.) SiN cavity.

Abstract not provided.

Abstract not provided.

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and in this work, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

Abstract not provided.

Abstract not provided.

Abstract not provided.