The magnetic-Rayleigh-Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh-Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.
This Technology Commercialization Fund (TCF) project was a complete success. As described in our original TCF proposal, TRL1 through TRL5 development of Sandia Cooler technology was carried out under EERE BTO and Sandia LDRD funding. The objective of this TCF project was for Sandia National Labs and commercial partner Wakefield-Vette to convert this TRL5 DOE technology to a market-ready, cost-competitive, electronics thermal management product (TRL 8). A further objective was to target an energy sector application with high-visibility that can effectively serve as an advertisement for DOE's Sandia Cooler technology.
We report measurements of K-shell fluorescence lines induced by fast electrons in ramp-compressed Co targets. The fluorescence emission was stimulated by fast electrons generated through short-pulse laser-solid interaction with an Al target layer. Compression up to 2.1× solid density was achieved while maintaining temperatures well below the Fermi energy, effectively removing the thermal effects from consideration. We observed small but unambiguous redshifts in the Kβ fluorescence line relative to unshifted Cu Kα. Redshifts up to 2.6 eV were found to increase with compression and to be consistent with predictions from self-consistent models based on density-functional theory.
Scott, Ethan A.; Hattar, Khalid M.; Braun, Jeffrey L.; Rost, Christina M.; Gaskins, John T.; Bai, Tingyu; Wang, Yekan; Ganski, Claire; Goorsky, Mark; Hopkins, Patrick E.
Despite the exceptional thermal and mechanical functionalities of diamond, its superlative properties are highly subject to the presence of point defects, dislocations, and interfaces. In this study, polycrystalline diamond is ion implanted with C3+, N3+, and O3+ ions at an energy of 16.5 MeV, producing an amorphous layer at the projected range and a damaged crystalline region between the surface and amorphous layer. Using time-domain thermoreflectance in combination with thermal penetration depth calculations based upon the multilayer heat diffusion equation, it is determined that reductions in the thermal conductivity can span nearly two orders of magnitude while still maintaining a polycrystalline structure within the regions thermally probed. Dynamical diffraction simulations of high-resolution x-ray diffraction measurements demonstrate the formation of a strained layer localized at the end of range, with much lower levels of strain near the surface. Furthermore, within the polycrystalline region above the amorphous layer, the average number of displacements-per-atom from the ion irradiation is found to be <1%, with mass impurity concentrations much less than 1%. These low defect concentrations within the thermally probed region demonstrate the remarkably large impact that dilute levels of defects from the ion implantation can have on the thermal conductivity of diamond.
Hexagonal and cubic crystals contain paired sets of internal planes that reflect the linearly polarized components of certain x rays into two separate, perpendicular directions. For the cubic crystals, two distinct crystal orientations provide the same polarization-splitting geometry. One of the orientations may have advantages for plasma spectroscopy by suppressing unwanted reflections. This paper demonstrates the two orientations with a germanium crystal and K characteristic lines from copper and zirconium.
Sandia National Laboratories continued evaluation of total system performance assessment (TSPA) computing systems for the previously considered Yucca Mountain Project. This was done to maintain the operational readiness of the computing infrastructure (computer hardware and software) and knowledge capability for total system performance assessment) type analysis, as directed by the National Nuclear Security Administration (NNSA), DOE 2010. The FY19 task included continued operation of the cluster; maintenance of the TSPA-LA models (with Gold Sim 9.60.300); preliminary assessment of the status of the Infiltration Model (a process model that feeds the TSPA-LA). In addition, precautionary actions were needed to extend the life of the cluster hardware. To do that, three new nodes were added to the cluster. In the event any of the original nodes fail they will be replaced with the new nodes, thereby maintaining the core capability. The 2014 cluster and supporting software systems are currently fully operational to support TSPA-LA type analysis.
A principal performance-enabling, or performance-limiting, component of Ground-Moving-Target-Indication (GMTI) radar systems is the antenna. Undesired clutter leakage into antenna sidelobes can be particularly problematic, generating undesired false alarms. GMTI system antennas can be designed with characteristics and features to allow discriminating and depressing/suppressing problematic sidelobe leakage of clutter and other undesired signals. We offer analysis and design guidelines for doing so.
As penetration of converter interfaced generators (CIGs) increases, the need for CIG frequency control participation increases. Traditionally, research in this area has been performed using positive sequence simulation software, which provides voltage magnitude and phase measurements, but not point-on-wave (POW) measurements. This means that the effect of frequency estimation algorithms cannot be accurately modeled, especially when the voltage waveform is distorted by faults or load connection events. This report serves as a user manual for an electromagnetic transient simulation testbed, which allows for accurate modeling of frequency estimation and control techniques.
Particle-laden turbulent flows subject to radiative heating are relevant in many applications, for example concentrated solar power receivers. Efficient and accurate simulations provide valuable insights and enable optimization of such systems. However, as there are many uncertainties inherent in such flows, uncertainty quantification is fundamental to improve the predictive capabilities of the numerical simulations. For large-scale, multi-physics problems exhibiting high-dimensional uncertainty, characterizing the stochastic solution presents a significant computational challenge as most strategies require a large number of high-fidelity solves. This requirement might result in an infeasible number of simulations when a typical converged high-fidelity simulation requires intensive computational resources. To reduce the cost of quantifying high-dimensional uncertainties, we investigate the application of a non-intrusive, bi-fidelity approximation to estimate statistics of quantities of interest associated with an irradiated particle-laden turbulent flow. This method exploits the low-rank structure of the solution to accelerate the stochastic sampling and approximation processes by means of cheaper-to-run, lower fidelity representations. The application of this bi-fidelity approximation results in accurate estimates of the quantities of interest statistics, while requiring a small number of high-fidelity model evaluations.
We expand the second-order fluid-structure coupling scheme of Farhat et al. (1998, "Load and Motion Transfer Algorithms for 19 Fluid/Structure Interaction Problems With Non-Matching Discrete Interfaces: Momentum and Energy Conservation, Optimal Discretization and Application to Aeroelasticity,"Comput. Methods Appl. Mech. Eng., 157(1-2), pp. 95-114; 2006, "Provably Second-Order Time-Accurate Loosely-Coupled Solution Algorithms for Transient Nonlinear Computational Aeroelasticity,"Comput. Methods Appl. Mech. Eng., 195(17), pp. 1973-2001) to structural acoustics. The staggered structural acoustics solution method is demonstrated to be second-order accurate in time, and numerical results are compared to a monolithically coupled system. The partitioned coupling method is implemented in the Sierra Mechanics software suite, allowing for the loose coupling of time domain acoustics in sierra/sd to structural dynamics (sierra/sd) or solid mechanics (sierra/sm). The coupling is demonstrated to work for nonconforming meshes. Results are verified for a one-dimensional piston, and the staggered and monolithic results are compared to an exact solution. Huang, H. (1969, "Transient Interaction of Plane Acoustic Waves With a Spherical Elastic Shell,"J. Acoust. Soc. Am., 45(3), pp. 661-670) sphere scattering problem with a spherically spreading acoustic load demonstrates parallel capability on a complex problem. Our numerical results compare well for a bronze plate submerged in water and sinusoidally excited (Fahnline and Shepherd, 2017, "Transient Finite Element/Equivalent Sources Using Direct Coupling and Treating the Acoustic Coupling Matrix as Sparse,"J. Acoust. Soc. Am., 142(2), pp. 1011-1024).
Sandia National Laboratories (SNL) has hosted the International Training Course on the Physical Protection of Nuclear Materials and Nuclear Facilities since 1978. This course is the flagship training course of the International Atomic Energy Agency (IAEA). On behalf of the National Nuclear Security Administration (NNSA), SNL manages, develops, and coordinates all course materials, and works closely with the IAEA to arrange all logistical details for the course. ITC-28 incorporated several new approaches based on feedback and experience with ITC-27 and earlier versions of the course. For ITC-28, an addition to the Integrated Security Facility (ISF) at SNL was a mock reactor hall with a mock reactor pool. Other facilities at the ISF include a mock processing facility, material receiving area, and central alarm station. The physical protection system at the ISF—an area that formerly housed Category I nuclear material—provides many opportunities for hands-on, real world training in the design and evaluation of a physical protection system (PPS). This document provides a brief description of ITC-28, including a summary of lessons learned and key recommendations for future development efforts.
Tungsten samples were exposed to He plasmas generated by an RF source (Γi= 3.5 x 1016 He cm-2s-1, ion energy = 92 eV.) The range of exposure conditions selected here is conducive to the growth of nearsurface He bubbles, and at higher fluence, the formation of W nanotendrils ranging between 50 — 100 nm in diameter. The evolution of these surface features was probed using a fixed-angle ellipsometer (280 — 1000 nm wavelength range) with direct line-of-sight to the sample. Over the parameter space explored here, changes in the two angles (p, S) that define the polarization of the reflected light followed a distinct trajectory with increasing plasma fluence. Ex-situ ellipsometry of 22 additional tungsten specimens tested at a wide range of plasma fluences and temperatures mapped onto these in-situ results well. We used helium ion microscopy and focused ion beam profiling to provide a direct calibration of the ellipsometry measurements. Our results indicate that for a reproducible process such as the growth helium-induced surface morphologies, ellipsometry is a practical in-situ diagnostic to study how fusion plasmas modify materials. To study more general effects of plasmas on surfaces, including co-deposition and sputtering, different approaches to modelling the optical properties of the exposed surfaces are also considered.
The V28 containment vessel was procured by the US Army Recovered Chemical Material Directorate (RCMD) as a replacement vessel for use on the P2 Explosive Destruction Systems. It is the fourth EDS vessel to be fabricated under Code Case 2564 of the ASME Boiler and Pressure Vessel Code, which provides rules for the design of impulsively loaded vessels. The explosive rating for the vessel, based on the Code Case, is nine (9) pounds TNT- equivalent for up to 637 detonations. This report documents the results of explosive tests that were done on the vessel at Sandia National Laboratories in Albuquerque New Mexico to qualify the vessel for explosive use. The primary qualification test consisted of six, 1.5 pounds charges of Composition C-4 (equivalent to 11.25 pounds TNT) distributed around the vessel in accordance with the User Design Specification. This test was repeated due to a lack of proper clamp settings. Two additional tests using less explosive were performed, one identical in configuration to a test performed in the V27 qualification series as a baseline for comparison, and one where the separation distance of the charges was increased to extend the V27 analysis of distributed load effects on the P2 vessel. All vessel acceptance criteria were met.
Performance results collected for EMPIRE plasma simulations on ATS-1, ATS-2 and Astra are shown, showing portability and areas for improvement during the FY20 L1 milestone.
The cover letter and enclosure provide the Endorsement SM-GOV Analytical Results Summary Table for the ABCWUA split sampling activities that occurred January 22, 2019, through January 29, 2019.
This letter represents the five-day written notification for elevated fluoride levels at the Building 858N, Acid Waste Neutralization (AWN) System in accordance with the ABCWUA Sewer Use and Wastewater Ordinance § 3-6-5(B), Reports of Potential Problems. This AWN System is covered under the ABCWUA Industrial Wastewater Permit No. 2069G. The immediate notification to the ABCWUA occurred at 21:42 on August 1, 2019.
This January 2020 monthly report is intended to communicate the status of North Slope Atmospheric Radiation Measurement (ARM) facilities managed by Sandia National Labs.
Sandia National Laboratories has tested and evaluated two variations of a new model of infrasound sensor, the Hyperion 5313A and 5119A. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, full scale, self-noise, dynamic range, response, passband, sensitivity variation due to changes in static pressure and temperature, and sensitivity to vertical acceleration. The Hyperion infrasound sensors are being evaluated for potential use in the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO).
We summarize here the double exponential and inverse double exponential approximations for two common EMP waveforms, the Bell Laboratories (Bell Labs) and the International- military standard (IEC-MIL-STD). Both models have been used frequently due to their relatively easy analytical expressions for both the time domain waveforms and their associated frequency domain spectra.
Characterizing the electromagnetic environment created by the recently upgraded HERMES III pulser is the primary objective addressed by this report. The pulser upgrade design was intended not only to provide higher electric field and magnetic field illumination amplitudes but also to provide a more collimated beam, focused along the center-line of the courtyard. For high-voltage pulsers to be useful as qualification systems, they must have the requisite peak values and uniform, predictable peak field spatial patterns. To assess the peak electric field values and uniformity of the spatial pattern of the electric fields for the upgraded HERMES system, a 1 6 full-power shot illumination of 9 electric and 9 magnetic field sensors in the courtyard of the HERMES III facility was conducted. The data acquired from this test showed an increase an increase in electric field amplitude, but still an asymmetry in the spatial distributions.
Photovoltaic energy prediction models include functions or modifiers to account for sun angle reflection losses. These functions may be known interchangeably as Angle of Incidence (AOI) or Incident Angle Modifier (IAM). While standards exist, there is no universally accepted single best practice for developing these functions. They can be generated through characterization of representative modules or single cells, in natural sunlight or indoors using simulated light sources. Repeatability of measurements and the viability of cross-laboratory comparisons are critical to confidence in validation of both methods. To investigate the differences between methods and labs, The Technical University of Denmark (DTU) initiated an international round-robin test comparison between several key test labs with AOI measurement capability. A total of six minimodules were provided in three different cell/interconnect/backsheet combinations. Sandia characterized these minimodules using methods developed over two decades specifically for the outdoor characterization of full-size photovoltaic modules. This report documents the characterization results, summarizes key observations and tabulates the processed data for comparison to results provided by other characterization labs.
This project is focused on developing advanced combustion strategies for mixing-controlled compressionignition (i.e., diesel-cycle) engines that are synergistic with renewable and/or unconventional fuels in a manner that enhances domestic energy security, economic competitiveness, and environmental quality. During this reporting period, the focus was on ducted fuel injection (DFI), a technology that differs from conventional diesel combustion (CDC) in that it involves injecting fuel along the axis of one or more small cylindrical ducts within the combustion chamber. Each duct performs a function similar to the tube on a Bunsen burner, helping to premix the fuel with the charge-gas before ignition, creating a stable flame that forms little to no soot. The purpose of the work conducted during Fiscal Year (FY) 2019 was to begin determining the extent to which the use of oxygenated fuels, when combined with DFI and charge-gas dilution, can simultaneously lower the soot and nitrogen-oxides (N0x) emissions from mixing-controlled compression-ignition engines, and what the corresponding impacts on other regulated emissions and efficiency are likely to be.
This report reviews and offers recommendations from Sandia National transportation of hazardous materials in the U.S. The risk criteria should be used with the results of a quantitative risk assessment (QRA) in risk acceptance decision-making. The QRA for transportation is fundamentally the same as a fixed facility. However, there are differences in calculations of both the probabilities of occurrence and location of hazards. Involuntary individual fatality risk is recommended to be acceptable for annual probabilities of less than 3 x 10-7 for any population, including vulnerable populations, and may be considered acceptable at the regulators discretion for non-sensitive/non-vulnerable populations if less than 5 x 10-5 and demonstrated to be as low as reasonably practicable (ALARP). Societal risk is recommended to be acceptable if the annual frequency of events that would result in N or more fatalities is less than 10-5/N events per year and may be considered acceptable at the regulators discretion if less than 10-3/N events per year and demonstrated to be ALARP. These criteria should be applied to the societal risk over the entire transportation route, not normalized per-distance. These values are adapted from the National Fire Protection Association (NFPA) 59A, a U.S. and international standard for liquefied natural gas (LNG) facility siting.
A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.
On February 11th, 2020 FRMAC participated in the 2020 RadResponder nationwide drill involving survey data review and approval. During this drill, dozens of organizations from across the country participated in uploading simulated survey results to their own RadResponder event and reviewing those data given some "best-practice" guidelines. FRMAC participated in this drill to learn how this process is carried out in RadResponder and to determine some areas of improvement that will make FRMACs use of RadResponder more effective and efficient. As drill participants, the four FRMAC staff each set up their own RadResponder events so they could work independently to work through each step. Each were given an import file with 60 survey measurements to import and review. At the end of the drill a hotwash was held where users were able to ask questions and receive feedback from Chainbridge and FRMAC participants. In summary, several areas for improvement have been identified and will be categorized as items that need immediate attention (Bugs), areas for improvement in the near-term (change requests), and ideas on how the overall process can be improved (feature requests). The items presented in this report by no means reflect all the required changes FRMAC Monitoring and Sampling and Assessment will need to carry out their procedures in RadResponder but do represent a step in the right direction. A full review and exercise of the FRMAC process will need to be conducted with more subject matter experts to discover all the required features FRMAC needs to carry out their mission.
Energy storage technologies are positioned to play a substantial role in power delivery systems. They are being touted as an effective new resource to maintain reliability and allow for increased penetration of renewable energy. However, due to their relative infancy, there is a lack of knowledge on how these resources truly operate over time. Data analysis can help ascertain the operational and performance characteristics of these emerging technologies. Rigorous testing and data analysis are important for all stakeholders to ensure a safe, reliable system that performs predictably on a macro level. Standardizing testing and analysis approaches to verifying the performance of energy storage devices, equipment, and systems when integrating them into the grid will improve the understanding and benefit of energy storage over time from technical and economic vantage points. Demonstrating the life-cycle value and capabilities of energy storage systems begins with the data the provider supplies for analysis. After review of energy storage data received from several providers, it has become clear that some of these data are inconsistent and incomplete, raising the question of their efficacy for robust analysis. This report reviews and proposes general guidelines such as sampling rates and data points that providers must supply for robust data analysis to take place. Consistent guidelines are the basis of the proper protocol and ensuing standards to (a) reduce the time it takes data to reach those who are providing analysis; (b) allow them to better understand the energy storage installations; and (c) enable them to provide high-quality analysis of the installations. This report is intended to serve as a starting point for what data points should be provided when monitoring. As battery technologies continue to advance and the industry expands, this report will be updated to remain current.
A Talbot-Lau X-ray Deflectometer (TXD) was implemented in the OMEGA EP laser facility to characterize the evolution of an irradiated foil ablation front by mapping electron densities >1022 cm-3 by means of Moiré deflectometry. The experiment used a short-pulse laser (30-100 J, 10 ps) and a foil copper target as an X-ray backlighter source. In the first experimental tests performed to benchmark the diagnostic platform, grating survival was demonstrated and X-ray backlighter laser parameters that deliver Moiré images were described. The necessary modifications to accurately probe the ablation front through TXD using the EP-TXD diagnostic platform are discussed.
JEDEC matrix trays are an industry standard for the safe handling, transport, and storage of electronics. Manufacturing tolerances of the trays (typical for injection-molded components) forces the creation of an envelope around components rather than a secure interface to prevent the trays themselves from damaging the electronics during use. However, this allows relative motion between the components, another potential damage source. This paper covers the design of tray features (flexures) than enable intentional, tuned contact preload to prevent relative motion, maximizing component safety while keeping the JEDEC tray form-factor. The target preload is balanced between gravity and shipment forces that could be experience during adverse handling or transport to prevent loss of contact or slip and keeping the load itself low enough to prevent component damage. Typical features/dimensions for JECEC trays were employed to maximize utility and minimize tray fabrication difficulty. The design also increases usability and component safety by making the chips visible while still sandwiched between trays. The point design described here exemplifies a simple, easy-to-manufacture tray-matrix flexure that significantly improves the security of supported components
Gamma-ray dose rates computed by two radiation source model-based methods and five measured spectrum-based methods are compared. With the exception of a Detective DX-100, which reported dose rates that are about 25% less than the mean, dose rates are generally within 10% of mean values. Methods for computation of dose rates employed by the Gamma Detector Response and Analysis Software (GADRAS) are described and critically evaluated. Results are presented for a feasibility evaluation that seeks to use spectra recorded inside vehicles to estimate the minimum dose outside the vehicles.
Three codes are used for comparisons in this report to evaluate the adequacy of MACCS in the nearfield, AERMOD, ARCON96 and QUIC. Test cases were developed to give a broad range of weather conditions, building dimensions and plume buoyancy. Based on the comparisons of MACCS with AERMOD, ARCON96 and QUIC across the test cases, the following observations are made: MACCS calculations configured with point-source, ground-level, nonbuoyant plumes provide nearfield results that bound the centerline, ground-level air concentrations from AERMOD, ARCON96, and QUIC. MACCS calculations with ground-level, nonbuoyant plumes that include the effects of the building wake (area source) provide nearfield results that bound the results from AERMOD and QUIC and the results from ARCON96 at distances >250 m. If using a point-source is too conservative and it is desired to bound the results from all three codes, another alternative is to use area source parameters in MACCS that are less than the standard values, i.e., an area source intermediate between the standard recommendation and a point source. All these options provide results from MACCS that are bounding for the test cases evaluated. Based on these observations, it appears that MACCS is adequate for use in nearfield calculations, given the correct parameterization.
An MHD mode with a frequency of <10 kHz has been identified near the inner strike point from various diagnostics, i.e., divertor Langmuir probes, magnetics sensors, and interferometers, but does not appear in the upstream and core diagnostics. This MHD mode is associated with magnetic oscillations of ≳5 G, has a long wavelength in the toroidal direction with toroidal mode number n = 1, but is localized in a narrow radial region of a few cm. The mode appears when the outer strike point remains attached and the inner strike point nearly detaches, grows with increasing density, and eventually weakens and vanishes as the outer strike point detaches. This mode results in particle flux with an order of magnitude higher than the background plasmas near the inner strike point. The mode characteristics are consistent with the Current-Convective-Instability theory prediction. Initial simulations based on experimental input have found oscillations with similar frequencies but weaker amplitude.
Weiss, Lukas; Wensing, Michael; Hwang, Joonsik; Pickett, Lyle M.; Skeen, Scott A.
Abstract: The method for direct injection of fuel in the cylinder of an IC engines is important to high-efficiency and low-emission performance. Optical spray diagnostics plays an important role in understanding plume movement and interaction for multi-hole injectors, and providing baseline understanding used for computational optimization of fuel delivery. Traditional planar or line-of-sight diagnostics fail to capture the liquid distribution because of optical thickness concerns. This work proposes a high-speed (67 kHz) extinction imaging technique at various injector rotations coupled to computed tomography (CT) for time-resolved reconstruction of liquid volume fraction in three dimensions. The number of views selected and processing were based on synthetic (modeled) liquid volume fraction data where extinction and CT adequately reconstructed each plume. The exercise showed that for an 8-hole, symmetric-design injector (ECN Spray G), only three different views are enough to reproduce the direction of each plume, and particularly the mean plume direction. Therefore, the number of views was minimized for experiments to save expense. Measurements applying this limited-view technique confirm plume–plume variations also detected with mechanical patternation, while providing better spatial and temporal resolution than achieved previously. Uncertainties due to the limited view within pressurized spray chambers, the droplet size, and optically thick regions are discussed. Graphic abstract: [Figure not available: see fulltext.].
We report on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (Itotal ∼5 × 1014 W/cm2, 1.5 ns square pulse, and Etotal ∼16.2 kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.
Collisions that induce melting and vaporization can have a substantial effect on the thermal and geochemical evolution of planets. However, the thermodynamics of major minerals are not well known at the extreme conditions attained during planet formation. We obtained new data at the Sandia Z Machine and use published thermodynamic data for the major mineral forsterite (Mg2SiO4) to calculate the specific entropy in the liquid region of the principal Hugoniot. We use our calculated specific entropy of shocked forsterite, and revised entropies for shocked silica, to determine the critical impact velocities for melting or vaporization upon decompression from the shocked state to 1 bar and the triple points, which are near the pressures of the solar nebula. We also demonstrate the importance of the initial temperature on the criteria for vaporization. Applying these results to N-body simulations of terrestrial planet formation, we find that up to 20% to 40% of the total system mass is processed through collisions with velocities that exceed the criteria for incipient vaporization at the triple point. Vaporizing collisions between small bodies are an important component of terrestrial planet formation.
This project plan gives a high-level description of the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Disposition (SFWD) campaign in situ borehole heater test project planned for the Waste Isolation Pilot Plant (WIPP) site, titled the Brine Availability Test in Salt (BATS). BATS is the first stage in a planned sequence of tests to bolster the technical basis for disposal of heat generating waste in salt. This plan provides an overview of the schedule and responsibilities of the parties involved. This project is a collaborative effort by Sandia, Los Alamos, and Lawrence Berkeley National Laboratories to execute a series of small-diameter borehole heater tests in salt for the DOE-NE SFWD campaign. Design of a heater test in salt at WIPP has evolved over several years. The experiment has begun in January 2020 and the first phase will continue for several months with the possibility for follow on testing. BATS comprises a suite of modular tests, which consist of a group of adjacent horizontal boreholes in the wall of drifts at WIPP. Each test is centered around a packer-isolated heated borehole (12.2 cm [4.81 diameter) containing equipment for water-vapor collection and borehole closure monitoring, surrounded by smaller-diameter (up to 5.3 cm [2.11 diameter) satellite observation boreholes. Observation boreholes contain grouted-in temperature sensors, electrical resistivity tomography (ERT) sensors, and fiber optics; packer-isolated tracer release and sampling intervals; and acoustic emission (AE) piezoelectric sensors. A larger-diameter (12.2 cm [4.81) satellite borehole includes sorel and salt cement plugs, as part of an engineered barrier sealing test. The first two tests, to be implemented in parallel, are heated (target borehole wall temperature of 120 °C) and unheated, with similar arrays of observation borehole monitoring changes. Follow-on tests will be designed using information gathered from the first two tests, and may be conducted at other borehole wall temperatures, use multiple observation boreholes, and may include different measurement types and test designs. This 2020 update of the original 2018 project plan satisfies DOE-NE Spent Fuel and Waste Science and Technology (SFWST) milestone M3SF-205N010303034, as part of the SNL "Salt Disposal R&D" work package.
This report describes research and development (R&D) activities conducted during fiscal year 2019 (FY19) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Eneregy (DOE). The R&D activities focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Genreric Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc.The FY19 EBS activities involved not only modeling and analysis work, but experimental work as well. The report documents the FY19 progress made in seven different research areas as follows: (1) thermal analysis for the disposal of dual purpose canisters (DPCs) in sedimentary host rock using the semianalytical method, (2) tetravalent uranium solubility and speciation, (3) modeling of high temperature, thermal-hydrologic-mechanical-chemical (THMC) coupled processes, (4) integration of coupled thermalhydrologic- chemical (THC) model with GDSA using a Reduced-Order Model, (5) studying chemical controls on montmorillonite structure and swelling pressure, (6) transmission x-ray microscope for in-situ nanotomography of bentonite and shale, and (7) in-situ electrochemical testing of uranium dioxide under anoxic conditions. The R&D team consisted of subject matter experts from Sandia National Laboratories, Lawrence Berkeley National Laboratory (LBNL), Los Alamos National Laboratory (LANL), Pacific Northwest National Laboratory (PNNL), the Bureau de Recherches Géologiques et Minières (BRGM), the University of California Berkeley, and Mississippi State University. In addition, the EBS R&D work leverages international collaborations to ensure that the DOE program is active and abreast of the latest advances in nuclear waste disposal. For example, the FY19 work on modeling coupled THMC processes at high temperatures relied on the bentonite properties from the Full-scale Engineered Barrier EXperiment (FEBEX) Field Test conducted at the Grimsel Test Site in Switzerland. Overall, significant progress has been made in FY19 towards developing the modeling tools and experimental capabilities needed to investigate the performance of EBS materials and the associated interactions in the drift and the surrounding near-field environment under a variety of conditions including high temperature regimes.
Security models for the bioeconomy have largely been developed on risk profiles borrowed from the financial industry and, in some cases, industrial control systems. The bioeconomy, however, has unique characteristics that require domain-specific knowledge of the risks to environmental, health and economic impact from directly targeted attacks. Key questions need to be assessed and answered for any facility engaged in bioproduction. These include: how much technical information must the attacker possess in order to target a specific process or facility? Which processes cause the largest economic impact if disrupted? Can an attacker lead companies down the wrong path of research, leading to irrecoverable losses of time, resources and capital? Can attackers disrupt venture capital strategies and affect financial returns? What are the resources required to attack key workflows, and subsequently what is the cost of defense? What strategies are effective against such attackers, and what are their cost? Can government provide an active role in assurance of material, process and results for critical bioeconomic infrastructure?
A recently completed LDRD project has provided a new facility at Sandia for testing effects of energetic neutrons on electronic components. 14 MeV neutrons are produced with a deuterium ion beam onto a thin-film tritide target. The goal of the project was to increase the neutron fluence to levels needed for radiation effects testing and qualification. This goal was achieved through two technical advances. First, a new multi-layer target concept was developed to reduce the rate of tritium loss from the target by isotope exchange, thereby reducing tritium usage and increasing target lifetime. The second advance was the construction of a new test chamber designed to maximize neutron flux at test locations. Together, these increased the available neutron fluence by several orders of magnitude. This new capability is being used in tests for Sandia nuclear weapon programs, evaluation of commercial parts such as highly-scaled CMOS SRAM lCs, and tests of new devices under development at Sandia such as lll-V HBTs, gallium nitride high-voltage diodes, and for fundamental studies of physical mechanisms of device failure.
Assess the evolving integration potential of light-duty (LDV) and heavy-duty vehicle (HDV) technologies, fuels, and infrastructure and their contributions to lowering emissions and petroleum consumption. Leverage existing LDV and build HDV ParaChoice capability to conduct parametric analyses that explore the trade-space for key factors that influence consumer choice and technology, fuel and infrastructure development. ParaChoice provides the unique capability to examine tipping points and tradeoffs, and can help quantify the effects of and mitigate uncertainty introduced by data sources and assumptions.
In the case for each of the tasks, implementation started in a private topic branch. That branch was later submitted as a merge request where the code was run through regression tests across multiple test platforms. The merge requests were also subjected to human reviewers for approval. After necessary modifications were made, the code was merged to VTK-m's master branch. Subsequently, documentation was written for the VTK-m User's Guide.
The goal of the DOE OE Energy Storage System Safety Roadmap is to foster confidence in the safety and reliability of energy storage systems. There are three interrelated objectives to support the realization of that goal: research, codes and standards (C/S) and communication/coordination. The objective focused on C/S is "To apply research and development to support efforts that refocused on ensuring that codes and standards are available to enable the safe implementation of energy storage systems in a comprehensive, non-discriminatory and science-based manner."
Sandia and SLAC staff evaluated Princeton Plasma Physics Laboratory (PPPL) engineering procedures. The report documents, team members, interviewees and documents consulted, lines of inquiry, strengths and best practices, weaknesses and risks, opportunities for improvement, and general comments.
Transition metal ions provide a rich set of optically active defect spins in wide bandgap semiconductors. Chromium (Cr4+) in silicon-carbide (SiC) produces a spin-1 ground state with a narrow, spectrally isolated, spin-selective, near-telecom optical interface. However, previous studies were hindered by material quality resulting in limited coherent control. In this work, we implant Cr into commercial 4H-SiC and show optimal defect activation after annealing above 1600 °C. We measure an ensemble optical hole linewidth of 31 MHz, an order of magnitude improvement compared to as-grown samples. An in-depth exploration of optical and spin dynamics reveals efficient spin polarization, coherent control, and readout with high fidelity (79%). We report T1 times greater than 1 s at cryogenic temperatures (15 K) with a T2* = 317 ns and a T2 = 81 μs, where spin dephasing times are currently limited by spin-spin interactions within the defect ensemble. Our results demonstrate the potential of Cr4+ in SiC as an extrinsic, optically active spin qubit.
Organic linkers in metal-organic framework (MOF) materials exhibit differences in hydrogen bonding (H-bonding), which can alter the geometric, electronic, and optical properties of the MOF. Density functional theory (DFT) simulations were performed on a photoluminescent Y-2,5-dihydroxyterephthalic acid (DOBDC) MOF with H-bonding concentrations between 0 and 100%; the H-bonds were located on both bidentate-and monodentate-bound DOBDC linkers. At 0% H-bond concentration in the framework, the lattice parameters contracted, the density increased, and simulated X-ray diffraction patterns shifted. Comparison with published experimental data identified that Y-DOBDC MOF structures must have a degree of H-bond concentration. The concentration of H-bonds in the system shifted the calculated band gap energy from 2.25 eV at 100% to 3.00 eV at 0%. The band gap energies also indicate a distinction of H-bonds formed on bidentate-coordinated linkers compared to those on monodentate linkers. Additionally, when the calculated optical spectra are compared with experimental data, the ligand-to-ligand charge-transfer luminescence in Y-DOBDC MOFs is expected to result from an average of 20-40% H-bonding with at least 50% of the bidentate linkers containing H-bonding. Therefore, the type of H-bonding within the DOBDC linker determines the electronic structure and the optical absorption of the MOF framework structure. Tuning of the H-bonding in rare-earth MOFs provides an opportunity to control the specific optical and adsorption properties of the MOF framework on the basis of reactions between the linker and the environment.
Bacterial tRNA have been found interrupted at various positions in the anticodon loop by group I introns, in four types. The primary bioinformatic tool for group I intron discovery is a covariance model that can identify conserved features in the catalytic core and can sometimes identify the typical uridine residue at the -1 position, preceding the 5-prime splice site, but cannot identify the typical guanidine residue at the omega position, preceding the 3-prime splice site, to achieve precise mapping. One approach to complete the automation of group I intron mapping is to focus instead on the exons, which is enabled by the regularity of tRNAs. We develop a software module, within a larger package (tFind) aimed at mapping bacterial tRNA and tmRNA genes precisely, that expands this list of four known classes of intron-interrupted tRNAs to 21 cases. A new covariance model for these introns is presented. The wobble base pair formed by the -1 uridine is considered a determinant of the 5-prime splice site, yet one reasonably large new type bears a cytidine nucleotide at that position.
All experiments involving the Pecos target chamber in calendar year (CY) 2019 were dedicated to Pre-Heat studies in the context of Magnetized Liner Inertial Fusion (MagLIF). Activities at the target area included actual laser shots but also diagnostic and maintenance, and preparatory work for experiments with cryogenically cooled targets, which are anticipated for CY 2020. Since the Z-Beamlet and Z-Petawatt lasers support multiple campaigns, we can only anticipate laser and laser-operator time for one out of three shot windows in a day, and typically not on every day of the week. For that reason, many of our non-shot activities need to be traded against shots as well.
Fatty alcohols (FOHs) are important feedstocks in the chemical industry to produce detergents, cosmetics, and lubricants. Microbial production of FOHs has become an attractive alternative to production in plants and animals due to growing energy demands and environmental concerns. However, inhibition of cell growth caused by intracellular FOH accumulation is one major issue that limits FOH titers in microbial hosts. In addition, identification of FOH-specific exporters remains a challenge and previous studies towards this end are limited. To alleviate the toxicity issue, we exploited nonionic surfactants to promote the export of FOHs in Rhodosporidium toruloides, an oleaginous yeast that is considered an attractive next-generation host for the production of fatty acid-derived chemicals. Our results showed FOH export efficiency was dramatically improved and the growth inhibition was alleviated in the presence of small amounts of tergitol and other surfactants. As a result, FOH titers increase by 4.3-fold at bench scale to 352.6 mg/L. With further process optimization in a 2-L bioreactor, the titer was further increased to 1.6 g/L. The method we show here can potentially be applied to other microbial hosts and may facilitate the commercialization of microbial FOH production.
Terry steam turbines are employed in the safety systems of many nuclear Boiling Water Reactors to drive pumps and provide cooling water to the nuclear reactor core. While the turbine efficiency is low, the more important feature is high reliability under off-normal conditions. An important aspect of reliability is the ability to function with two-phase steam-water injection into the turbine, as most likely occurred in the Fukushima Dai-ichi nuclear accidents. This study investigates the characteristics of a Terry turbine during air-water injection with gas mass fractions ranging from 1 (dry gas) to 0.05 (wet gas), to better understand the Terry turbine's true operational capabilities and provide justification for extended Terry turbine use for reactor safety. Other parameters investigated are the inlet pressure, the exhaust backpressure and the turbine's rotational speed. The turbine performance is presented in terms of dynamometer loading and pump performance change as functions of the gas mass fraction.
This chapter first describes the traditional view of “context” in systems engineering and identifies challenges to this view related to “the Fourth Industrial Revolution”. It then explores gaps in traditional views, introduces nontraditional approaches to context for systems, and provides more detail on the “context of use” concept for advanced systems engineering. In response to technological evolution(s), advanced systems engineering should seek to more clearly and comprehensively describe operating environments - to include accounting for contextual descriptions consisting of the interrelated human behavior, social, and organizational factors that impact system performance and success. Three academic literatures - systems theory, organization science, and engineering systems - offer insights to better understand and incorporate context into advanced systems engineering. To further make the case for including the context of use in advanced systems engineering, the chapter explores improving systems engineering approaches for security at high consequence facilities.
Simultaneous control of nanoscale surface morphology and composition remains a challenge in preparing bimetallic catalysts, particularly at the large scale required for industrial application and with high-surface-area substrates. Atomic layer electroless deposition (ALED) is a scalable approach to prepare surface-modified metal powders in which elements more noble than the surface hydrides of the substrate metal are deposited layer-by-layer in a surface-limited fashion. Herein we demonstrate that high-surface-area Pt powder is a viable substrate for controlled deposition of Pd adlayers using this technique, with the potential for large-scale preparation, for use in electrocatalytic and catalytic applications such as fuel cells and functionalization of petrochemical feedstocks. Two different growth mechanisms have been proposed based on bulk and surface Pd atomic fractions obtained from atomic absorption spectroscopy and X-ray photoelectron spectroscopy, respectively. Further, spectral simulations were performed to strengthen the proposed growth mechanisms, favoring conformal growth in initial deposition followed by island formation in subsequent cycles. Observation of multiple pathways suggests a means of controlling adlayer surface morphology of ALED materials, in which an initial cycle of deposition sets the fractional coverage and subsequent cycles tune adlayer thickness.
Determine what the baseline radiological concentrations are at SNL/NM for water discharge to the sanitary sewer system and provide limits for other analyzed radionuclides.
This report documents the work performed in characterizing the vacuum conductance of porous structures made by conventional sintering (purchased as commercially available products from Mott Corporation ranging from 20 to 100 media grade) and additive manufacturing (powder bed fusion). The additively manufactured structures described in this report were originally intended to be the first iteration of several in an effort to produce desirable conductance characteristics. While resources were not available to link the experimental results to a modeling effort to better understand why certain characteristics were observed, the author hopes that this report may provide a useful set of data for future use, especially as sintered porous structures are not uncommonly procured from Mott corporation for research and development purposes. For that reason, all of the raw data is tabulated in the Appendix: it is possible to reproduce all of the figures shown in this report independently. For a quick order-of magnitude scan of the conductances, refer to Figures 5, 6, 8 and 9.
With the dawn of the exascale era, computer scientists and engineers are faced with tremendous challenges across all facets of the HPC system - scalability, performance, reliability, and power consumption. In particular, the power-performance benefit from one processor generation to the next is seeing ever-diminishing returns and will require fundamental changes in the way we approach computation. In fact, it is likely that different applications will require different types of accelerators in order to meet power, performance, and reliability requirements at scale. One potential type of accelerator, a dataflow architecture, diverges from the traditional sequentially executed instruction model into one that reflects the inherent instruction-level parallelism in a program. This work presents the initial steps toward a tool that can extract the control-dataflow graph from an application.
The performance of Nuclear Deterrence (ND) systems could be improved by adopting advanced commercial CMOS technologies. However, the radiation hardness of these ND systems must be assured. This project quantified dose rate upset thresholds, allowing us to evaluate whether advanced commercial technologies can be used in rad hard applications. This project also evaluated the susceptibility of advanced commercial technologies to neutron displacement damage and single event effects, and developed hardness assurance methods. As a result of this work, Sandia is now collaborating more closely with DoD agencies and their contractors to understand and improve the radiation hardness of advanced commercial technologies. Sandia Capabilities (SPHINX, lon Beam Lab, FPGA test capabilities) were developed and are being utilized by other programs. Staff were trained to do radiation survivability testing and developed proposals that were funded. Internal and external collaborations and partnerships were developed, including with Georgia Tech. This project received the Best Paper award at the 2019 Hardened Electronics And Radiation Technology (HEART) Conference.
The high-energy upgrade of the Z-Petawatt (ZPW) laser system involved a shot hiatus starting in October 2018 to allow the installation of full-aperture optics in the chain, including in the main amplifier section as well as in the branches associated with the ZPW compressor and co-injection (the narrowband nanosecond operating mode that is co-bored with Z-Beamlet). Laser shots resumed May 1st, with most shots since being in the ns mode to reduce complexity (while conditioning optics and calibrating diagnostics) as well as to support experimental requests. To facilitate experimental needs while a final full-sized co-injection optic was prepared, a sub-apertured beam has been used for experimental applications and a full-apertured beam has been used to validate the system. In the latter case, the beam is directed to a large calorimeter placed before branching occurs to the ZPW compressor or co-injection. As of September 2019, full-aperture co-injection shots had ramped up to 960J in 2.3ns, with many cross-calibrations being used to map out a relative humidity dependency in the leaky diagnostic mirror. Several such cross-calibration shots have been done using a chirped pulse seed as well, verifying that the diagnostic calibrations are valid in both broadband and narrow-band modes. Summarizing, since May 1st of 2019 when shots resumed, 93 full system shots were performed with the ZPW laser. This included 2 co-injection shots with ZBL into Z (in support of MagLIF), 26 Pecos gas cell shots (in support of MagLIF), 11 Conchas gas cell shots (in support of UXI applications development), 6 Jemez shots (in support of backlighter development), and 48 calorimeter shots (in support of validating the system performance). Of these 93 shots, 45 shots were in conjunction with ZBL.
Fast charging of lithium ion batteries is a critical enabler for mass EV adoption. Sandia National Laboratories (SNL) has been working with the University of Michigan to develop graphite anodes with novel 3D structures that facilitate faster charging while avoiding lithium plating, a main danger of unaided fast charging. SNL is using its unique high precision cycling capability, developed through ARPA-E funding, to characterize the ability of improved anodes to withstand fast charge and resist lithium plating, and the danger of lithium plating in present-day batteries.
Barlos, Fotis; Skinner, Anna; Peeke, Richard; Mohan, Robert; Kelic, Andjelka; Beyeler, Walter; Homan, Rossitza; Starz, James; Hoffman, Mark; Hofmann, Martin O.; Guo, Katherine; Van Brackle, David; Pioch, Nicholas; Alonso, Rafael; Brown, Kerry; Miller, Scott; Diehl, Michael; Ma, Laura P.; Basu, Prithwish; Wright, William; Shellman, Stephen M.; Hazard, Chris; Brown, Matthew; Fang, Fei; Shen, Weiran; Geib, Christopher
An emergent type of geopolitical warfare in recent years has been coined "Gray Zone competition," "competition short of armed conflict", or simply "competition," because it sits in a nebulous area between peace and conventional conflict. As such, the Gray Zone (GZ) "is a conceptual space between peace and war, occurring when actors use instruments of power to achieve political-security objectives ... [that] ... fall below the level of large-scale direct military conflict". It's not openly declared or defined; it's slower and is prosecuted more subtly using social, psychological, religious, information, cyber and other means to achieve physical or cognitive objectives, with or without violence. The lack of clarity of intent in competition activity makes it challenging to detect, characterize, and counter an enemy fighting this way.
Graph neural networks (GNNs) are an emerging model for learning graph embeddings and making predictions on graph structured data. However, robustness of graph neural networks is not yet well-understood. In this work, we focus on node structural identity predictions, where a representative GNN model is able to achieve near-perfect accuracy. We also show that the same GNN model is not robust to addition of structural noise, through a controlled dataset and set of experiments. Finally, we show that under the right conditions, graph-augmented training is capable of significantly improving robustness to structural noise.
The stakes have been raised with the onslaught of Information and Communications technological disruptive change coupled with wholesale integration of Information Technologies (IT) with Operational Technologies (OT). Industrial Control Systems (ICS) and Supervisory Control and Data Acquisition (SCADA) systems are now virtualized, automated, computerized and connected. Examples include Smart Grid and Smart Metering in the Energy sector; Automated Self Driving Vehicles in the Transportation Systems sector; ICS/SCADA modernization in the Water and Wastewater Systems sector; and Telemedicine in the Healthcare and Public Health sector. Vast broadband expansion brings new connectivity to rural and tribal communities with impacts to business, education, anchor institutions, and first responders.
The increasingly large payloads of Unmanned Aircraft Systems (UASs) are exponentially increasing the threat to the nuclear enterprise. Current mitigation using RF interference is effective, but not feasible for fully autonomous systems and is prohibited in many areas. A new approach to UAS threat mitigation is needed that does not create radio interference but is effective against any type of vehicle. At the present time there is no commercial counter-UAS system that directly assaults the mems gyros and accelerometers in the Inertial Measurement Unit on the aircraft. But lab testing has revealed resonances in some IMUs that make them susceptible to moderate amplitude acoustic monotones. Sandia's energetic materials facility has enabled a quick and thorough exploration of UAS vulnerability to directed acoustic energy by using intense acoustic impulses to destabilize or down a UAS. We have: 1) detonated/deflagrated explosive charges of various sizes; 2) accurately measured impulse pressure and pulse duration; 3) determined what magnitude of acoustic insult to the IMU disrupts flight and for how long and; 4) determined if the air blast/shock wave on aircraft/propellers disrupts flight.
Sandia's work in micro and nanosystems/microelectronics is critical for DOE and other national security missions. Collocated REtD and production Is foundational to current and future capabilities as are partnerships, requiring continuous engagement in the following areas: Advanced REtD that deepens our scientific understanding of technologies to meet national needs and address emerging threats; External partnerships that strengthen the US semiconductor science and industrial base and leverage US and global supply chains; and Trusted and low-volume production for NNSA and synergistic national security missions that can't, won't, or shouldn't be supplied by industry. Microsystems are the primary enabler for agile, responsive weapon development and advanced functionality in the non-nuclear domain, and Sandia's MESA Complex is the only remaining US foundry with the proven ability to deliver the custom trusted and strategically radiation-hardened (TSRH) microelectronic components required to sustain the nation's nuclear deterrent. Four key strategies will advance MESA's technological leadership against future mission challenges.
This report characterizes Mach number effects in the hypersonic range from Mach 3 to Mach 7 for an arbitrary, computationally generated fragment. Force and moment coefficients in the x, y, z direction were successfully obtained and resulting trends were compared with theoretical expectations. Three orientations were tested in the form of unit vectors: [1,0,0], [0,1,0], and [0,0,1]; these represent the i, j, k principle axes. The final results ultimately showcased a trend very similar to that shown in previous literature — drag coefficients for a given orientation decreased with increasing Mach number by a very small amount in the hypersonic range. It was therefore concluded that the aerodynamic quantities used to obtain the fragment trajectory at low hypersonic Mach numbers (such as Mach 3) can be used to characterize the aerodynamic qualities at higher hypersonic Mach numbers (Mach 5, 7) with a reasonably small margin. A sensitivity study was also conducted which ultimately defined drag effects as a result of changing grid spacing and normal extrusion parameters for the unstructured tetrahedral mesh used in the simulation. Further work can be done in the form of experimental validation, specifically with regards to wind tunnels and ballistic range testing. Error can be reduced by testing more mesh variabilities and capturing a larger amount of fragment orientations. The results for force and moment characteristics dependent on meshing parameters and high Mach numbers are satisfactory and consistent with expected trends.
The Integrated Military Systems (IMS) Center at Sandia National Laboratories is a technical systems integrator of high-risk, high-impact systems. It provides cutting-edge technologies, subsystems, solutions, and products that provide a competitive advantage at every point of engagement. Building on more than six decades of experience, IMS seeks to pursue and provide solutions to the most challenging national security issues. We leverage technical capabilities from the Laboratories' core mission and innovatively apply technology to defense problems. IMS' diverse capabilities produce revolutionary technologies in the areas such as precision strike and defense, autonomy, NG&C, directed energy, signal processing/ATR, lethality and threat assessment, and systems engineering. Our talented engineers and staff take a problem from initial design all the way through testing and fielding.
Abuse tests are designed to determine the safe operating limits of HEV\PHEV energy storage devices. Testing is intended to achieve certain worst-case scenarios to yield quantitative data on cell\module\pack response, allowing for failure mode determination and guiding developers toward improved materials and designs. Standard abuse tests with defined start and end conditions are performed on all devices to provide comparison between technologies. New tests and protocols are developed and evaluated to more closely simulate real-world failure conditions. While robust mechanical models for vehicles and vehicle components exist, there is a gap for mechanical modeling of EV batteries. The challenge with developing a mechanical model for a battery is the heterogeneous nature of the materials and components (polymers, metals, metal oxides, liquids). This year saw the stand up of a new drop tower tester capable of providing dynamic mechanical test results.
We propose a bilevel optimization approach for the determination of parameters in nonlocal image denoising. We consider both spatial weights in front of the fidelity term, as well as weights within the kernel of the nonlocal operator. In both cases we investigate the differentiability of the solution operator in function spaces and derive a first order optimality system that characterizes local minima. For the numerical solution of the problems, we propose a second-order optimization algorithm in combination with a finite element discretization of the nonlocal denoising models and a computational strategy for the solution of the resulting dense linear systems. Several experiments are run in order to show the suitability of our approach.
With current lithium ion batteries optimized for performance under relatively low charge rate conditions, implementation of XFC has been hindered by drawbacks including Li plating, kinetic polarization, and heat dissipation. This project will utilize model-informed design of 3-D hierarchical electrodes to tune key XFCrelated variables like 1) bulk porosity/tortuosity 2) vertical pore diameter, spacing, and lattice 3) crystallographic orientation of graphite particles relative to exposed surfaces 4) interfacial chemistry of the graphite surfaces through "artificial sEr' formation using ALD 5) current collector surface roughness (aspect ratio, roughness factor, etc.). A key aspect of implementing novel electrodes is characterizing them in relevant settings. For this project, ultimately led out of University of Michigan by Neil Dasgupta, that includes both coin cell and 2+ Ah pouch cell testing, as well as comparison testing against baselines. Sandia National Labs will be conducting detailed cell characterization on iterative versions/improvements of the model-based hierarchical electrodes, as well as COTS cells for baseline comparisons. Key metrics include performance under fast charge conditions, as well as the absence or degree of lithium plating. Sandia will use their unique high precision cycling and rapid EIS capabilities to accurately characterize performance and any lithium plating during 6C charging and beyond, coupling electrochemical observations with cell teardown. Sandia will also design custom fixturing to cool cells during rapid charge, to decouple any kinetic effects brought about by cell heating and allow comparisons between different cells and charge rates. Using these techniques, Sandia will assess HOH electrodes from the University of Michigan, as well as aiding in iterative model and electrode design.
Neuromorphic computing is known for its integration of algorithms and hardware elements that are inspired by the brain. Conventionally, this nontraditional method of computing is used for many neural or learning inspired applications. Unfortunately, this has resulted in the field of neuromorphic computing being relatively narrow in scope. In this paper we discuss two research areas actively trying to widen the impact of neuromorphic systems. The first is Fugu, a high-level programming interface designed to bridge the gap between general computer scientists and those who specialize in neuromorphic areas. The second aims to map classical scientific computing problems onto these frameworks through the example of random walks. This elucidates a class of scientific applications that are conducive to neuromorphic algorithms.
This co-op was performed at Sandia National Laboratories (SNL) within the Weapon Surety Engineering 11 Department (9435), which is a customer-focused organization that provides engineering and system analysis expertise to ensure a safe deterrent throughout the entire weapon lifecycle. Engineers in 9435 develop a deep knowledge in both component and system function to support rigorous evaluation and mitigation of potential failure modes. They operate across multiple engineering disciplines (mechanical, electrical, computer, materials and systems engineering) with opportunities to develop crossdisciplinary capability. The team supports a variety of stockpile systems and is part of the product teams during development, design and production, as well as advanced concept and exploratory efforts. 9435 engineers use their extensive component and system knowledge, as well as experimental data and modeling, to ensure products function as intended with assured safety. Members of the department learn the function of multiple components, how those components integrate within the system and apply surety principles and processes to ensure a safe, secure, reliable, and effective U.S. nuclear deterrent.
The purpose of this work is to determine MTD effectiveness and cost for protecting non-IP C2 networks, experiment on networks used in satellite systems, and randomize features (e.g. device address) to prevent adversary from conducting reconnaissance.
The National Technology & Engineering Solutions of Sandia, LLC (NTESS) Prime Contract department manages an organizational conflicts of interest (OCI) program that is designed to provide NTESS employees an accessible means to assist in identifying potential OCI risk and providing development support services to avoid, neutralize, or mitigate the risk. The OCI program is responsible for maintaining a regulatory infrastructure, such as corporate OCI policy and training content, to comply with NTESS's obligations set forth in prime contract clauses H-4: ORGANIZATIONAL CONFLICT OF INTEREST (OCI) — SPECIAL PROVISION and Section I Department of Energy Acquisition Regulation 952.209-72, Organizational Conflicts of Interest, Alternate I, paragraph (c)(1), Disclosure After Award. Risk identification is a key OCI program priority. Providing employees with a straight-forward, accessible and navigable mechanism to identify potential OCI risk will assist the OCI program in meeting this priority.
Lithographic quantum dots (QDs) are highly controllable few-level quantum systems created in semiconductor nanoelectronic devices, with a variety of scientific applications. These include technologically-driven applications like quantum computing and more fundamental applications in which they serve as a platform for exploring basic many-body physics. This document is a brief summary of my Ph.D. research so far and the directions with which I intend to continue it. Highlights include theoretical efforts to understand and design qubits in germanium hole QDs, as well as explorations of the possibility of using QDs coupled to nearby baths for analog simulation of quantum impurity models.
Molybdenum-99 (Moly 99) is a critical raw material for Technetium (Tc) 99m, a radioactive isotope most widely used in nuclear medicine procedures. Moly 99 has a short half-life of about six hours, which means it cannot be stockpiled. When Moly 99 decays, it turns into Tc 99m, which has a half-life of 214,000 years. The photon energy emitted from the decay of Moly 99 is used in a variety of nuclear imaging technologies such as gamma cameras. Furthermore, radiopharmaceutical manufacturers use the photon energy emitted from the decay of Moly 99 to produce generators for hospitals, clinics, and radiopharmacies. Once Moly 99 decays to Tc 99m, it is used to make individual patient doses for a variety of diagnostic imaging procedures. The Moly 99 Reactor Design is a conceptual blueprint for Moly 99 production that does not use weapon-grade uranium; instead, the reactor has a target core of low-enriched uranium. The reactor design is small, reaching a foot-and-a-half in height and diameter and consumes less than two megawatts of power. The reactor's only purpose is for medical isotope production and with every fission, Moly 99 is produced.
Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.
The mechanical properties of additively manufactured metals tend to show high variability, due largely to the stochastic nature of defect formation during the printing process. This study seeks to understand how automated high throughput testing can be utilized to understand the variable nature of additively manufactured metals at different print conditions, and to allow for statistically meaningful analysis. This is demonstrated by analyzing how different processing parameters, including laser power, scan velocity, and scan pattern, influence the tensile behavior of additively manufactured stainless steel 316L utilizing a newly developed automated test methodology. Microstructural characterization through computed tomography and electron backscatter diffraction is used to understand some of the observed trends in mechanical behavior. Specifically, grain size and morphology are shown to depend on processing parameters and influence the observed mechanical behavior. In the current study, laser-powder bed fusion, also known as selective laser melting or direct metal laser sintering, is shown to produce 316L over a wide processing range without substantial detrimental effect on the tensile properties. Ultimate tensile strengths above 600 MPa, which are greater than that for typical wrought annealed 316L with similar grain sizes, and elongations to failure greater than 40% were observed. It is demonstrated that this process has little sensitivity to minor intentional or unintentional variations in laser velocity and power.
Silicon photonics is a platform that enables densely integrated photonic components and systems and integration with electronic circuits. Depletion mode modulators designed on this platform suffer from a fundamental frequency response limit due to the mobility of carriers in silicon. Lithium niobate-based modulators have demonstrated high performance, but the material is difficult to process and cannot be easily integrated with other photonic components and electronics. In this manuscript, we simultaneously take advantage of the benefits of silicon photonics and the Pockels effect in lithium niobate by heterogeneously integrating silicon photonic-integrated circuits with thin-film lithium niobate samples. We demonstrate the most CMOS-compatible thin-film lithium niobate modulator to date, which has electro-optic 3 dB bandwidths of 30.6 GHz and half-wave voltages of 6.7 V×cm. These modulators are fabricated entirely in CMOS facilities, with the exception of the bonding of a thin-film lithium niobate sample post fabrication, and require no etching of lithium niobate.
Lithium-ion battery technology is rapidly being adopted in transportation applications and energy storage industries. Safety concerns, in particular, fire and explosion hazards, are threatening widespread adoption. In some failure events, lithium-ion cells can undergo thermal runaway, which can result in the release of flammable gases that pose fire and explosion hazards for the compartment housing the cells. However, there is little available information characterizing the flammability properties of the gases released after cell thermal runaway. In this paper, analytical and modeling methods to estimate explosion characteristics, such as lower flammability limit, laminar flame speed, and maximum over-pressure are evaluated for use in quantifying the effect of cell chemistry, state-of-charge and other parameters on the overall explosion hazard potential for confined cells.
Rempe, Susan R.; Basdogan, Yasemin; Groenenboom, Mitchell C.; Henderson, Ethan; De, Sandip; Keith, John A.
Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1-, and 2-, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.
Molecular dynamics simulations confirm recent extensional flow experiments showing ring polymer melts exhibit strong extension-rate thickening of the viscosity at Weissenberg numbers Wi « 1. Thickening coincides with the extreme elongation of a minority population of rings that grows with Wi. The large susceptibility of some rings to extend is due to a flow-driven formation of topological links that connect multiple rings into supramolecular chains. Links form spontaneously with a longer delay at lower Wi and are pulled tight and stabilized by the flow. Once linked, these composite objects experience larger drag forces than individual rings, driving their strong elongation. The fraction of linked rings depends non-monotonically on Wi, increasing to a maximum when Wi 1 before rapidly decreasing when the strain rate approaches 1/Te.
Hannah, Walter M.; Jones, Christopher R.; Hillman, Benjamin H.; Norman, Matthew R.; Bader, David C.; Taylor, Mark A.; Leung, Lai-Yung R.; Pritchard, Michael S.; Branson, Mark D.; Lin, Guangxing; Pressel, Kyle G.; Lee, Jungmin M.
Results from the new DOE super-parameterized (SP) Energy Exascale Earth System Model (SP-E3SM) are analyzed and compared to the traditionally parameterized E3SMv1 and previous studies using SP models. SP-E3SM is unique in that it utilizes GPU hardware acceleration, CRM mean-state acceleration, and reduced radiation to dramatically increase the model throughput and allow decadal experiments at 100-km external resolution. It also differs from other SP models by using a spectral element dynamical core on a cubed sphere grid and a finer vertical grid with a higher model top. Despite these differences, SP-E3SM generally reproduces the behavior of other super-parameterized models. Tropical wave variability is improved relative to E3SM, including the emergence of a Madden-Julian Oscillation and a realistic slowdown of Moist Kelvin Waves. However, the distribution of precipitation exhibits an unrealistically large variance, and while the timing of diurnal rainfall shows modest improvements the signal is not as coherent as observations. A notable grid imprinting bias is identified in the precipitation field and attributed to a unique feedback associated with the interactions between explicit convection and the spectral element grid structure. Spurious zonal mean column water tendencies due to grid imprinting are quantified – while negligible for the conventionally parameterized E3SM, they become large with super-parameterization, approaching 10% of the physical tendencies. The implication is that finding a remedy to grid imprinting will become especially important as spectral element dynamical cores begin to be combined with explicitly resolved convection.
An improved meshing scheme allows simulations to proceed further than before, though convergence difficulty is still encountered upon full opening. Changing the initial shape of the window does not greatly affect the rate of opening. Changing the initial window stress greatly affects the initial snap-back upon opening, but the long-term trendline as full opening nears is less strongly affected. At 60 psi pressure, the opening time is approximately 4-4.5 μs for the configurations tested in this study.
It is well established that the variability in mechanical response and ultimate failure of additively manufactured metals correlates to uncertainties introduced in the build process, among which include internal void structure and residual stresses. Here, we quantify the aforementioned variabilities in 316L stainless steels by conducting simulations in Sierra/SM of the specimens/geometries used in Sandia's third fracture challenge (SFC3). We leverage the simulations and experimental work presented in 6 to construct a statistical representation of the internal void structure of the tension specimen used for material parameter calibration as well as the "challenge" geometry. Voided mesh samples of both specimens are generated given a set of statistical variables, and the physics simulations are conducted for multiple sets of realization to determine the effects of void structure on variability in the fracture paths and displacement-to-failure. Lastly, a series of simulations are presented which highlight the effect of the powder bed fusion additive manufacturing process on the formation of residual stresses in the as-built geometries.
Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions (SOAP), the Spectral Neighbor Analysis Potential (SNAP) bispectrum components, and moment tensors --- using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model, and consequently computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.
Energy storage systems with Li-ion batteries are increasingly deployed to maintain a robust and resilient grid and facilitate the integration of renewable energy resources. However, appropriate selection of cells for different applications is difficult due to limited public data comparing the most commonly used off-the-shelf Li-ion chemistries under the same operating conditions. This article details a multi-year cycling study of commercial LiFePO4 (LFP), LiNixCoyAl1-x-yO2 (NCA), and LiNixMnyCo1-x-yO2 (NMC) cells, varying the discharge rate, depth of discharge (DOD), and environment temperature. The capacity and discharge energy retention, as well as the round-trip efficiency, were compared. Even when operated within manufacturer specifications, the range of cycling conditions had a profound effect on cell degradation, with time to reach 80% capacity varying by thousands of hours and cycle counts among cells of each chemistry. The degradation of cells in this study was compared to that of similar cells in previous studies to identify universal trends and to provide a standard deviation for performance. All cycling files have been made publicly available at batteryarchive.org, a recently developed repository for visualization and comparison of battery data, to facilitate future experimental and modeling efforts.
The heat generated during a single cell failure within a high energy battery system can force adjacent cells into thermal runaway, creating a cascading propagation effect through the entire system. This work examines the response of modules of stacked pouch cells after thermal runaway is induced in a single cell. The prevention of cascading propagation is explored on cells with reduced states of charge and stacks with metal plates between cells. Reduced states of charge and metal plates both reduce the energy stored relative to the heat capacity, and the results show how cascading propagation may be slowed and mitigated as this varies. These propagation limits are correlated with the stored energy density. Results show significant delays between thermal runaway in adjacent cells, which are analyzed to determine intercell contact resistances and to assess how much heat energy is transmitted to cells before they undergo thermal runaway. A propagating failure of even a small pack may stretch over several minutes including delays as each cell is heated to the point of thermal runaway. This delay is described with two new parameters in the form of gap-crossing and cell-crossing time to grade the propensity of propagation from cell to cell.
The mechanical degradation of all-solid-state Li-ion batteries (ASSLBs) is expected to be more severe than that in traditional Li-ion batteries with liquid electrolytes due to the additional mechanical constraints imposed by the solid electrolyte on the deformation of electrodes. Cracks and fractures could occur both inside the solid electrolyte (SE) and at the SE/electrode interfconce. A coupled electrochemical-mechanical model was developed and solved by the Finite Element Method (FEM) to evaluate the stress development in ASSLBs. Two sources of volume change were considered, namely the expansion/shrinkage of electrodes due to lithium concentration change and the interphase formation at the SE/electrode interface due to the decomposition of SEs. The most plausible solid electrolyte decomposition reactions and their associated volume change were predicted by density functional theory (DFT) calculations. It was found that the stress associated with a volume change due to solid electrolyte decomposition can be much more significant than that of electrode volumetric changes associated with Li insertion/extraction. This model can be used to design 3D ASSLB architectures to minimize their internal stress generation.
On 12/2/2019 at approximately 1747 hours, a vehicle headed east along East Avenue (LLNL property) made contact with a Sandia subcontractor on foot headed north through the marked crosswalk. The Sandia subcontractor required hospitalization as a result of the incident.
Li deposition at the graphitic anode is widely reported in literature as one of the leading causes of capacity fade in lithium-ion batteries (LIBs). Previous literature has linked Li deposition resulting from low-temperature ageing to diminished safety characteristics, however no current research has probed the effects of Li deposition on the abuse response of well-characterized cells. Using overtemperature testing, a relationship between increased concentrations of Li deposition and exacerbated abuse response in 1 Ah pouch cells has been established. A novel Li deposition technique is also investigated, where cells with n:p < 1 (anode-limiting) have been cycled at a high rate to exploit Li+ diffusion limitations at the anode. Scanning Electron Microscopy of harvested anodes indicates substantial Li deposition in low n:p cells after 20 cycles, with intricate networks of Li(s) deposits which hinder Li+ intercalation/deintercalation. Peak broadening and decreased amplitude of differential capacity plots further validates a loss of lithium inventory to Li+ dissolution, and Powder X-ray Diffraction indicates Li+ intercalation with staging in anode interstitial sites as the extent of Li deposition increases. A cradle-to-grave approach is leveraged on cell fabrication and testing to eliminate uncertainty involving the effects of cell additives on Li deposition and other degradation mechanisms.
The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). The dry cask simulator (DCS) was previously built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the validity of the modeling used to determine cladding temperatures in modern vertical dry casks. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. In this study, a model validation exercise was carried out using the data obtained from dry cask simulator testing in the vertical, aboveground configuration. Five modeling institutions – Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), Centro de Investigaciones Energéticas, MedioAmbientales y Tecnológicas (CIEMAT), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) in collaboration with Universidad Politécnica de Madrid (UPM) – were granted access to the input parameters from SAND2017-13058R, “Materials and Dimensional Reference Handbook for the Boiling Water Reactor Dry Cask Simulator”, and results from the vertical aboveground BWR dry cask simulator tests reported in NUREG/CR-7250, “Thermal-Hydraulic Experiments Using A Dry Cask Simulator”. With this information, each institution was tasked to calculate minimum, average, and maximum fuel axial temperature profiles for the fuel region as well as the axial temperature profiles of the DCS structures. Transverse temperature profiles and air mass flow rates within the dry cask simulator were also calculated. These calculations were done using modeling codes (ANSYS FLUENT, STARCCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, four test cases of the vertical, aboveground dry cask simulator were considered, defined by two independent variables – either 0.5 kW or 5 kW fuel assembly decay heat, and either 100 kPa or 800 kPa internal helium pressure. For the results in this report, each model was assigned a model number. Three of the models used porous media model representations of the fuel, two models used explicit fuel representations, and one model used an explicit subchannel representation of the fuel. Even numbers were assigned to explicit fuel models and odd numbers were assigned to porous media models. The plots provided in Chapter 3 of this report show the axial and transverse temperature profiles obtained from the dry cask simulator experiments in the aboveground configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), average fuel temperatures across six axial levels, transverse temperatures across the PCT locations for the four test cases, and air mass flow rates. The peak cladding temperature is typically the most important target variable for cask performance, and all models capture the PCT within 5% RMS error. Two models show comparable fits to experimental results when considering the combined RMS error of all target variables. Since one uses a porous media representation of the fuel while the other uses an explicit fuel representation, it can be concluded that the porous media fuel representation can achieve modeling calculation results of peak cladding temperatures, average fuel temperatures, transverse temperatures, and air mass flow rates that are comparable to explicit fuel representation modeling results.
Nuclear pore complexes (NPCs) regulate all cargo traffic across the nuclear envelope. The transport conduit of NPCs is highly enriched in disordered phenylalanine/glycine-rich nucleoporins (FG-Nups), which form a permeability barrier of still elusive and highly debated molecular structure. Here we present a microfluidic device that triggered liquid-to-liquid phase separation of FG-Nups, which yielded droplets that showed typical properties of a liquid state. On the microfluidic chip, droplets were perfused with different transport-competent or -incompetent cargo complexes, and then the permeability barrier properties of the droplets were optically interrogated. We show that the liquid state mimics permeability barrier properties of the physiological nuclear transport pathway in intact NPCs in cells: that is, inert cargoes ranging from small proteins to large capsids were excluded from liquid FG-Nup droplets, but functional import complexes underwent facilitated import into droplets. Collectively, these data provide an experimental model of how NPCs can facilitate fast passage of cargoes across an order of magnitude in cargo size.
Exposure testing was performed on CoCrFeMnNi equiatomic high entropy alloy (HEA) produced via directed energy deposition additive manufacturing in NaNO3-KNO3 (60-40 wt%) molten salt at 500 °C for 50 h to evaluate the corrosion performance and oxide film chemistry of the HEA. Potentiodynamic electrochemical corrosion testing, scanning electron microscopy, focused ion beam milling coupled with energy dispersive spectroscopy, Raman spectroscopy, X-ray diffraction, X-ray photoelectron spectroscopy, and inductively coupled plasma optical emission spectroscopy were used to analyze the corrosion behavior and chemistry of the HEA/nitrate molten salt system. The CoCrFeMnNi HEA exhibited a higher passive current density during potentiodynamic polarization testing than steel alloys SS316L and 4130 and the high-Ni alloy 800 H in identical conditions. The oxide film was primarily composed of a (Mn,Co,Ni)Fe2O4 spinel with a vertical plate-like morphology at the surface. Cr and Ni were found to be totally depleted at the outer surface of the oxide and dissolved in high concentrations in the molten salt. While Cr was expected to dissolve into the molten salt, the high concentration of dissolved Ni has not been observed with traditional alloys, suggesting that Ni is less stable in the spinel when Mn and Co are present.
Rechargeable alkaline batteries containing zinc anodes suffer from redistribution of active material due to the high solubility of ZnO in the electrolyte, limiting achievable capacity and lifetime. Here, we investigate pre-saturating the KOH electrolyte with ZnO as a strategy to mitigate this issue, utilizing rechargeable Ni-Zn cells. In contrast to previous reports featuring this approach, we use more practical limited-electrolyte cells and systematically study ZnO saturation at different levels of zinc depth-of-discharge (DODZn), where the pre-dissolved ZnO is included in the total system capacity. Starting with 32 wt. % KOH, cells tested at 14%, 21%, and 35% DODZn with ZnO-saturated electrolyte exhibit 191%, 235%, and 110% longer cycle life respectively over identically tested cells with ZnO-free electrolyte, with similar energy efficiency and no voltage-related energy losses. Furthermore, anodes cycled in ZnO-saturated electrolyte develop more favorable compact zinc deposits with less overall mass loss. The effect of initial KOH concentration was also studied, with ZnO saturation enhancing cycle life for 32 wt % and 45 wt % KOH but not for 25 wt % KOH, likely due to cell failure by passivation rather than shorting. The simplicity of ZnO addition and its beneficial effect at high zinc utilization make it a promising means to make secondary alkaline zinc batteries more commercially viable.
The U.S. Navy is investing in the development of new technologies that broaden warship capabilities and maintain U.S. naval superiority. Specifically, Naval Sea Systems Command (NAVSEA) is supporting the development of power systems technologies that enable the Navy to realise an all-electric warship. A challenge to fielding an all-electric power system architecture includes minimising the size of energy storage systems (ESS) while maintaining the response times necessary to support potential pulsed loads. This work explores the trade-off between energy storage size requirements (i.e. mass) and performance (i.e. peak power, energy storage, and control bandwidth) in the context of a power system architecture that meets the needs of the U.S. Navy. In this work, the simulated time domain responses of a representative power system were evaluated under different loading conditions and control parameters, and the results were considered in conjunction with sizing constraints of and estimated specific power and energy densities of various storage technologies. The simulation scenarios were based on representative operational vignettes, and a Ragone plot was used to illustrate the intersection of potential energy storage sizing with the energy and power density requirements of the system. Furthermore, the energy storage control bandwidth requirements were evaluated by simulation for different loading scenarios. Two approaches were taken to design an ESS: one based only on time domain power and energy requirements from simulation and another based on bandwidth (specific frequency) limitations of various technologies.
An open question in the metal hydride community is whether there are simple, physics-based design rules that dictate the thermodynamic properties of these materials across the variety of structures and chemistry they can exhibit. While black box machine learning-based algorithms can predict these properties with some success, they do not directly provide the basis on which these predictions are made, therefore complicating the a priori design of novel materials exhibiting a desired property value. In this work we demonstrate how feature importance, as identified by a gradient boosting tree regressor, uncovers the strong dependence of the metal hydride equilibrium H2 pressure on a volume-based descriptor that can be computed from just the elemental composition of the intermetallic alloy. Elucidation of this simple structure-property relationship is valid across a range of compositions, metal substitutions, and structural classes exhibited by intermetallic hydrides. This permits rational targeting of novel intermetallics for high-pressure hydrogen storage (low-stability hydrides) by their descriptor values, and we predict a known intermetallic to form a low-stability hydride (as confirmed by density functional theory calculations) that has not yet been experimentally investigated.
The effects of 14-MeV neutron displacement damage (DD) on waveguide (WG)-integrated germanium-on-silicon p-i-n photodiodes (PDs) for silicon photonics have been investigated up to the fluences of 7.5× 1012 n/cm2 (14 MeV) or 1.4× 1013 n1-MeVeq/cm2(Si). This article includes the measurements of dark current-voltage characteristics across temperature from 150 to 375 K, measurements of PD junction capacitance, spectral response measurements from 1260 to 1360 nm, and frequency-response measurements. The devices are found to be susceptible to DD-induced carrier removal effects; however, they also continue to operate without meaningful impact to performance for the DD dose levels examined. Since the PD test chips include silicon photonic integrated grating couplers and WGs, which carry the optical signal to the PD, some assessment of the impact of DD on these passive devices can also be inferred. This article does not examine the short-term annealing or transient behavior of the DD, and instead, it has only considered the lasting damage that remains after any initial period of room-temperature annealing.
The US Department of Energy Office of Nuclear Energy is conducting a brine availability heater test to characterize the thermal, mechanical, hydrological and chemical response of salt at elevated temperatures. In the heater test, brines will be collected and analyzed for chemical compositions. In order to support the geochemical modeling of chemical evolutions of the brines during the heater test, we are recalibrating and validating the solubility models for the mineral constituents in salt formations up to 100°C, based on the solubility data in multiple component systems as well as simple systems from literature. In this work, we systematically compare the model-predicted values based on the various solubility models related to the constituents of salt formations, with the experimental data. As halite is the dominant constituent in salt formations, we first test the halite solubility model in the Na-Mg-Cl dominated brines. We find the existing halite solubility model systematically over-predict the solubility of halite. We recalibrate the halite model, which can reproduce halite solubilities in Na-Mg-Cl dominated brines well. As gypsum/anhydrite in salt formations controls the sulfate concentrations in associated brines, we test the gypsum solubility model in NaCl solutions up to 5.87 mol•kg-1 from 25°C to 50°C. The testing shows that the current gypsum solubility model reproduces the experimental data well when NaCl concentrations are less than 1 mol•kg-1. However, at NaCl concentrations higher than 1, the model systematically overpredicts the solubility of gypsum. In the Na - Cl - SO4 - CO3 system, the validation tests up to 100°C demonstrate that the model excellently reproduces the experimental data for the solution compositions equilibrated with one single phase such as halite (NaCl) or thenardite (Na2SO4), with deviations equal to, or less than, 1.5 %. The model is much less ideal in reproducing the compositions in equilibrium with the assemblages of halite + thenardite, and of halite + thermonatrite (Na2CO3•H2O), with deviations up to 31 %. The high deviations from the experimental data for the multiple assemblages in this system at elevated temperatures may be attributed to the facts that the database has the Pitzer interaction parameters for Cl - CO3 and SO4 - CO3 only at 25°C. In the Na - Ca - SO4 - HCO3 system, the validation tests also demonstrate that the model reproduces the equilibrium compositions for one single phase such as gypsum better than the assemblages of more than one phase.
Kinetic simulations of plasma phenomena during and after formation of the conductive plasma channel of a nanosecond pulse discharge are analyzed and compared to existing experimental measurements. Particle-in-cell with direct simulation Monte Carlo collisions (PIC-DSMC) modeling is used to analyze a discharge in helium at 200 Torr and 300 K over a 1 cm gap. The analysis focuses on physics that would not be reproduced by fluid models commonly used at this high number density and collisionality, specifically non-local and stochastic phenomena. Similar analysis could be used to improve the predictive capability of lower fidelity or reduced order models. First, the modeling results compare favorably with experimental measurements of electron number density, temperature, and 1D electron energy distribution function at the same conditions. Second, it is shown that the ionization wave propagates in a stochastic, stepwise manner, dependent on rare, random ionization events ahead of the ionization wave when the ionization fraction in front of the ionization wave is very low, analagous to the stochastic branching of streamers in 3D. Third, analysis shows high-energy runaway electrons accelerated in the cathode layer produce electron densities in the negative glow region over an order of magnitude above those in the positive column. Future work to develop reduced order models of these two phenomena would improve the accuracy of fluid plasma models without the cost of PIC-DSMC simulations.
Classical molecular dynamics (MD) simulations were performed to provide a conceptual understanding of the amorphous-crystalline interface for a candidate negative thermal expansion (NTE) material, ZrW2O8. Simulations of pressure-induced amorphization at 300 K indicate that an amorphous phase forms at pressures of 10 GPa and greater, and this phase persists when the pressure is subsequently decreased to 1 bar. However, the crystalline phase is recovered when the slightly distorted 5 GPa phase is relaxed to 1 bar. Simulations were also performed on a two-phase model consisting of the high-pressure amorphous phase in direct contact with the crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the crystalline-amorphous interface. Differences in local atomic structure at the interface are quantified from the simulation trajectories.
Pyroshock events from the actuation of separation devices in satellites and launch vehicles are potentially damaging, very short, high intensity events with high frequency content. The pyroshock damage risk is mitigated somewhat by the fact that the shock intensity is attenuated by the spacecraft structure. The NASA and MIL standards, developed from extensive tests performed in the 1960’s, provide pyroshock attenuation guidelines for various structures common to spacecraft and launch vehicles. In this paper, we present the results from a numerical investigation of pyroshock attenuation in cylindrical shell structures. Pyroshock events were modeled using Sandia National Laboratories’ engineering mechanics simulation codes, specifically Sierra/SD. Upon verifying the numerical simulation results against a NASA-HDBK-7005 curve, various structural features were added and design variables were varied to investigate their effects on pyroshock wave propagation and attenuation. The results showed that current numerical simulation tools, given appropriate tuning parameters, are capable of modeling pyroshock events in a simple cylindrical geometry at a reasonable cost. The numerical simulations showed that the presence of geometric features had greater attenuating effects than previously understood. However, shock attenuation levels were less sensitive to design variables of the structural features than expected.
Porosity in additively manufactured metals can reduce material strength and is generally undesirable. Although studies have shown relationships between process parameters and porosity, monitoring strategies for defect detection and pore formation are still needed. In this paper, instantaneous anomalous conditions are detected in-situ via pyrometry during laser powder bed fusion additive manufacturing and correlated with voids observed using post-build micro-computed tomography. Large two-color pyrometry data sets were used to estimate instantaneous temperatures, melt pool orientations and aspect ratios. Machine learning algorithms were then applied to processed pyrometry data to detect outlier images and conditions. It is shown that melt pool outliers are good predictors of voids observed post-build. With this approach, real time process monitoring can be incorporated into systems to detect defect and void formation. Alternatively, using the methodology presented here, pyrometry data can be post processed for porosity assessment.
A hybrid femtosecond/picosecond CARS instrument probed the Q-branch of molecular hydrogen in the multiphase plume of an aluminized solid propellant burn. A single 50 fs regenerative amplifier pumped an OPA and etalon, providing the Stokes and probe pulses respectively. The spectra were recorded at 1 kHz and fit to synthetic spectra to infer the gas rotational temperature. Recorded spectra required dynamic background corrections due to the intense emission of the propellant plume. Two different days of propellant burns were studied, with the lessons learned from nonresonant background issues with the first test applied to the second. For the second attempt, three burns were examined, with mean temperatures differing only by 30 K with a combined mean of 2574 K.
Two techniques have extended the effective frequency limits of postage-stamp PIV, in which a pulse-burst laser and very small fields of view combine to achieve high repetition rates. An interpolation scheme reduced measurement noise, raising the effective frequency response of previous 400-kHz measurements from about 120 kHz to 200 kHz. The other technique increased the PIV acquisition rate to very nearly MHz rates (990 kHz) by using a faster camera. Charge leaked through the camera shift register at these framing rates but this was shown not to bias the measurements. The increased framing rate provided oversampled data and enabled use of multi-frame correlation algorithms for a lower noise floor, increasing the effective frequency response to 240 kHz where the interrogation window size begins to spatially filter the data. Good agreement between the interpolation technique and the MHz-rate PIV measurements was established. The velocity spectra suggest turbulence power-law scaling in the inertial subrange steeper than the theoretical-5/3 scaling, attributed to an absence of isotropy.
This paper presents a practical methodology for propagating and combining the effects of random variations of several continuous scalar quantities and several random-function quantities affecting the failure pressure of a heated pressurized vessel. The random functions are associated with stress-strain curve test-to-test variability in replicate material strength tests (uniaxial tension tests) on nominally identical material specimens. It is demonstrated how to effectively propagate the curve-to-curve discrete variations and appropriately account for the small sample size of functional data realizations. This is coordinated with the propagation of aleatory variability described by uncertainty distributions for continuous scalar quantities of pressure-vessel wall thickness, weld depth, and thermal-contact factor. Motivated by the high expense of the pressure vessel simulations of heating, pressurization, and failure, a simple dimension-and order-adaptive polynomial response surface approach is used to propagate effects of the random variables and enable uncertainty estimates on the error contributed by using the surrogate model. Linear convolution is used to aggregate the resultant aleatory uncertainty from the parametrically propagated random variables with an appropriately conservative probability distribution of aleatory effects from propagating the multiple stress-strain curves for each material. The response surface constructions, Monte Carlo sampling of them for uncertainty propagation, and linear sensitivity analysis and convolution procedures, are demonstrated with standard EXCEL spreadsheet functions (no special software needed).
Bench-top tests are conducted to characterize Femtosecond Laser Electronic Excitation Tagging (FLEET) in static low pressure (35 mTorr-760 Torr) conditions, and to measure the acoustic disturbance caused by the resulting filament as a function of tagging wavelength and energy. The FLEET line thickness as a function of pressure and delay is described by a simple diffusion model. Initial FLEET measurements in a Mach 8 flow show that gate times of ≥ 1µs can produce visible smearing of the FLEET emission and challenge the traditional Gaussian fitting methods used to find the line center. To minimize flow perturbations and uncertainty of the final line position, several recommendations are offered: using third harmonic FLEET at 267 nm for superior signal levels with lower energy deposition than both 800 nm and 400 nm FLEET, and short camera delays and exposure times to reduce fitting uncertainty. This guidance is implemented in a Mach 8 test condition and results are presented.
We demonstrate simultaneous monitoring of temperature and pressure using a hybrid femtosecond/picosecond pure-rotational CARS technique in a one-dimensional line-imaging configuration. The method employs two detection channels and two 60-ps-duration probe laser beams with independently adjustable time delays from the broadband 35-fs pump/Stokes pulse. Simultaneous temperature and pressure monitoring is demonstrated along the centerline of a canonical underexpanded compressible air jet flow emanating from a choked, sonic nozzle. Temperature is measured almost independently of pressure by analyzing CARS spectra obtained with a probe pulse near zero time delay for nearly collision-free acquisition. Pressure is obtained from spectra acquired with long probe time delays to sample the impact of gas-phase collisions. The CARS measurements were obtained in both time-averaged and single-laser-shot mode with 67 µm spatial resolution along the jet axis along a nominally 6-mm line. The measurements span a temperature and pressure range of T = 70-300 K and P = 0.05-1.2 atm.
Multi-Input-Multi-Output (MIMO) vibration testing is considered more representative of the true loading environment (flight or wind induced vibration) where the inputs are not through a single point. The derivation of N inputs for testing typically involves matching the response at M locations (outputs). This involves inversion of a N × M Transfer Functions (TRF) matrix corresponding to the N input and M output locations. The matrix inversion is affected by both mathematical and physical parameters (ill-conditioned matrix, structural modes, signal noise). Tikhonov regularization is commonly used in inverting an ill-conditioned N × M matrix. A low value of the Tikhonov regularization parameter minimizes the distortion of the original equations while a higher value can minimize error. In practice this introduces an interesting dilemma where obtaining realistic input loads and maintaining accuracy of output are often pitted against each other. A study was conducted using data synthesized from a simply-supported plate structure with known vibration modes with added noise at outputs. The objective of the study was to understand how noise or errors in the output and the Transfer function affect the input. This leads to a more judicious choice of the Tikhonov parameter that can achieve a balance between reducing input loads while preserving desired accuracy of output vibration.
We discuss techniques for efficient local detection of silent data corruption in parallel scientific computations, leveraging physical quantities such as momentum and energy that may be conserved by discretized PDEs. The conserved quantities are analogous to “algorithm-based fault tolerance” checksums for linear algebra but, due to their physical foundation, are applicable to both linear and nonlinear equations and have efficient local updates based on fluxes between subdomains. These physics-based checksums enable precise intermittent detection of errors and recovery by rollback to a checkpoint, with very low overhead when errors are rare. We present applications to both explicit hyperbolic and iterative elliptic (unstructured finite-element) solvers with injected memory bit flips.
It has been predicted that several changes will occur when premixed flames are subjected to the extreme levels of turbulence that can be found in practical combustors. This paper is a review of recent experimental and DNS results that have been obtained for the range of extreme turbulence, and it includes a discussion of cases that agree or disagree with predictions. “Extreme turbulence” is defined to correspond to a turbulent Reynolds number (ReΤ, based on integral scale) that exceeds 2800 or a turbulent Karlovitz number that exceeds 100, for reasons that are discussed in Section 2.1. Several data bases are described that include measurements made at Lund University, the University of Sydney, the University of Michigan and the U.S. Air Force Research Lab. The data bases also include DNS results from Sandia National Laboratory, the University of New South Wales, Newcastle University, the California Institute of Technology and the University of Cambridge. Several major observations are: (a) DNS now can be achieved for a realistic geometry (of the Lund University jet burner) even for extreme turbulence levels, (b) state relations (conditional mean profiles) from DNS and experiments do tend to agree with laminar profiles, at least for methane-air and hydrogen-air reactants that are not preheated, and (c) regime boundaries have been measured and they do not agree with predicted boundaries. These findings indicate that the range of conditions for which flamelet models should be valid is larger than what was previously believed. Additional parameters have been shown to be important; for example, broken reactions occur if the “back-support” is insufficient due to the entrainment of cold gas into the product gas. Turbulent burning velocity measurements have been extended from the previous normalized turbulence levels (u’/SL) of 24 up to a value of 163. Turbulent burning velocities no longer follow the trend predicted by Shchelkin but they tend to follow the trend predicted by Damköhler. The boundary where flamelet broadening begins was measured to occur at ReTaylor = 13.8, which corresponds to an integral scale Reynolds number (ReT) of 2800. This measured regime boundary can be explained by the idea that flame structure is altered when the turbulent diffusivity at the Taylor scale exceeds a critical value, rather than the idea that changes occur when Kolmogorov eddies just fit inside a flamelet. A roadmap to extend DNS to complex chemistry and to higher Reynolds numbers is discussed. Exascale computers, machine learning, adaptive mesh refinement and embedded DNS show promise. Some advances are reviewed that have extended the use of line and planar PLIF and CARS laser diagnostics to studies that consider complex hydrocarbon fuels and harsh environments.
Sikeridis, Dimitrios; Bidram, Ali; Devetsikiotis, Michael; Reno, Matthew J.
Distribution and transmission protection systems are considered vital parts of modern smart grid ecosystems due to their ability to isolate faulted segments and preserve the operation of critical loads. Current protection schemes increasingly utilize cognitive methods to proactively modify their actions according to extreme power system changes. However, the effectiveness and robustness of these information-driven solutions rely entirely on the integrity, authenticity, and confidentiality of the data and control signals exchanged on the underlying relay communication networks. In this paper, we outline a scalable adaptive protection platform for distribution systems, and introduce a novel blockchain-based distributed network architecture to enhance data exchange security among the smart grid protection relays. The proposed mechanism utilizes a tiered blockchain architecture to counter the current technology limitations providing low latency with better scalability. The decentralized nature removes singular points of failure or contamination, enabling direct secure communication between smart grid relays. We also present a security analysis that demonstrates how the proposed framework prohibits any alterations on the blockchain ledger providing integrity and authenticity of the exchanged data (e.g., realtime measurements/relay settings). Finally, the performance of the proposed approach is evaluated through simulation on a blockchain benchmarking framework with the results demonstrating a promising solution for secure smart grid protection system communication.
Fe-Co-2V is a soft ferromagnetic alloy used in electromagnetic applications due to excellent magnetic properties. However, the discontinuous yielding (Luders bands), grain-size-dependent properties (Hall-Petch behavior), and the degree of order/disorder in the Fe-Co-2V alloy makes it difficult to predict the mechanical performance, particularly in abnormal environments such as elevated strain rates and high/low temperatures. Thus, experimental characterization of the high strain rate properties of the Fe-Co-2V alloy is desired, which are used for material model development in numerical simulations. In this study, the high rate tensile response of Fe-Co-2V is investigated with a pulse-shaped Kolsky tension bar over a wide range of strain rates and temperatures. Effects of temperature and strain rate on yield stress, ultimate stress, and ductility are discussed.
[Abstract] To support reduced order modeling of heat transfer for reentry bodies we develop an approximate solution method is identified that provides good estimates for the local wall derivative (and thereby the skin friction and Nusselt numbers) for a wide range of self-similar laminar formulations. These formulations include: Blasius flow, axisymmetric and planar stagnation flows and the Faulkner-Skan flows. The approach utilized is simply an extension of the classical Weyl formulation for the Blasius equation. Using this solution form estimates that naturally represent combined flow behaviors are represented without post-solution interpolation. An important example, namely axisymmetric stagnation equally combined with laminar zero pressure gradient (flat plate) flow, shows a difference of 10% between the pre-solution combination developed here and s simple post-solution arithmetic average. Clearly, the nonlinearity inherent to these solutions prevails in terms of these simple solutions. Compressible extensions to the basic incompressible result are achieved by including a near wall Chapman-Rubesin term making these solutions suitable for adiabatic wall problems. Direct comparison of the wall gradient estimation procedure developed here demonstrates excellent agreement with empirically fit blunt body heat transfer models such as the asymptotically consistent model of Kemp et. al. which are deemed more appropriate than the classical stagnation point scaling approaches.
In this work, we investigated microstructural features of elastomeric foam with the goal of identifying descriptors other than porosity that have a significant effect on the macroscale mechanical response. X-ray computed tomography (XCT) provided three-dimensional images of several flexible polyurethane foam samples prior to mechanical testing. The samples were then compressed to approximately 80% engineering strain. Stereo digital image correlation was used to measure the three-dimensional surface displacement data, from which strain was determined. The strain data, which were calculated with respect to the undeformed coordinates, were then overlaid on the corresponding surface generated from XCT. Heterogeneities in the strain-field were cross-correlated with topological quantities such as pore size distribution. A statistically significant correlation was identified between the distance transform of the pore phase and strain fluctuations.
Recent work has demonstrated that block preconditioning can scalably accelerate the performance of iterative solvers applied to linear systems arising in implicit multiphysics PDE simulations. The idea of block preconditioning is to decompose the system matrix into physical sub-blocks and apply individual specialized scalable solvers to each sub-block. It can be advantageous to block into simpler segregated physics systems or to block by discretization type. This strategy is particularly amenable to multiphysics systems in which existing solvers, such as multilevel methods, can be leveraged for component physics and to problems with disparate discretizations in which scalable monolithic solvers are rare. This work extends our recent work on scalable block preconditioning methods for structure-preserving discretizatons of the Maxwell equations and our previous work in MHD system solvers to the context of multifluid electromagnetic plasma systems. We argue how a block preconditioner can address both the disparate discretization, as well as strongly-coupled off-diagonal physics that produces fast time-scales (e.g. plasma and cyclotron frequencies). We propose a block preconditioner for plasma systems that allows reuse of existing multigrid solvers for different degrees of freedom while capturing important couplings, and demonstrate the algorithmic scalability of this approach at time-scales of interest.
This article describes a parallel implementation of a two-level overlapping Schwarz preconditioner with the GDSW (Generalized Dryja–Smith–Widlund) coarse space described in previous work [12, 10, 15] into the Trilinos framework; cf. [16]. The software is a significant improvement of a previous implementation [12]; see Sec. 4 for results on the improved performance.
Concern about memory errors has been widespread in high-performance computing (HPC) for decades. These concerns have led to significant research on detecting and correcting memory errors to improve performance and provide strong guarantees about the correctness of the memory contents of scientific simulations. However, power concerns and changes in memory architectures threaten the viability of current approaches to protecting memory (e.g., Chipkill). Returning to less protective error-correcting codes (ECC), e.g., single-error correction, double-error detection (SECDED), may increase the frequency of memory errors, including silent data corruption (SDC). SDC has the potential to silently cause applications to produce incorrect results and mislead domain scientists. We propose an approach for exploiting unnecessary bits in pointer values to support encoding the pointer with a Reed-Solomon code. Encoding the pointer allows us to provides strong capabilities for correcting and detecting corruption of pointer values. In this paper, we provide a detailed description of how we can exploit unnecessary pointer bits to store Reed-Solomon parity symbols. We evaluate the performance impacts of this approach and examine the effectiveness of the approach against corruption. Our results demonstrate that encoding and decoding is fast (less than 45 per event) and that the protection it provides is robust (the rate of miscorrection is less than 5% even for significant corruption). The data and analysis presented in this paper demonstrates the power of our approach. It is fast, tunable, requires no additional per-pointer storage resources, and provides robust protection against pointer corruption.
We discuss techniques for efficient local detection of silent data corruption in parallel scientific computations, leveraging physical quantities such as momentum and energy that may be conserved by discretized PDEs. The conserved quantities are analogous to “algorithm-based fault tolerance” checksums for linear algebra but, due to their physical foundation, are applicable to both linear and nonlinear equations and have efficient local updates based on fluxes between subdomains. These physics-based checksums enable precise intermittent detection of errors and recovery by rollback to a checkpoint, with very low overhead when errors are rare. We present applications to both explicit hyperbolic and iterative elliptic (unstructured finite-element) solvers with injected memory bit flips.
Productivity and Sustainability Improvement Planning (PSIP) is a lightweight, iterative workflow that allows software development teams to identify development bottlenecks and track progress to overcome them. In this paper, we present an overview of PSIP and how it compares to other software process improvement (SPI) methodologies, and provide two case studies that describe how the use of PSIP led to successful improvements in team effectiveness and efficiency.
This work presents the application of stabilized continuation to the planar unpowered hypersonic trajectory generation problem using indirect optimal control methods. This scheme is non-iterative and guaranteed to terminate within a finite number of floating point operations-thereby making it well-suited for onboard autonomous operations. This algorithm involves using the stabilized continuation method in multiple stages starting with a “loose” integration tolerance and subsequently ramping up toward a “strict” integration tolerance. An important outcome of this approach is that even when the underlying optimal control problem governing the planar hypersonic trajectory becomes numerically stiff, our studies indicate that the stabilized continuation scheme terminates successfully with a converged solution.
The design, construction, and testing of a high-magnification, long working-distance plenoptic camera is reported. A plenoptic camera uses a microlens array to enable resolution of the spatial and angular information of the incoming light field. Instantaneous images can be numerically refocused and perspective shifted in post-processing to enable threedimensional (3D) resolution of a scene. Prior to this work, most applications of plenoptic imaging were limited to relatively low magnifications (1× or less) or small working distances. Here, a unique system is developed with enables 5× magnification at a working distance of over a quarter meter. Experimental results demonstrate ~25 µm spatial resolution with 3D imaging capabilities. This technology is demonstrated for 3D imaging of the shock structure in a underexpanded, Mach 3.3 free air jet.
Truly predictive numerical simulations can only be obtained by performing Uncertainty Quantification. However, many realistic engineering applications require extremely complex and computationally expensive high-fidelity numerical simulations for their accurate performance characterization. Very often the combination of complex physical models and extreme operative conditions can easily lead to hundreds of uncertain parameters that need to be propagated through high-fidelity codes. Under these circumstances, a single fidelity uncertainty quantification approach, i.e. a workflow that only uses high-fidelity simulations, is unfeasible due to its prohibitive overall computational cost. To overcome this difficulty, in recent years multifidelity strategies emerged and gained popularity. Their core idea is to combine simulations with varying levels of fidelity/accuracy in order to obtain estimators or surrogates that can yield the same accuracy of their single fidelity counterparts at a much lower computational cost. This goal is usually accomplished by defining a priori a sequence of discretization levels or physical modeling assumptions that can be used to decrease the complexity of a numerical model realization and thus its computational cost. Less attention has been dedicated to low-fidelity models that can be built directly from a small number of available high-fidelity simulations. In this work we focus our attention on reduced order models (ROMs). Our main goal in this work is to investigate the combination of multifidelity uncertainty quantification and ROMs in order to evaluate the possibility to obtain an efficient framework for propagating uncertainties through expensive numerical codes. We focus our attention on sampling-based multifidelity approaches, like the multifidelity control variate, and we consider several scenarios for a numerical test problem, namely the Kuramoto-Sivashinsky equation, for which the efficiency of the multifidelity-ROM estimator is compared to the standard (single-fidelity) Monte Carlo approach.
Truly predictive numerical simulations can only be obtained by performing Uncertainty Quantification. However, many realistic engineering applications require extremely complex and computationally expensive high-fidelity numerical simulations for their accurate performance characterization. Very often the combination of complex physical models and extreme operative conditions can easily lead to hundreds of uncertain parameters that need to be propagated through high-fidelity codes. Under these circumstances, a single fidelity uncertainty quantification approach, i.e. a workflow that only uses high-fidelity simulations, is unfeasible due to its prohibitive overall computational cost. To overcome this difficulty, in recent years multifidelity strategies emerged and gained popularity. Their core idea is to combine simulations with varying levels of fidelity/accuracy in order to obtain estimators or surrogates that can yield the same accuracy of their single fidelity counterparts at a much lower computational cost. This goal is usually accomplished by defining a priori a sequence of discretization levels or physical modeling assumptions that can be used to decrease the complexity of a numerical model realization and thus its computational cost. Less attention has been dedicated to low-fidelity models that can be built directly from a small number of available high-fidelity simulations. In this work we focus our attention on reduced order models (ROMs). Our main goal in this work is to investigate the combination of multifidelity uncertainty quantification and ROMs in order to evaluate the possibility to obtain an efficient framework for propagating uncertainties through expensive numerical codes. We focus our attention on sampling-based multifidelity approaches, like the multifidelity control variate, and we consider several scenarios for a numerical test problem, namely the Kuramoto-Sivashinsky equation, for which the efficiency of the multifidelity-ROM estimator is compared to the standard (single-fidelity) Monte Carlo approach.
As computer architectures are rapidly evolving (e.g. those designed for exascale), multiple portability frameworks have been developed to avoid new architecture-specific development and tuning. However, portability frameworks depend on compilers for auto-vectorization and may lack support for explicit vectorization on heterogeneous platforms. Alternatively, programmers can use intrinsics-based primitives to achieve more efficient vectorization, but the lack of a gpu back-end for these primitives makes such code non-portable. A unified, portable, Single Instruction Multiple Data (simd) primitive proposed in this work, allows intrinsics-based vectorization on cpus and many-core architectures such as Intel Knights Landing (knl), and also facilitates Single Instruction Multiple Threads (simt) based execution on gpus. This unified primitive, coupled with the Kokkos portability ecosystem, makes it possible to develop explicitly vectorized code, which is portable across heterogeneous platforms. The new simd primitive is used on different architectures to test the performance boost against hard-to-auto-vectorize baseline, to measure the overhead against efficiently vectroized baseline, and to evaluate the new feature called the “logical vector length” (lvl). The simd primitive provides portability across cpus and gpus without any performance degradation being observed experimentally.