Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. We also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.

Presentation of damage detection and experimental design approaches for structures under frequency-domain dynamics

Abstract not provided.

Abstract not provided.

Abstract not provided.

This work details the reconfiguration of the 4.5 m Gigahertz Transverse Electromagnetic test facility at Sandia National Laboratories to operate in accordance with the RS105 (radiated susceptibility) test from MIL-STD-461 representing a high-altitude electromagnetic pulse. This reconfiguration involved removal of the existing continuous wave source and connecting both a high voltage feed and a coaxial feed housing the Marx bank pulser. Marx control settings were calibrated for several voltage levels across two pulsers, and position-dependent measurements of the peak electric field were taken throughout the test volume for each pulser. The results showed field uniformity and purity across the test volume comparable to continuous wave operations, and field peaks were measured from 1.63 kV/m to 54.8 kV/m, with maximum capabilities expected to exceed 100 kV/m. Some challenges in consistent pulser operations at lower Marx bank voltages and high frequency reflections in the system were identified for future capability improvements.

Abstract not provided.

Overview of work completed as a summer intern in the CCD.

Abstract not provided.

Abstract not provided.

The U.S. Department of Energy/National Nuclear Security Administration (DOE/NNSA) and National Technology & Engineering Solutions of Sandia, LLC (NTESS), the management and operating contractor for Sandia National Laboratories/California (SNL/CA), has prepared this addendum to Soil Sampling Results for Closure of a Portion of Solid Waste Management Unit #16 to report the results of additional soil sampling relating to the closure of a portion of Solid Waste Management Unit (SWMU) #16. This additional sampling was in response to a request by the San Francisco Bay Regional Water Quality Control Board (SFRWQCB) in their letters dated February 16 and August 18, 2022 relating to the detection of the benzidine above the defined project action level in a soil sample collected adjacent to the sanitary sewer line in borehole BH-056 (SFRWQCB, 2022A; 2022b).

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Z Fundamental Science (ZFS) Program is intended to provide access to the Z machine and its diagnostics for high energy density (HED) experiments in collaboration with a broad community of academic, industrial, and national laboratory research interests. ZFS experiments on the Z Facility focus on conducting fundamental research in HED science and help provide research experience necessary to maintain and grow the HED community, especially through involvement of researchers from academia. This report serves as an executive summary of the ZFS Program and provides a succinct synopsis of the history of the ZFS Program, metrics and impacts of the Program, as well as a brief list of the most impactful publications that have resulted from the various ZFS Projects relevant to laboratory astrophysics, plasma physics, and planetary physics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Recent work on the development of integrated thermographic phosphors and digital image correlation (TP+DIC) for combined thermal–mechanical measurements has revealed the need for a flexible, stretchable phosphor coating for metal surfaces. Herein, we coat stainless steel substrates with a polymer-based phosphor ink in a DIC speckle pattern and demonstrate that the ink remains well bonded under substrate deformation. In contrast, a binderless phosphor DIC coating produced via aerosol deposition (AD) partially debonded from the substrate. DIC calculations reveal that the strain on the ink coating matches the strain on the substrate within 4% error at the highest substrate loads (0.05 mm/mm applied substrate strain), while the strain on the AD coating remains near 0 mm/mm as the substrate deforms. Spectrally resolved emission from the phosphor is measured through the transparent binder throughout testing, and the ratio method is used to infer temperature with an uncertainty of 1.7 °C. This phosphor ink coating will allow for accurate, non-contact strain and temperature measurements of a deforming surface.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Micro ribonucleic acids (miRNA) give our immune systems the ability to recognize viruses and other pathogens by their complementary single-stranded RNA (ssRNA) produced in the reproduction of the pathogen in our cells. When miRNA of a specific sequence is detected in a cell sample, it can be assumed that the immune system is activated and attempting to track down the infection. This pathway can be utilized to diagnose infection from a pathogen before the individual even develops symptoms, aiding in early disease detection and proper treatment. One of the ways that we can detect miRNA is through an assay of clustered regularly interspaced short palindromic repeats or “CRISPR” and the bacterial protein Cas13a. This report details discoveries made while attempting to optimize this assay for miRNA detection. After looking at several different factors within the assay, it was determined that some factors, such as reporter type and metallic ion concentration, are more impactful on the overall assay sensitivity than other factors, such as the overall concentration of Cas13a, CRISPR RNA (crRNA), or ssRNA reporter. It was also discovered that different sequences with different lengths require renewed optimization efforts, as each target has a unique binding affinity determined by the sequence length and composition. This information is crucial in the development of point of care molecular detection devices as they become sensitive enough to identify pathogens before they spread.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Co-Author Voon Hui Lai, Australian National University, will be attending IUGG and presenting the work. The session of interested is "S05 - Advances in Earthquake and Explosion Monitoring Using Distributed Acoustic Sensing​"

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present the characterization of several atmospheric aerosol analogs in a tabletop chamber and an analysis of how the concentration of NaCl present in these aerosols influences their bulk optical properties. Atmospheric aerosols (e.g., fog and haze) degrade optical signal via light–aerosol interactions causing scattering and absorption, which can be described by Mie theory. This attenuation is a function of the size distribution and number concentration of droplets in the light path. These properties are influenced by ambient conditions and the droplet’s composition, as described by Köhler theory. It is therefore possible to tune the wavelength-dependent bulk optical properties of an aerosol by controlling droplet composition. We present experimentation wherein we generated multiple microphysically and optically distinct atmospheric aerosol analogs using salt water solutions with varying concentrations of NaCl. The results demonstrate that changing the NaCl concentration has a clear and predictable impact on the microphysical and optical properties of the aerosol

Abstract not provided.

Abstract not provided.

Abstract not provided.

Dynamic compression studies have been used to study the nucleation kinetics of water to ice VII for decades. Diagnostics such as photon Doppler velocimetry, transmission loss, and imaging have been used to measure pressure/density, and phase fraction, while temperature has remained the difficult thermodynamic property to quantify. In this work, we measured pressure/density and implemented a diagnostic to measure the temperature. In doing so the temperature shows quasi-isentropically compressed liquid water forms ice at pressures below the previously defined metastable limit, and the liquid phase is not hypercoooled as previously thought above that limit. Instead, the latent heat raises the temperature to the liquid-ice-VII melt line, where it remains with increasing pressure. We propose a hypothesis to corroborate these results with previous work on dynamic compression freezing. These results provide constraints for nucleation models, and suggest this technique be used to investigate phase transitions in other materials.

This report describes research and development (R&D) activities conducted during Fiscal Year 2023 (FY23) in the Advanced Fuels and Advanced Reactor Waste Streams Strategies work package in the Spent Fuel Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). This report is focused on evaluating and cataloguing Advanced Reactor Spent Nuclear Fuel (AR SNF) and Advanced Reactor Waste Streams (ARWS) and creating Back-end Nuclear Fuel Cycle (BENFC) strategies for their disposition. The R&D team for this report is comprised of researchers from Sandia National Laboratories and Enviro Nuclear Services, LLC.

This is an investigation on two experimental datasets of laminar hypersonic flows, over a double-cone geometry, acquired in Calspan—University at Buffalo Research Center’s Large Energy National Shock (LENS)-XX expansion tunnel. These datasets have yet to be modeled accurately. A previous paper suggested that this could partly be due to mis-specified inlet conditions. The authors of this paper solved a Bayesian inverse problem to infer the inlet conditions of the LENS-XX test section and found that in one case they lay outside the uncertainty bounds specified in the experimental dataset. However, the inference was performed using approximate surrogate models. Here in this paper, the experimental datasets are revisited and inversions for the tunnel test-section inlet conditions are performed with a Navier–Stokes simulator. The inversion is deterministic and can provide uncertainty bounds on the inlet conditions under a Gaussian assumption. It was found that deterministic inversion yields inlet conditions that do not agree with what was stated in the experiments. An a posteriori method is also presented to check the validity of the Gaussian assumption for the posterior distribution. This paper contributes to ongoing work on the assessment of datasets from challenging experiments conducted in extreme environments, where the experimental apparatus is pushed to the margins of its design and performance envelopes.

Work accomplished: Collected and compared historic data for the 1993 Rock Valley earthquake sequence; Compared preliminary and prior location work from different location algorithms, phase pick sets, station constellations, and velocity models; Selected a common set of stations that could be used across all location methods for consistency; Reviewed 8 different sets of phase picks and converged on a single, reviewed set of picks for all common stations; Evaluated four pre-existing regional velocity models and incorporated new and preliminary results for five new velocity models that provide information on the very shallow (< 2km) structure near station RTPP; Compared location results from different methods while using the common sets of picks, stations, and velocity models

Real-time monitoring of a research nuclear reactor, a system in which all generated power is dissipated to the environment, can be performed via analysis of the heat rejection from the cooling system. Given an inlet water temperature and flow rate, the reactor power can be well-approximated from the outlet water temperature; however, the instrumentation to measure outlet conditions may not be robust or accurate. If we know how a cooling tower performs from historical data, but cannot measure the outlet temperature, a mathematical representation of the system can be inverted to obtain the outlet water temperature that describes the cooling capacity. Unfortunately, model inversion processes are computationally expensive. To address this, an artificial neural network (ANN) is implemented to assess the performance of a multi-cell cooling tower for a nuclear reactor. This approach leverages the Merkel model to obtain an extensive data set describing performance of the cooling tower cells throughout a wide array of potential operating conditions. The Merkel model is expressed as a function of four parameters: the inlet and outlet water temperatures, inlet air wet bulb temperature, and ratio of liquid-to-gas mass flow rates (L/G), which together provide a non-dimensional number indicative of cooling tower performance, called the Merkel integral. Computing a 4-dimensional data structure that describes finite combinations of the Merkel integral, an inverse model is then generated using an ANN to determine the cell outlet water temperature from the other three model parameters along with the computed Merkel integral. Compared to traditional model inversion methods, the ANN reduces the computational time by approximately 4 orders of magnitude, with effectively no sacrifice to solution accuracy, and could be applied for different cooling towers in the event the performance curve is known. Finally, three use cases of the ANN are then reviewed: (1) determining the cell outlet water temperatures when gas flow at rated conditions (GFRC) is known, (2) performing the prior case without knowledge of the GRFC, and (3) assessing performance differences between the individual tower cells.

Electrothermal instability plays an important role in applications of current-driven metal, creating striations (which seed the magneto-Rayleigh-Taylor instability) and filaments (which provide a more rapid path to plasma formation). However, the initial formation of both structures is not well understood. Simulations show for the first time how a commonly occurring isolated defect transforms into the larger striation and filament, through a feedback loop connecting current and electrical conductivity. Simulations have been experimentally validated using defect-driven self-emission patterns.

Using three-dimensional (3D) magnetohydrodynamic simulations, we study how a pit on a metal surface evolves when driven by intense electrical current density j. Redistribution of j around the pit initiates a feedback loop: j both reacts to and alters the electrical conductivity σ, through Joule heating and hydrodynamic expansion, so that j and σ are constantly in flux. Thus, the pit transforms into larger striation and filament structures predicted by the electrothermal instability theory. Both structures are important in applications of current-driven metal: Here, the striation constitutes a density perturbation that can seed the magneto-Rayleigh-Taylor instability, while the filament provides a more rapid path to plasma formation, through 3D j redistribution. Simulations predict distinctive self-emission patterns, thus allowing for experimental observation and comparison.

We present a highly diagonal “split-well resonant-phonon” (SWRP) active region design for GaAs/Al_0.3Ga_0.7As terahertz quantum cascade lasers (THz-QCLs). Negative differential resistance is observed at room temperature, which indicates the suppression of thermally activated leakage channels. The overlap between the doped region and the active level states is reduced relative to that of the split-well direct-phonon (SWDP) design. The energy gap between the lower laser level (LLL) and the injector is kept at 36 meV, enabling a fast depopulation of the LLL. Within this work, we investigated the temperature performance and potential of this structure.

Gas intercalation into clay interlayers may result in hydrogen loss in the geological storage of hydrogen; a phenomenon that has not been fully understood and quantified. Here we use metadynamics molecular simulations to calculate the free energy landscape of H2 intercalation into montmorillonite interlayers and the H2 solubility in the confined water; in comparison with results obtained for CO2. The results indicate that H2 intercalation into hydrated interlayers is thermodynamically unfavorable while CO2 intercalation can be favorable. H2 solubility in hydrated clay interlayers is in the same order of magnitude as that in bulk water and therefore no over-solubility effect due to nanoconfinement is observed - in striking contrast with CO2. These results indicate that H2 loss and leakage through hydrated interlayers due to intercalation in a subsurface storage system, if any, is limited.

Sapphire (Al₂O₃) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350-500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. We also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

Previous studies of the cantilevered pipeline conveying fluid system have included motion-limiting constraints in the form of trilinear springs. While this is desirable in experimental scenarios, it may not be representative of real-world applications. Therefore, here, this study focuses on multi-segmented motion-limiting constraints. As this type of motion-limiting constraint has not been investigated with a cantilevered pipeline system, a wide variety of outer and inner constraint stiffness and constraint gap sizes are investigated in this study to gain a comprehensive understanding of how the multi-segmented constraints affect the dynamics of the cantilevered pipeline. In this effort, bifurcation diagrams, phase portraits, Poincare maps, time histories, and power spectra are used to investigate the dynamics of the system, and the fluid flow speeds where dynamic characteristics are considered. In general, it is found that critical flow speeds like when the pipe sticks in the constraints are reduced as the constraint stiffnesses are increased. Additionally, the sticking flow speed occurred at lower flow speeds as the gap sizes of the inner and outer constraints decrease, and a larger constraint offset results in a smaller inner gap size leading to critical behaviors occurring at earlier flow speeds.

ALEGRA is a multiphysics finite-element shock hydrodynamics code, under development at Sandia National Laboratories since 1990. Fully coupled multiphysics capabilities include transient magnetics, magnetohydrodynamics, electromechanics, and radiation transport. Importantly, ALEGRA is used to study hypervelocity impact, pulsed power devices, and radiation effects. The breadth of physics represented in ALEGRA is outlined here, along with simulated results for a selected hypervelocity impact experiment.

Fluorinated graphite materials are of interest for an assortment of applications and can be synthesized under a variety of synthetic conditions from many different types of carbon. Due to such variations, structural disorders in the form of defects and polymorphism are often present. Here, we investigate the impact of local structural variations on the C-F bond dissociation energies (BDEs) in carbon-based fluoride materials using density functional theory (DFT) computational methods. Employing fluorographene (FG) cluster models, we determine the impact of different C-F bonding configurations in the core of each platelet on the equilibrium BDEs for each C-F bond. The introduction of structural disorder decreases the first C-F BDE by approximately 1 eV compared to the canonical arrangement of axial C-F bonds ordered as in a network of cyclohexane “chairs”. Variability of calculated BDEs among the different polymorphs decreases upon subsequent F removal. Common structural tendencies of the adiabatic defluorination pathways for each polymorph are identified. Our analysis suggests that at F/C ratios near 1.0, disorder in the local structure can play a significant role in the energetics of the initial carbon fluoride defluorination and that the influence of this configurational disorder diminishes with decreasing F/C ratios.

Out-of-equilibrium electrochemical reaction mechanisms are notoriously difficult to characterize. However, such reactions are critical for a range of technological applications. For instance, in metal-ion batteries, spontaneous electrolyte degradation controls electrode passivation and battery cycle life. Here, to improve our ability to elucidate electrochemical reactivity, we for the first time combine computational chemical reaction network (CRN) analysis based on density functional theory (DFT) and differential electrochemical mass spectroscopy (DEMS) to study gas evolution from a model Mg-ion battery electrolyte-magnesium bistriflimide (Mg(TFSI)2) dissolved in diglyme (G2). Automated CRN analysis allows for the facile interpretation of DEMS data, revealing H2O, C2H4, and CH3OH as major products of G2 decomposition. These findings are further explained by identifying elementary mechanisms using DFT. While TFSI-is reactive at Mg electrodes, we find that it does not meaningfully contribute to gas evolution. The combined theoretical-experimental approach developed here provides a means to effectively predict electrolyte decomposition products and pathways when initially unknown.

Bayesian analysis enables flexible and rigorous definition of statistical model assumptions with well-characterized propagation of uncertainties and resulting inferences for single-shot, repeated, or even cross-platform data. This approach has a strong history of application to a variety of problems in physical sciences ranging from inference of particle mass from multi-source high-energy particle data to analysis of black-hole characteristics from gravitational wave observations. The recent adoption of Bayesian statistics for analysis and design of high-energy density physics (HEDP) and inertial confinement fusion (ICF) experiments has provided invaluable gains in expert understanding and experiment performance. In this Review, we discuss the basic theory and practical application of the Bayesian statistics framework. We highlight a variety of studies from the HEDP and ICF literature, demonstrating the power of this technique. Due to the computational complexity of multi-physics models needed to analyze HEDP and ICF experiments, Bayesian inference is often not computationally tractable. Two sections are devoted to a review of statistical approximations, efficient inference algorithms, and data-driven methods, such as deep-learning and dimensionality reduction, which play a significant role in enabling use of the Bayesian framework. We provide additional discussion of various applications of Bayesian and machine learning methods that appear to be sparse in the HEDP and ICF literature constituting possible next steps for the community. We conclude by highlighting community needs, the resolution of which will improve trust in data-driven methods that have proven critical for accelerating the design and discovery cycle in many application areas.

The Source Physics Experiments (SPE) were designed to improve our physics-based understanding of explosion sources for the purposes of nuclear test monitoring. Phase I consisted of 6 chemical explosions in the Climax Stock Granite of the Nevada National Security Site (NNSS), while Phase II consisted of 4 explosions in a contrasting dry alluvium geology (DAG) in Yucca Flat, providing essential data in various media and emplacement conditions to further modeling efforts. For Phase III, the Rock Valley Direct Comparison (RVDC) seeks to directly compare earthquake and explosion source types. An unusually shallow series of events in 1993 along the Rock Valley Fault Zone in the southeastern portion of the NNSS has been targeted for this direct comparison. Depth ranges for the events, previously estimated to be less than 3 km, is achievable by modern drilling techniques and accessibility to the epicentral locations would require minimal improvements to the infrastructure. The events providing this unique opportunity for direct comparison are the focus of this report.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The atomization, mixing, combustion and emissions characteristics of aviation fuels were measured using a novel approach based on a non-continuous injection scheme called the single-hole atomizer (SHA). High-speed microscopy revealed differences between fuels in terms of evaporation and mixing regimes over conditions relevant to modern and next generation aero-engine combustors. Measurements of liquid and vapor penetration, mixing fields, combustion and emissions metrics (ignition delay, lift-off length, PAH formation, soot mass) highlighted the effects of fuels and combustor conditions. The experimental results are being leveraged to adjust and validate chemical and CFD models. Detailed analysis of sampled soot showed subtle differences in soot morphology between fuels. The results revealed the presence of contaminants potentially affecting surface chemistry and the nucleation propensity of water droplets on particles. Chemical mechanisms for NJFCP A-2, C-1 and C-4 showed good performance over a large parameter space. Spray breakup at relight conditions is vastly different from the atomization observed at high pressure. CFD simulations of the SHA target conditions confirmed the good behavior of the C-1 kinetic mechanism. The simulations support the strong relationship between low and high temperature reactions. New altitude chamber facility to enable detailed characterization of the heterogeneous nucleation process of water on aerosol particles.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Marine hydrokinetic devices, such as wave energy converters (WECs), can unlock untapped energy from the ocean's currents and waves. Acoustic impact assessments are required to ensure that the noise these devices generate will not negatively impact marine life, and accurate modeling of noise provides an a priori means to viably perform this assessment. We present a case study of the PacWave South site, a WEC testing site off the coast of Newport, Oregon, demonstrating the use of ParAcousti, an open-source hydroacoustic propagator tool, to model noise from an array of 28 WECs in a 3-dimensional (3-D) realistic marine environment. Sound pressure levels are computed from the modeled 3-D grid of pressure over time, which we use to predict marine mammal acoustic impact metrics (AIMs). We combine two AIMs, signal to noise ratio and sensation level, into a new metric, the effective signal level (ESL), which is a function of propagated sound, background noise levels, and hearing thresholds for marine species and is evaluated across 1/3 octave frequency intervals. The ESL model can be used to predict and quantify the potential impact of an anthropogenic signal on the health and behavior of a marine mammal species throughout the 3-D simulation area.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In order to impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

Abstract not provided.

Abstract not provided.

In recent years, the Engine Combustion Network (ECN) has developed as a worldwide reference for understanding and describing engine combustion processes, successfully bringing together experimental and numerical efforts. Since experiments and numerical simulations both target the same boundary conditions, an accurate characterization of the stratified environment that is inevitably present in experimental facilities is required. The difference between the core-, and pressure-derived bulk-temperature of pre-burn combustion vessels has been addressed in various previous publications. Additionally, thermocouple measurements have provided initial data on the boundary layer close to the injector nozzle, showing a transition to reduced ambient temperatures. The conditions at the start of fuel injection influence physicochemical properties of a fuel spray, including near nozzle mixing, heat release computations, and combustion parameters. To address the temperature stratification in more detail, thermocouple measurements at larger distances from the spray axis have been conducted. Both the temperature field prior to the pre-combustion event that preconditions the high-temperature, high-pressure ambient, as well as the stratification at the moment of fuel injection were studied. To reveal the cold boundary layer near the injector with a better spatial resolution, Rayleigh scattering experiments and thermocouple measurements at various distances close to the nozzle have been carried out. The impact of the boundary layers and temperature stratification are illustrated and quantified using numerical simulations at Spray A conditions. Next to a reference simulation with a uniform temperature field, six different stratified temperature distributions have been generated. These distributions were based on the mean experimental temperature superimposed by a randomized variance, again derived from the experiments. The results showed that an asymmetric flame structure arises in the computed results when the temperature stratification input is used. In these predictions, first-stage ignition is advanced by 24μs, while second-stage ignition is delayed by 11μs. At the same time a lift-off length difference between the top and the bottom of up to 1.1 mm is observed. Furthermore, the lift-off length is less stable over time. Given the shown dependency, the temperature data is made available along with the vessel geometry data as a recommended basis for future numerical simulations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Discharge of sodium coolant into containment from a sodium-cooled fast reactor vessel can occur in the event of a pipe leak or break. In this situation, some of the liquid sodium droplets discharged from the coolant system will react with oxygen in the air before reaching the containment. This phase of the event is normally termed the sodium spray fire phase. Unreacted sodium droplets pool on the containment floor where continued reaction with containment atmospheric oxygen occurs. This phase of the event is normally termed the sodium pool fire phase. Both phases of these sodium-oxygen reactions (or fires) are important to model because of the heat addition and aerosol generation that occur. Any fission products trapped in the sodium coolant may also be released during this progression of events, which if released from containment could pose a health risk to workers and the public. The paper describes progress of an international collaborative research in the area of the sodium fire modeling in the sodium-cooled fast reactors between the United States and Japan under the framework of the Civil Nuclear Energy Research and Development Working Group. In this collaboration between Sandia National Laboratories and Japan Atomic Energy Agency, the validation basis for and modeling capabilities of sodium spray and pool fires in MELCOR of Sandia National Laboratories and SPHINCS of Japan Atomic Energy Agency are being enhanced. This study documents MELCOR and SPHINCS sodium pool fire model validation exercises against the JAEA's sodium pool fire experiments, F7-1 and F7-2. The proposed enhancement of the sodium pool fire models in MELCOR through addition of thermal hydraulic and sodium spreading models that enable a better representation of experimental results is also described.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Effective diversion of surge currents is vital to prevent unwanted damage to sensitive electronics. Among the most successful and efficient strategies relies on a dielectric stimulated arc breakdown mechanism with high permittivity ceramic granules in a spark-gap geometry. Although generally regarded as a self-healing process, substantial energy deposition may occur that, over time, diminishes the ability to withstand repeated electrical assaults. We investigate the susceptibility of lead–magnesium–niobate–lead titanate (PMN–PT) granule microstructure and composition changes following many exposures to high voltage impulses resulting in arc breakdown. Scanning electron microscopy and energy-dispersive spectroscopy mapping reveal a broad range of thermal and mechanical defects entailing thermal reduction of constituent PMN–PT metal ions and recasting due to rapid eruption of volatile species. Additionally, evidence of local melting and microcracking are apparent that can have deleterious impact on the proper function of the granules, namely, the ability to concentrate electric fields across air gaps to establish and sustain discharge pathways. We propose that the localized nature of damage and stochasticity associated with the dielectric stimulated breakdown mechanism may allow granules to maintain functionality provided no permanent conduction paths are established.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Accurate measurements of the arrival times of seismic waves are crucial for seismological analyses such as robust locations of earthquakes, characterization of seismic sources, and high-fidelity imaging of the Earth’s interior. However, these travel-time measurements can sometimes be contaminated by timing errors at the stations which record this data. In this study, we apply a classical approach, based on identifying time-dependence in measured body wave arrival times, to identify these timing errors in a dataset on the order of 107 individual measurements. We find timing deviations at a subset of the stations in our dataset and document the temporal location, extent, and severity of these errors, finding errors at 83 stations, and impacting ~100,000 measurements. This catalog of deviations may enable future investigators to obtain a more accurate dataset through the implementation of quality control measures to eliminate the contaminated data we have identified.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Sandia National Laboratories, in California (SNL/CA) is a research and development facility, owned by the U.S. Department of Energy’s National Nuclear Security Administration agency (DOE/NNSA). The laboratory is located in the City of Livermore (the City) and is comprised of approximately 410 acres. The SNL/CA facility is operated by National Technology and Engineering Solutions of Sandia, LLC (NTESS) under a contract with the DOE/NNSA. The DOE/ NNSA’s Sandia Field Office (SFO) oversees the operations of the site. North of the SNL/CA facility is the Lawrence Livermore National Laboratory (LLNL), in which SNL/CA’s sewer system combines with before discharging to the City’s Publicly Owned Treatment Works (POTW) for final treatment and processing. The City’s POTW authorizes the wastewater discharge from SNL/CA via the assigned Wastewater Discharge Permit #1251 (the Permit), which is issued to the DOE/NNSA’s main office for Sandia National Laboratories, located in New Mexico (SNL/NM). The Monitoring and Reporting Condition 2.B of the Permit requires compliance with the semiannual reporting requirements contained in federal categorical pretreatment standards regulations (40 CFR 403.12). These regulations set numerical limits on the concentration of pollutants allowed to discharge from certain categories of industrial processes. This report is submitted to the City to satisfy this reporting requirement.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The notion of plane shock waves is a macroscopic, very fruitful idealization of near discontinuous disturbance propagating at supersonic speed. Such a picture is comparable to the picture of shorelines seen from a very high altitude. When viewed at the grain scale where the structure of solids is inherently heterogeneous and stochastic, features of shock waves are non-laminar and field variables, such as particle velocity and pressure, fluctuate. This paper reviews select aspects of such fluctuating nonequilibrium features of plane shock waves in solids with focus on grain scale phenomena and raises the need for a paradigm change to achieve a deeper understanding of plane shock waves in solids.

Abstract not provided.