Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.8 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

Abstract not provided.

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

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.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In total, light-duty vehicles in the United States travel on the order of 3 trillion miles annually, providing tremendous societal and personal benefits. However, the environmental burden is excessive, prompting Co-Optimization of Fuel and Engines (Co-Optima) program efforts to provide the science needed to increase engine efficiency and produce non-fossil fuels with reduced greenhouse gas emissions. Boosted spark-ignition (SI) engines provide high power density by offering high loads and engine speeds, making them light-weight and attractive for light-duty vehicles. Unfortunately, the engine efficiency drops off at lower loads and speeds, where the engine spends most time during typical driving. Multimode SI engines can use a more efficient advanced lean combustion mode at lower loads and speeds, while reverting to boosted SI under high-load conditions. Within Co-Optima, multiple advanced lean combustion modes have been explored; these include stratified-charge SI, pre-chamber lean SI, and advanced compression ignition (ACI) techniques such as spark-assisted compression ignition (SACI). For these combustion modes, focus has been on determining fuel properties that enable higher engine efficiency, clean and stable combustion, and effective exhaust aftertreatment. This report highlights recent efforts funded by the Vehicle Technologies Office at multiple National Laboratories that supported the multimode project in Co-Optima. It also includes a brief summary of biofuel production research funded by the Bioenergy Technologies Office.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Aims: High-resolution information on soils’ vulnerability to climate-induced soil organic carbon (SOC) loss can enable environmental scientists, land managers, and policy makers to develop targeted mitigation strategies. This study aims to estimate baseline and decadal changes in continental US surface SOC stocks under future emission scenarios. Location: Continental United States. Time period: 2014–2100. Methods: We used recent SOC field observations (n = 6,213 sites), environmental factors (n = 32), and an ensemble machine learning (ML) approach to estimate baseline SOC stocks in surface soils across the continental United States at 100-m spatial resolution, and decadal changes under the projected climate scenarios of Coupled Model Intercomparison Project Phase Six (CMIP6) earth system models (ESMs). Results: Baseline SOC projections from ML approaches captured more than 50% of variability in SOC observations, whereas ESMs represented only 6–16% of observed SOC variability. ML estimates showed a mean total loss of 1.8 Pg C from US surface soils under the high-emission scenario by 2100, whereas ESMs showed no significant change in SOC stocks with wide variation among ESMs. Both ML and ESM predictions agree on the direction of SOC change (net emissions or sequestration) across 46–51% of continental US land area. These differences are attributable to the high-resolution site-specific data used in the ML models compared to the relatively coarse grid represented in CMIP6 ESMs. Main conclusions: Our high-resolution estimates of baseline SOC stocks, identification of key environmental controllers, and projection of SOC changes from US land cover types under future climate scenarios suggest the need for high-resolution simulations of SOC in ESMs to represent the heterogeneity of SOC. We found that the SOC change is sensitive to key soil related factors (e.g. soil drainage and soil order) that have not been historically considered as input parameters in ESMs, because currently more than 95% variability in the SOC of CMIP6 ESMs is controlled by net primary productivity, temperature, and precipitation. Using additional environmental factors to estimate the baseline SOC stocks and predict the future trajectory of SOC change can provide more accurate results.

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.

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.

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.

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.

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.

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.

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.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This manual gives usage information for the Charon semiconductor device simulator. Charon was developed to meet the modeling needs of Sandia National Laboratories and to improve on the capabilities of the commercial TCAD simulators; in particular, the additional capabilities are running very large simulations on parallel computers and modeling displacement damage and other radiation effects in significant detail. The parallel capabilities are based around the MPI interface which allows the code to be ported to a large number of parallel systems, including linux clusters and proprietary “big iron” systems found at the national laboratories and in large industrial settings.

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.

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.

Worldwide growth in electric vehicle use is prompting new installations of private and public electric vehicle supply equipment (EVSE). EVSE devices support the electrification of the transportation industry but also represent a linchpin for power systems and transportation infras-tructures. Cybersecurity researchers have recently identified several vulnerabilities that exist in EVSE devices, communications to electric vehicles (EVs), and upstream services, such as EVSE vendor cloud services, third party systems, and grid operators. The potential impact of attacks on these systems stretches from localized, relatively minor effects to long-term national disruptions. Fortunately, there is a strong and expanding collection of information technology (IT) and operational technology (OT) cybersecurity best practices that may be applied to the EVSE environment to secure this equipment. In this paper, we survey publicly disclosed EVSE vulnerabilities, the impact of EV charger cyberattacks, and proposed security protections for EV charging technologies.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Presented in this document are the theoretical aspects of capabilities contained in the Sierra/SM code. This manuscript serves as an ideal starting point for understanding the theoretical foundations of the code. For a comprehensive study of these capabilities, the reader is encouraged to explore the many references to scientific articles and textbooks contained in this manual. It is important to point out that some capabilities are still in development and may not be presented in this document. Further updates to this manuscript will be made as these capabilities come closer to production level.

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.

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.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.