Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, in this paper, we present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni (111) symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.
Khan, Rabia; Miller, Michelle; Valerio, Analise; Willhite, Isabella; Hernandez, Zachary
A need for a standardized Division 10000 onboarding program for virtual hires was identified by management to formalize the way employees and interns are onboarded and trained into Division 10000. This white paper provides effective short and long-term suggestions in the efforts of improving virtual onboarding. Data suggests that remote work is going to become the forefront of many industry practices, which indicates the need of a standardized virtual onboarding practices. With our research, gap assessments, benchmarking, and conducting interviews both internally and externally, we found that clarity, culture, and connection proved to be the strongest solutions in order to maintain Sandia’s competitive edge and sustain workers both remote and in-person.
Monti, Juan M.; Misra, Deepankar; Weck, Philippe F.; Rivarola, Roberto D.; Tribedi, Lokesh C.
Absolute double differential cross sections (DDCS) of electrons emitted from uracil and 5-bromouracil (BrU) in collisions with protons of energy 200 keV have been measured for various forward and backward emission angles over wide range of electron energies. The measured DDCS are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) calculations. The optimized structure of the BrU was estimated along with the population analysis of all the occupied orbitals using a self-consistent field density. A comparison between the measured DDCS data for the two molecules show that the cross section of low energy electrons emitted from BrU is substantially larger than that for uracil. The BrU-to-uracil DDCS ratios obtained from the present measurements indicate an enhancement of the electron emission by a factor which is as large as 2.0 to 2.5. These electrons being the major agent for damaging the DNA/RNA of the malignant tissues, the present results are expected to provide an important input for the radiosensitization effect in hadron therapy. It is noteworthy to mention that the CDW-EIS calculations for Coulomb ionization cannot predict such enhancement. A large angular asymmetry is observed for uracil with a broad structure, which is absent in case of BrU.
An analytical expression for turbulent kinematic viscosity (vt), based solely on the hydraulic Reynolds number (Re), was derived and evaluated. The analytical expression is valid for the fast estimation of vt for internal, isotropic, fully developed flows. The expression was compared with experimental and simulation data for air, water, and liquid sodium, and was shown to provide reasonable values for 2100 ≤ Re ≤ 3.6 × 106 and Prandtl number (Pr) range of 0.0107 ≤ Pr ≤ 5.65. In addition, new expressions suitable for the central portion of internal flows, away from the wall, were derived for the turbulent Reynolds number (ReT), showing its relationship to Re, as well as to the ratio of vt and the molecular kinematic viscosity (v).
Ducted fuel injection (DFI) is a novel combustion strategy that has been shown to significantly attenuate soot formation in diesel engines. While previous studies have used optical diagnostics and optical filter smoke number methods to show that DFI reduces in-cylinder soot formation and engine-out soot emissions, respectively, this is the first study to measure solid particle number (PN) emissions in addition to particle mass (PM). Furthermore, this study quantitatively evaluates the use of transient particle instruments for measuring particles from skip-fired operation in an optical single cylinder research engine (SCRE). Engine-out PN was measured using an engine exhaust particle sizer following a catalytic stripper, and PM was measured using a photoacoustic analyzer. The study improves on earlier preliminary emissions studies by clearly showing that DFI reduces overall PM by 76%–79% and PN for particles larger than 23 nm by 77% relative to conventional diesel combustion at a 1200-rpm, 13.3-bar gross indicated mean effective pressure operating condition. The degree of engine-out PM reduction with DFI was similar across both particulate measurement instruments used in the work. Through the use of bimodal distribution fitting, DFI was also shown to reduce the geometric mean diameter of accumulation mode particles by 26%, similar to the effects of increased injection pressure in conventional diesel combustion systems. This work clearly shows the significant solid particulate matter reductions enabled by DFI while also demonstrating that engine-out PN can be accurately measured from an optical SCRE operating in a skip-fired mode. Based on these results, it is believed that DFI has the potential to enable fuel savings when implemented in multi-cylinder engines, both by lowering the required frequency of active diesel particulate filter regeneration, and by reducing the backpressure imposed by exhaust filtration systems.
We present an adaptive algorithm for constructing surrogate models of multi-disciplinary systems composed of a set of coupled components. With this goal we introduce “coupling” variables with a priori unknown distributions that allow surrogates of each component to be built independently. Once built, the surrogates of the components are combined to form an integrated-surrogate that can be used to predict system-level quantities of interest at a fraction of the cost of the original model. The error in the integrated-surrogate is greedily minimized using an experimental design procedure that allocates the amount of training data, used to construct each component-surrogate, based on the contribution of those surrogates to the error of the integrated-surrogate. The multi-fidelity procedure presented is a generalization of multi-index stochastic collocation that can leverage ensembles of models of varying cost and accuracy, for one or more components, to reduce the computational cost of constructing the integrated-surrogate. Extensive numerical results demonstrate that, for a fixed computational budget, our algorithm is able to produce surrogates that are orders of magnitude more accurate than methods that treat the integrated system as a black-box.
Achieving practical, high-energy-density calcium batteries requires controlling the stability of Ca2+electrolytes during calcium metal cycling. Because of the highly reactive nature of calcium, most typical electrolyte constituents are unstable, leading to electrode passivation and low Coulombic efficiency. Among various commercially available salts, calcium bis(trifluoromethylsulfonyl)imide (Ca(TFSI)2) is attractive because of its oxidative stability and high solubility in a variety of solvents. However, this salt does not allow for calcium metal plating, and it has been proposed that TFSI-instability induced by Ca2+coordination is to blame. In this work, we test the ability of strongly coordinating Ca2+cosalts such as halides and borohydrides to displace TFSI-from the first coordination shell of Ca2+and thereby stabilize TFSI-based electrolytes to enable calcium plating. Through spectroscopic analysis, we find that the effectiveness of these cosalts at displacing the TFSI-anion is dependent on the solvent's coordination strength toward Ca2+. Surprisingly, electrochemical calcium deposition behavior is not correlated to the population of bound or free TFSI-. Instead, the nature of the coordination interaction between Ca2+and the cosalt anion is more important for determining stability. Our findings indicate that TFSI-anions are inherently unstable during calcium deposition even in the nominally free state. Therefore, strategies aimed at eliminating the interactions of these anions with the electrode surface via interface/interphase design are required.
A series of experiments will be performed to test the integral effects of molybdenum on the reactivity of a critical system. These experiments will use the 7uPCX assembly with the 1.55 cm triangular pitch grid plates. Molybdenum sleeves, consisting of 19.6 inch long 0.5-inch nominal outside diameter molybdenum tubes with 0.031-inch nominal wall thickness and centering hardware, will be placed on some of the fuel rods in the array. The purpose of this analysis is to examine two configurations of the 7uPCX using the 1.55 cm triangular pitch grid plates in fully-reflected approach-to-critical experiments with the number of fuel rods in the array as the approach parameter. This document presents the results of the analysis that was done to allow completion of the 7uPCX Configuration Checklist from Appendix A of SPRF-AP-005 [SNL 2020] for the cores noted above. The checklists for these cores are shown in Appendix A.
The inverse methods team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of weapons systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.
Thornton, Chappel S.; Tuttle, Blair; Turner, Emily; Law, Mark E.; Pantelides, Sokrates T.; Wang, George T.; Jones, Kevin S.
A recently discovered, enhanced Ge diffusion mechanism along the oxidizing interface of Si/SiGe nanostructures has enabled the formation of single-crystal Si nanowires and quantum dots embedded in a defect-free, single-crystal SiGe matrix. Here, we report oxidation studies of Si/SiGe nanofins aimed at gaining a better understanding of this novel diffusion mechanism. A superlattice of alternating Si/Si0.7Ge0.3layers was grown and patterned into fins. After oxidation of the fins, the rate of Ge diffusion down the Si/SiO2interface was measured through the analysis of HAADF-STEM images. The activation energy for the diffusion of Ge down the sidewall was found to be 1.1 eV, which is less than one-quarter of the activation energy previously reported for Ge diffusion in bulk Si. Through a combination of experiments and DFT calculations, we propose that the redistribution of Ge occurs by diffusion along the Si/SiO2interface followed by a reintroduction into substitutional positions in the crystalline Si.
Physics-informed machine learning (PIML) has emerged as a promising new approach for simulating complex physical and biological systems that are governed by complex multiscale processes for which some data are also available. In some instances, the objective is to discover part of the hidden physics from the available data, and PIML has been shown to be particularly effective for such problems for which conventional methods may fail. Unlike commercial machine learning where training of deep neural networks requires big data, in PIML big data are not available. Instead, we can train such networks from additional information obtained by employing the physical laws and evaluating them at random points in the space-time domain. Such PIML integrates multimodality and multifidelity data with mathematical models, and implements them using neural networks or graph networks. Here, we review some of the prevailing trends in embedding physics into machine learning, using physics-informed neural networks (PINNs) based primarily on feed-forward neural networks and automatic differentiation. For more complex systems or systems of systems and unstructured data, graph neural networks (GNNs) present some distinct advantages, and here we review how physics-informed learning can be accomplished with GNNs based on graph exterior calculus to construct differential operators; we refer to these architectures as physics-informed graph networks (PIGNs). We present representative examples for both forward and inverse problems and discuss what advances are needed to scale up PINNs, PIGNs and more broadly GNNs for large-scale engineering problems.
Abbott, Benjamin W.; Brown, Michael; Carey, Joanna C.; Ernakovich, Jessica; Frederick, Jennifer M.; Guo, Laodong; Lee, Raymond M.; Loranty, Michael M.; Macdonald, Robie; Mann, Paul J.; Natali, Susan M.; Olefeldt, David; Pearson, Pam; Rec, Abigail; Robards, Martin; Salmon, Verity G.; Sayedi, Sayedeh S.; Schadel, Christina; Schuur, Edward A.G.; Shakil, Sarah; Shogren, Arial J.; Strauss, Jens; Tank, Suzanne E.; Thornton, Brett F.; Treharne, Rachael; Voigt, Carolina; Wright, Nancy; Yang, Yuanhe; Zarnetske, Jay P.; Zhang, Qiwen; Zolkos, Scott
Climate change is an existential threat to the vast global permafrost domain. The diverse human cultures, ecological communities, and biogeochemical cycles of this tenth of the planet depend on the persistence of frozen conditions. The complexity, immensity, and remoteness of permafrost ecosystems make it difficult to grasp how quickly things are changing and what can be done about it. Here, we summarize terrestrial and marine changes in the permafrost domain with an eye toward global policy. While many questions remain, we know that continued fossil fuel burning is incompatible with the continued existence of the permafrost domain as we know it. If we fail to protect permafrost ecosystems, the consequences for human rights, biosphere integrity, and global climate will be severe. The policy implications are clear: the faster we reduce human emissions and draw down atmospheric CO2, the more of the permafrost domain we can save. Emissions reduction targets must be strengthened and accompanied by support for local peoples to protect intact ecological communities and natural carbon sinks within the permafrost domain. Some proposed geoengineering interventions such as solar shading, surface albedo modification, and vegetation manipulations are unproven and may exacerbate environmental injustice without providing lasting protection. Conversely, astounding advances in renewable energy have reopened viable pathways to halve human greenhouse gas emissions by 2030 and effectively stop them well before 2050. We call on leaders, corporations, researchers, and citizens everywhere to acknowledge the global importance of the permafrost domain and work towards climate restoration and empowerment of Indigenous and immigrant communities in these regions.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.
Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process ows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.
The present study has used a variety of characterization techniques to determine the products and reaction pathways involved in the rechargeable Li-FeS2 system. We revisit both the initial lithiation and subsequent cycling of FeS2 employing an ionic liquid electrolyte to investigate the intermediate and final charge products formed under varying thermal conditions (room temperature to 100 °C). The detection of Li2S and hexagonal FeS as the intermediate phases in the initial lithiation and the electrochemical formation of greigite, Fe3S4, as a charge product in the rechargeable reaction differ significantly from previous reports. The conditions for Fe3S4 formation are shown to be dependent on both the temperature (∼60 °C) and the availability of sulfur to drive a FeS to Fe3S4 transformation. Upon further cycling, Fe3S4 transforms to a lower sulfur content iron sulfide phase, a process which coincides with the loss of sulfur based on the new reaction pathways established in this work. The connection between sulfur loss, capacity fade, and charge product composition highlights the critical need to retain sulfur in the active material upon cycling.
Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.
SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.
A method is developed to calculate the length into a sample to which a Frequency Domain Thermoreflectance (FDTR) measurement is sensitive. Sensing depth and sensing radius are defined as limiting cases for the spherically spreading FDTR measurement. A finite element model for FDTR measurements is developed in COMSOL multiphysics and used to calculate sensing depth and sensing radius for silicon and silicon dioxide samples for a variety of frequencies and laser spot sizes. The model is compared to experimental FDTR measurements. Design recommendations for sample thickness are made for experiments where semi-infinite sample depth is desirable. For measurements using a metal transducer layer, the recommended sample thickness is three thermal penetration depths, as calculated from the lowest measurement frequency.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.
The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.
As part of the development process, scaled testing of wave energy converter devices are necessary to prove a concept, study hydrodynamics, and validate control system approaches. Creating a low-cost, small, lightweight data acquisition system suitable for scaled testing is often a barrier for wave energy converter developers’ ability to test such devices. This paper outlines an open-source solution to these issues, which can be customized based on specific needs. Furthermore, this will help developers with limited resources along a path toward commercialization.
New and novel HPC platforms provide interesting challenges and opportunities. Analysis of these systems can provide a better understanding of both the specific platform being studied as well as large-scale systems in general. Arm is one such architecture that has been explored in HPC for several years, however little is still known about its viability for supporting large-scale production workloads in terms of system reliability. The Astra system at Sandia National Laboratories was the first public peta-FLOPS Arm-based system on the Top500 and has been successfully running production HPC applications for a couple of years. In this paper, we analyze memory failure data collected from Astra while the system was in production running unclassified applications. This analysis revealed several interesting contributions related to both the Arm platform and to HPC systems in general. First, we outline the number of components replaced due to reliability issues in standing-up this first-of-its-kind, large-scale HPC system. We show the distribution differences between correctable DRAM faults and errors on Astra, showing that, not properly accounting for faults can lead to erroneous conclusions. Additionally, we characterize DRAM faults on the system and show contrary to existing work that memory faults are uniformly distributed across CPU socket, DRAM column, bank and rack region, but are not uniform across node, DIMM rank, DIMM slot on the motherboard, and system rack: some racks, ranks and DIMM slots experience more faults than others. Similarly, we show the impact of temperature and power on DRAM correctable errors. Finally, we make a detailed comparison of results presented here with the positional affects found in several previous large-scale reliability studies. The results of this analysis provide valuable guidance to organizations standing-up first-in- class platforms in HPC, organizations using Arm in HPC, and the entire large-scale HPC community in general.
The decomposition of higher-order joint cumulant tensors of spatio-temporal data sets is useful in analyzing multi-variate non-Gaussian statistics with a wide variety of applications (e.g. anomaly detection, independent component analysis, dimensionality reduction). Computing the cumulant tensor often requires computing the joint moment tensor of the input data first, which is very expensive using a naïve algorithm. The current state-of-the-art algorithm takes advantage of the symmetric nature of a moment tensor by dividing it into smaller cubic tensor blocks and only computing the blocks with unique values and thus reducing computation. We propose a refactoring of this algorithm by posing its computation as matrix operations, specifically Khatri-Rao products and standard matrix multiplications. An analysis of the computational and cache complexity indicates significant performance savings due to the refactoring. Implementations of our refactored algorithm in Julia show speedups up to 10x over the reference algorithm in single processor experiments. We describe multiple levels of hierarchical parallelism inherent in the refactored algorithm, and present an implementation using an advanced programming model that shows similar speedups in experiments run on a GPU.
There is a need to understand materials exposed to overlapping extreme environments such as high temperature, radiation, or mechanical stress. When these stressors are combined there may be synergistic effects that enable unique microstructural evolution mechanisms to activate. Understanding of these mechanisms is necessary for the input and refinement of predictive models and critical for engineering of next generation materials. The basic physics and underlying mechanisms require advanced tools to be investigated. The in situ ion irradiation transmission electron microscope (I³TEM) is designed to explore these principles. To quantitatively probe the complex dynamic interactions in materials, careful preparation of samples and consideration of experimental design is required. Particular handling or preparation of samples can easily introduce damage or features that obfuscate the measurements. There is no one correct way to prepare a sample; however, many mistakes can be made. The most common errors and things to consider are highlighted within. The I³TEM has many adjustable variables and a large potential experimental space, therefore it is best to design experiments with a specific scientific question or questions in mind. Experiments have been performed on large number of sample geometries, material classes, and with many irradiation conditions. The following are a subset of examples that demonstrate unique in situ capabilities utilizing the I3TEM. Au nanoparticles prepared by drop casting have been used to investigate the effects of single ion strikes. Au thin films have been used in studies on the effects of multibeam irradiation on microstructure evolution. Zr films have been exposed to irradiation and mechanical tension to examine creep. Ag nanopillars were subjected to simultaneous high temperature, mechanical compression, and ion irradiation to study irradiation induced creep as well. These results impact fields including: structural materials, nuclear energy, energy storage, catalysis, and microelectronics in space environments.
Power Spectrum Analysis (PSA) is a Sandia-developed, non-intrusive, electrical technique that captures distinct frequency-domain signatures of microelectronics devices using an innovative, unconventional biasing scheme (off-normal biasing). PSA can identify subtle differences in devices and is applicable in various areas such as device screening, counterfeit identification, reliability assurance, and trust authentication. From October 2020 to April 2021, Sandia worked with entrepreneurs from a new start-up company, Chiplytics, to commercialize PSA technology through NNSA-sponsored FedTech Program. In September 2021, Sandia received funding through Covid-19 Technical Assistance Program (CTAP) to provide technical assistance to Chiplytics for commercialization. Under the CTAP Statement of Work, Sandia was tasked with providing technical assistance to Chiplytics in PSA pilot testing for Naval Surface Warfare Center (NSWC) at Crane and other pilot participants. Sandia was also tasked with assisting Chiplytics in hardware development and evaluation of Chiplytics prototype system.
The magneto-Rayleigh-Taylor instability (MRTI) plays an essential role in astrophysical systems and in magneto-inertial fusion, where it is known to be an important degradation mechanism of confinement and target performance. In this Letter, we show for the first time experimental evidence of mode mixing and the onset of an inverse-cascade process resulting from the nonlinear coupling of two discrete preseeded axial modes (400- and 550-μm wavelengths) on an Al liner that is magnetically imploded using the 20-MA, 100-ns rise-time Z Machine at Sandia National Laboratories. Four radiographs captured the temporal evolution of the MRTI. We introduce a novel unfold technique to analyze the experimental radiographs and compare the results to simulations and to a weakly nonlinear model. We find good quantitative agreement with simulations using the radiation magnetohydrodynamics code hydra. Spectral analysis of the MRTI time evolution obtained from the simulations shows evidence of harmonic generation, mode coupling, and the onset of an inverse-cascade process. The experiments provide a benchmark for future work on the MRTI and motivate the development of new analytical theories to better understand this instability.
The E3 transition in irradiated GaAs observed in deep level transient spectroscopy (DLTS) was recently discovered in Laplace-DLTS to encompass three distinct components. The component designated E3c was found to be metastable, reversibly bleached under minority carrier (hole) injection, with an introduction rate dependent upon Si doping density. It is shown through first-principles modeling that the E3c must be the intimate Si-vacancy pair, best described as a Si sitting in a divacancy Sivv. The bleached metastable state is enabled by a doubly site-shifting mechanism: Upon recharging, the defect undergoes a second site shift rather returning to its original E3c-active configuration via reversing the first site shift. Identification of this defect offers insights into the short-time annealing kinetics in irradiated GaAs.
Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for individual components in binary fluid mixtures. The ANNs were tested on an experimental database of 4328 self-diffusion constants from 131 mixtures containing 75 unique compounds. The presence of strong hydrogen bonding molecules may lead to clustering or dimerization resulting in non-linear diffusive behavior. To address this, self- and binary association energies were calculated for each molecule and mixture to provide information on intermolecular interaction strength and were used as input features to the ANN. An accurate, generalized ANN model was developed with an overall average absolute deviation of 4.1%. Forward input feature selection reveals the importance of critical properties and self-association energies along with other fluid properties. Additional ANNs were developed with subsets of the full input feature set to further investigate the impact of various properties on model performance. The results from two specific mixtures are discussed in additional detail: one providing an example of strong hydrogen bonding and the other an example of extreme pressure changes, with the ANN models predicting self-diffusion well in both cases.
A direct comparison between electron transparent transmission electron microscope (TEM) samples prepared with gallium (Ga) and xenon (Xe) focused ion beams (FIBs) is performed to determine if equivalent quality samples can be prepared with both ion species. We prepared samples using Ga FIB and Xe plasma focused ion beam (PFIB) while altering a variety of different deposition and milling parameters. The samples' final thicknesses were evaluated using STEM-EELS t/λ data. Using the Ga FIB sample as a standard, we compared the Xe PFIB samples to the standard and to each other. We show that although the Xe PFIB sample preparation technique is quite different from the Ga FIB technique, it is possible to produce high-quality, large area TEM samples with Xe PFIB. We also describe best practices for a Xe PFIB TEM sample preparation workflow to enable consistent success for any thoughtful FIB operator. For Xe PFIB, we show that a decision must be made between the ultimate sample thickness and the size of the electron transparent region.
In this work we present the concept of g'clipping', scheduling receive events for wireless transmissions only on receivers within some distance of the transmitter. Combined with spatial indexing, this technique enables faster simulation of large-scale wireless networks containing tens of thousands or even hundreds of thousands of wireless nodes. We detail our additions and changes to ns-3 to implement this feature, demonstrate how it yields a 2 × speedup for a complex 5G scenario with minimal impact on simulation fidelity, and show how under special circumstances a speedup of over 40 × is achievable while producing identical results.
This is an addendum to the Sierra/SolidMechanics 5.8 User’s Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State’s International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 5.8 User’s Guide should be referenced for most general descriptions of code capability and use.
X-ray tomography is capable of imaging the interior of objects in three dimensions non-invasively, with applications in biomedical imaging, materials science, electronic inspection, and other fields. The reconstruction process can be an ill-conditioned inverse problem, requiring regularization to obtain satisfactory results. Recently, deep learning has been adopted for tomographic reconstruction. Unlike iterative algorithms which require a distribution that is known a priori, deep reconstruction networks can learn a prior distribution through sampling the training distributions. In this work, we develop a Physics-assisted Generative Adversarial Network (PGAN), a two-step algorithm for tomographic reconstruction. In contrast to previous efforts, our PGAN utilizes maximum-likelihood estimates derived from the measurements to regularize the reconstruction with both known physics and the learned prior. Compared with methods with less physics assisting in training, PGAN can reduce the photon requirement with limited projection angles to achieve a given error rate. The advantages of using a physics-assisted learned prior in X-ray tomography may further enable low-photon nanoscale imaging.
Civil infrastructure is made primarily of concrete structures or components and therefore understanding durability and fracture behavior of concrete is of utmost importance. Concrete contains an interfacial transition zone (ITZ), a porous region surrounding the aggregates, that is often considered to be the weakest region in the concrete. The ITZ is poorly characterized and property estimates for the ITZ differ considerably. In this simulation study, representative concrete mesostructures are produced by packing coarse aggregates with realistic geometries into a mortar matrix. A meshless numerical method, peridynamics, is utilized to simulate the mechanical response including fracture under uniaxial compression and tension. The sensitivity of the stiffness and fracture toughness of the samples to the ITZ properties is computed, showing strong relationships between the ITZ properties and the effective modulus and effective yield strength of the concrete. These results provides insight into the influence of the poorly characterized ITZ on the stiffness and strength of concrete. This work showcases the applicability of peridynamics to concrete systems, matching experimental strength and modulus values. Additionally, relationships between the ITZ's mechanical properties and the overall concrete strength and stiffness are presented to enable future design decisions.
The effect of psychosocial factors after a disaster is well acknowledged. They can manifest in many ways including policy decisions, increased medical problems due to stress, avoidance behavior, and lessened trust. These can have direct economic effects, including increased medical costs, decreased tourism, de facto boycotts of product, and delated remedial actions due to delayed decision making.
This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSMs) from Orano (formerly Transnuclear Inc.) and provides a description of the surface characterization activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. The goal of the testing is to collect highly defensible and detailed dust deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion crack (CISCC) research. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the decay heat and thermal hydraulic environment. Test and heater design is supported by detailed computational fluid dynamics modeling. Instrumentation throughout the canister, storage module, and environment will provide extensive information about thermal-hydraulic behavior. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer a high-fidelity quantification of the conditions experienced in a harsh yet realistic environment. Functional testing of the finalized heater assemblies and test apparatus is set to begin in December 2022. The proposed delivery of the canisters to the host test site is June/July 2023, which is well ahead of when the AHSM installations would be completed.
Laplace, Thibault A.; Goldblum, Bethany L.; Brown, Joshua A.; Leblanc, Glenn; Li, Tianyue; Manfredi, Juan J.; Brubaker, Erik B.
Recoil nuclei produce high ionization and excitation densities in organic scintillators leading to reduced light yield via ionization quenching. To improve understanding of the relationship between organic scintillator specific luminescence and the characteristics of the recoil particle, this work evaluates proton and carbon light yield data using ionization quenching models over an energy range of tens of keV to several MeV for protons and 1-5 MeV for carbon ions. Previously-measured proton and carbon light yield data were examined for a variety of commercial and novel organic scintillating media: EJ-309, a liquid with pulse shape discrimination (PSD) properties; EJ-204, a fast plastic; EJ-276, a PSD-capable plastic; and a custom organic glass scintillator developed by Sandia National Laboratories. The canonical model of Birks did not adequately describe the ionization quenching behavior. Models proposed by Yoshida et al. and Voltz et al. provided a reasonable description of the proton light yield of a variety of organic scintillators over a broad energy range, but additional work is needed to extend the models to carbon ions. The impact of stopping power data was also investigated by comparing model predictions using SRIM and PSTAR/MSTAR libraries, and the results show a significant discrepancy for carbon ions. This work enhances understanding of ionization quenching and facilitates the accurate modeling of scintillator-based neutron detection systems relevant for medical physics, nuclear security and nonproliferation, and basic science studies.
Conversion cathodes represent a viable route to improve rechargeable Li+battery energy densities, but their poor electrochemical stability and power density have impeded their practical implementation. Here, we explore the impact cell fabrication, electrolyte interaction, and current density have on the electrochemical performance of FeS2/Li cells by deconvoluting the contributions of the various conversion and intercalation reactions to the overall capacity. By varying the slurry composition and applied pressure, we determine that the capacity loss is primarily due to the large volume changes during (de)lithiation, leading to a degradation of the conductive matrix. Through the application of an external pressure, the loss is minimized by maintaining the conductive matrix. We further determine that polysulfide loss can be minimized by increasing the current density (>C/10), thus reducing the sulfur formation period. Analysis of the kinetics determines that the conversion reactions are rate-limiting, specifically the formation of metallic iron at rates above C/8. While focused on FeS2, our findings on the influence of pressure, electrolyte interaction, and kinetics are broadly applicable to other conversion cathode systems.
The characterization of the three-dimensional arrangement of dislocations is important for many analyses in materials science. Dislocation tomography in transmission electron microscopy is conventionally accomplished through intensity-based reconstruction algorithms. Although such methods work successfully, a disadvantage is that they require many images to be collected over a large tilt range. Here, we present an alternative, semi-automated object-based approach that reduces the data collection requirements by drawing on the prior knowledge that dislocations are line objects. Our approach consists of three steps: (1) initial extraction of dislocation line objects from the individual frames, (2) alignment and matching of these objects across the frames in the tilt series, and (3) tomographic reconstruction to determine the full three-dimensional configuration of the dislocations. Drawing on innovations in graph theory, we employ a node-line segment representation for the dislocation lines and a novel arc-length mapping scheme to relate the dislocations to each other across the images in the tilt series. We demonstrate the method for a dataset collected from a dislocation network imaged by diffraction-contrast scanning transmission electron microscopy. Based on these results and a detailed uncertainty analysis for the algorithm, we discuss opportunities for optimizing data collection and further automating the method.
Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this work we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.
Gao, Xian; Knueven, Bernard; Siirola, John D.; Miller, David C.; Dowling, Alexander W.
Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. In this paper, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).
Rechargeable Zn metal batteries (RZMBs) may provide a more sustainable and lower-cost alternative to established battery technologies in meeting energy storage applications of the future. However, the most promising electrolytes for RZMBs are generally aqueous and require high concentrations of salt(s) to bring efficiencies toward commercially viable levels and mitigate water-originated parasitic reactions including hydrogen evolution and corrosion. Electrolytes based on nonaqueous solvents are promising for avoiding these issues, but full cell performance demonstrations with solvents other than water have been very limited. To address these challenges, we investigated MeOH as an alternative electrolyte solvent. These MeOH-based electrolytes exhibited exceptional Zn reversibility over a wide temperature range, with a Coulombic efficiency > 99.5% at 50% Zn utilization without cell short-circuit behavior for > 1,800 h. More important, this remarkable performance translates well to Zn jj metal-free organic cathode full cells, supporting < 6% capacity decay after > 800 cycles at 240 °C.
This quick note outlines what we found after our conversion with you and your team. As suggested, we loaded 1547-2003 source requirements document (SRD) and then went back and loaded 1547-2018 SRD. This did result in implementing the new 1547-2018 settings. This short report focuses on the frequency-watt function and shows a couple of screen shots of the parameter settings via the Mojave HMI interface and plots of the results of the inverter with FW function enabled in both default and most aggressive settings response to frequency events. The first screen shot shows the 1547-2018 selected after selecting 1547-2003.
We propose a new scheme for simulation of collisions with multiple possible outcomes in variable-weight DSMC computations. The scheme is applied to a 0-D ionization rate coefficient computation, and 1-D electrical breakdown simulation. We show that the scheme offers a significant (up to an order of magnitude) improvement in the level of stochastic noise over the usual acceptance-rejection algorithm, even when controlling for the slight additional computational costs. Furthermore, the benefits and performance of the scheme are analyzed in detail, and possible extensions are proposed.
Attia, Peter M.; Bills, Alexander; Brosa Planella, Ferran; Dechent, Philipp; Dos Reis, Goncalo; Dubarry, Matthieu; Gasper, Paul; Gilchrist, Richard; Greenbank, Samuel; Howey, David; Liu, Ouyang; Khoo, Edwin; Preger, Yuliya P.; Soni, Abhishek; Stefanopoulou, Anna G.; Sulzer, Valentin
Lithium-ion batteries can last many years but sometimes exhibit rapid, nonlinear degradation that severely limits battery lifetime. In this work, we review prior work on “knees” in lithium-ion battery aging trajectories. We first review definitions for knees and three classes of “internal state trajectories” (termed snowball, hidden, and threshold trajectories) that can cause a knee. We then discuss six knee “pathways”, including lithium plating, electrode saturation, resistance growth, electrolyte and additive depletion, percolation-limited connectivity, and mechanical deformation—some of which have internal state trajectories with signals that are electrochemically undetectable. Additionally, we also identify key design and usage sensitivities for knees. Finally, we discuss challenges and opportunities for knee modeling and prediction. Our findings illustrate the complexity and subtlety of lithium-ion battery degradation and can aid both academic and industrial efforts to improve battery lifetime.
Thermodynamic modeling has been used to predict chemical compositions of brines formed by the deliquescence of sea salt aerosols. Representative brines have been mixed, and physical and chemical properties have been measured over a range of temperatures. Brine properties are discussed in terms of atmospheric corrosion of austenitic stainless steel, using spent nuclear fuel dry storage canisters as an example. After initial loading with spent fuel, during dry storage, the canisters cool over time, leading to increased surface relative humidities and evolving brine chemistries and properties. These parameters affect corrosion kinetics and damage distributions, and may offer important constraints on the expected timing, rate, and long-term impacts of canister corrosion.
Here, we report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3$\textit{d}$ orbitals are either in the core ($d^0$) or explicitly in the valence set ($d^{10}$). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability–a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the $d^0$ and $d^{10}$ pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the $d^0$ lattice constant with experiment suggests that the more computationally economical $d^0$ pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.
Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ∼80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1-2 orders of magnitude and decreased the modulus by 2-3 orders of magnitude due to surface hydroxylation. The strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.
Michael, Joseph R.; Perry, Daniel L.; Cummings, Damion P.; Walraven, Jeremy A.; Jordan, Matthew B.
Indium (In) and other low melting point metals are used as interconnects in a variety of hybridized circuits and a full understanding of the metallurgy of these interconnects is important to the reliability and performance of the devices. This paper shows that room temperature focused ion beam (FIB) preparation of cross-sections, using Ga+ or Xe+ can result in artifacts that obscure the true In microbump structure. The use of modified milling strategies to minimize the increased local sample temperature are shown to produce cross-sections that are representative of the In bump microstructure in some sample configurations. Furthermore, cooling of the sample to cryogenic temperatures is shown to reliably eliminate artifacts in FIB prepared cross-sections of In bumps allowing the true bump microstructure to be observed.
Yip, Ho L.; Srna, Ales S.; Wehrfritz, Armin; Kook, Sanghoon; Hawkes, Evatt R.; Chan, Qing N.
This study examines the flame evolution of autoigniting H2 jets with high-speed schlieren and OH∗ chemiluminescence optical methods in a constant-volume combustion chamber over a wide range of simulated compression-ignition engine conditions. Parametric variations include the injector nozzle orifice diameter (0.31–0.83 mm), injection reservoir pressure (100–200 bar), ambient temperature (1000–1140 K), density (12.5–24 kg/m3) and O2 concentration (10–21 vol.%). The jet ignition delay was found to be highly sensitive to changes in ambient temperature while all other parameter variations resulted in minor ignition delay changes. Optical imaging reveals that in most cases, the reaction front of the H2 jet initiates from a localised kernel, before engulfing the entire jet volume downstream and recessing towards the nozzle. The flames attach to the nozzle, except at the lowest ambient oxygen condition of 10 vol.% O2 for which a lifted flame is observed. The H2 diffusion flame length shows a dependence on both the mass flow rate and the level of O2 entrainment that follows the same correlations as previously established for atmospheric H2 jet flames.
Multiferroic materials are an interesting functional material family combining two ferroic orderings, e.g., ferroelectric and ferromagnetic orderings, or ferroelectric and antiferromagnetic orderings, and find various device applications, such as spintronics, multiferroic tunnel junctions, etc. Coupling multiferroic materials with plasmonic nanostructures offers great potential for optical-based switching in these devices. Here, we report a novel nanocomposite system consisting of layered Bi1.25AlMnO3.25 (BAMO) as a multiferroic matrix and well dispersed plasmonic Au nanoparticles (NPs) and demonstrate that the Au nanoparticle morphology and the nanocomposite properties can be effectively tuned. Specifically, the Au particle size can be tuned from 6.82 nm to 31.59 nm and the 6.82 nm one presents the optimum ferroelectric and ferromagnetic properties and plasmonic properties. Besides the room temperature multiferroic properties, the BAMO-Au nanocomposite system presents other unique functionalities including localized surface plasmon resonance (LSPR), hyperbolicity in the visible region, and magneto-optical coupling, which can all be effectively tailored through morphology tuning. This study demonstrates the feasibility of coupling single phase multiferroic oxides with plasmonic metals for complex nanocomposite designs towards optically switchable spintronics and other memory devices.
Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.
Once limited to chain-growth polymerizations, fine control over polymerization-induced phase separation (PIPS) has recently been demonstrated in rubber-toughened thermoset materials formed through step-growth polymerizations. The domain length scales of these thermoset materials can be elegantly tuned by utilizing a binary mixture of curing agents (CAs) that individually yield disparate morphologies. Importantly, varying the composition of the binary mixture affects characteristics of the materials such as glass transition temperature and tensile behavior. Here, we establish a full phase diagram of PIPS in a rubber-toughened epoxy system tuned by a binary CA mixture to provide a robust framework of phase behaviour. X-Ray scattering in situ and post-PIPS is employed to elucidate the PIPS mechanism whereby an initial polymerization-induced compositional fluctuation causes nanoscale phase separation of rubber and epoxy components prior to local chain crosslinking and potential macrophase separation. We further demonstrate the universality of this approach by alternatively employing binary epoxy or binary rubber mixtures to achieve broad variations in morphology and glass transitions.
Edgell, Dana H.; Colaitis, Arnaud; Guardalben, Mark J.; Kalb, Adam; Katz, Joe; Kwiatkowski, Joe; Mannion, Owen M.; Shvydky, Alex; Stoeckl, Christian; Turnbull, David; Froula, Dustin H.
Three dimensional polarization-dependent CBET plus beam balance and pointing are required to model nonuniformity in direct-drive implosions on OMEGA. Uniform laser energy absorption is essential for successful laser-direct-drive inertial confinement fusion but a growing body of evidence suggests OMEGA implosions are more asymmetric than predictions. By measuring the intensity and polarization of light scattered from individual beams, we have identified OMEGA’s polarization smoothing via distributed polarization rotators (DPR’s) as one previously unrealized source of nonuniformity. Polarization-dependent CBET along with beam energy balance, and beam pointing require three-dimensional modeling. Laser absorption mode 1 predictions from a fully three-dimensional CBET model correlate well with the observed direction of the core flow.
This study presents component-level testing of carbon fiber sandwich beams and the effect of carbon fiber material nonlinearity in its strain response in bending. A simple material model is presented and validated that accurately captures the carbon fiber longitudinal nonlinearity in both the tensile and compressive response. This material model is implemented in a finite element model of the BAR-DRC reference wind blade, a downwind 100-meter rotor blade, and the effects of the nonlinearity on ultimate limit states of the blade are analyzed. The material nonlinearity has negligible effect on the deflection, and material failure predictions. The buckling analysis revealed significant reductions in buckling load factor in the controlling flap direction caused by the material nonlinearity, revealing the importance of including this material model for buckling analyses of wind blade with carbon fiber reinforced spar caps.
This study presents component-level testing of carbon fiber sandwich beams and the effect of carbon fiber material nonlinearity in its strain response in bending. A simple material model is presented and validated that accurately captures the carbon fiber longitudinal nonlinearity in both the tensile and compressive response. This material model is implemented in a finite element model of the BAR-DRC reference wind blade, a downwind 100-meter rotor blade, and the effects of the nonlinearity on ultimate limit states of the blade are analyzed. The material nonlinearity has negligible effect on the deflection, and material failure predictions. The buckling analysis revealed significant reductions in buckling load factor in the controlling flap direction caused by the material nonlinearity, revealing the importance of including this material model for buckling analyses of wind blade with carbon fiber reinforced spar caps.
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.
Rezaul Karim, Mohammad; Narasimhachary, Santosh; Radaelli, Francesco; Amann, Christian; Dayal, Kaushik; Silling, Stewart A.; Germann, Timothy C.
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.
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.
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.
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.
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.
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.
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.
Keefe, Kevin; Alhajaji, Hassam; Brubaker, Erik B.; Druetzler, Andrew; Learned, John; Maggi, Paul; Manfredi, Juan J.; Nishimura, Kurtis; Souza, Bejamin P.; Steele, John T.; Sweany, Melinda; Takahashi, Eric
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.
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.