Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We have studied the electron emission from one of the polycyclic aromatic hydrocarbon (PAH) molecules namely, fluorene (C13H10), upon 3.5 MeV/u Si8+ ion impact. The experimentally measured absolute double differential cross sections (DDCS) are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) model and the first Born approximation including correct boundary conditions (CB1). The measurements are carried out in the ejected e-energy range of 1 eV-400 eV and in the angular range of 20 -160 . We have obtained the single differential and the total cross sections (TCSs) of e-emission as well. The CB1 calculation largely underestimates the data. The CDW-EIS model, which is applied for the PAH molecule for the first time, provides an overall better agreement with the double differential, single differential and TCS data. The DDCS data for fluorene has also been compared with that for CH4 molecule, at a few angles. The forward-backward angular asymmetry shows a relatively flatter distribution compared to the theoretical predictions. The contribution due to the giant plasmon resonance could not be clearly observed except a mild indication in the asymmetry parameter. The angular distribution of the carbon KLL Auger electron cross section shows certain variations. The study of the KLL hyper-satellite component indicates the double K-ionization cross section is about 8.6% of the single K-ionization one.

Several Department of Energy (DOE) facilities have materials stored in robust, welded, stainless - steel containers with presumed fire - induced pressure response behaviors. Lack of test data related to fire exposure requires conservative safety analysis assumptions for container response at these facilities. This conservatism can in turn result in the implementation of challenging operational restrictions with costly nuclear safety controls. To help address this issue for sites that store DOE 3013 stainless - steel containers, a series of ten tests were undertaken at Sandia National Laboratories. The goal of this test series was to obtain the response behavior for various configurations of DOE 3013 containers with various payload compositions when exposed to one of two ASTM fire conditions. Key parameters measured in the test series included identification of failure - specific characteristics such as pressure, temperature, and whether or not a vessel was breached during a test . Numerous failure - specific characteristics were identified from the ten tests. This report describes the implementation and execution of the test series performed to identify these failure - specific characteristics. Discussions on the test configurations, payload compositions, thermal insults, and experimental setups are presented. Test results in terms of pressurization and vessel breach (or no - breach) are presented along with corresponding discussions for each test.

Abstract not provided.

A collection of x-ray computed tomography scans of specimens from the Museum of Southwestern Biology.

Abstract not provided.

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.

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.

We describe the opportunities and challenges we faced when developing SIGNAL, an experimental wargame that was deployed as a distributed wargame.

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.

Laboratory experiments typically test opacity models by measuring spectrally resolved transmission of a sample using bright backlight radiation. A potential problem is that any unaccounted background signal contaminating the spectrum will artificially reduce the inferred opacity. Methods developed to measure background signals in opacity experiments at the Sandia Z facility are discussed. Preliminary measurements indicate that backgrounds are 9%-11% of the backlight signal at wavelengths less than 10 Å. Background is thus a relatively modest correction for all Z opacity data published to date. Future work will determine how important background is at longer wavelengths.

Transition metal dichalcogenides, such as MoS2, are known to undergo a structural phase transformation as well as a change in the electronic conductivity upon Li intercalation. These properties make them candidates for charge-tunable ion-insertion materials that could be used in electrochemical devices. In this paper, we study the phase stability and electronic structure of H and T' Li-intercalated MoX2 bilayers with X = S, Se, or Te. Using first-principles calculations in combination with classical and machine learning modeling approaches, we find that the H phase is more stable at low Li concentration for all X and the critical Li concentration at which the T'→H transition occurs decreases with increasing mass of X. Furthermore, the relative free energy of the two phases becomes less sensitive to Li insertion with increasing atomic mass of the chalcogen atom X. While the electronic conductivity increases with increasing ion concentration at low concentrations, we do not observe a (positive) conductivity jump at the phase transition from H to T'.

Sandia National Laboratories’ Global Chemical and Biological Security program (SNL/GCBS) has an on-going mission to enhance laboratory biosafety and biosecurity. SNL/GCBS has developed robust biosafety and biosecurity risk assessment methodologies and tools to aid laboratories seeking to implement biosecurity as advocated in the recently released World Health Organization’s “Biorisk Management: Laboratory Biosecurity Guidance.” BioRAM Security (BioRAM) is one of these tools. This tool was designed to complement the “Laboratory Biosecurity Handbook” written by Ren Salerno and Jennifer Gaudioso (Salerno, 2007). BioRAM was originally developed under the Sandia National Laboratories Laboratory Directed Research and Development (LDRD) process and has been updated with support from the US Department of Health Centers for Disease Control Select Agent Program, the US Department of State Biosecurity Engagement Program, and by funding from Sandia National Laboratories’ Associated Laboratories Directors.

The development of reduced-order models remains an active research area, despite advances in computational resources. The present work develops a novel order-reduction approach that is designed to incorporate isolated regions that contain, for example, nonlinearitites or accumulating damage. The approach is designed to use global modes of the overall system response, which are then naturally coupled to the response in the isolated region of interest. Two examples are provided to demonstrate both the accuracy and the computational efficiency of the proposed approach. The performance of this approach is compared to the exact response corresponding to a finite element simulation for the chosen problems. In addition, the accuracy and computational efficiency are shown relative to a standard Galerkin reduction based on the linear normal modes. It is found that the proposed reduction offer computational efficiency comparable to a Galerkin reduction, but more accurately represents the response of the system when both are compared to the finite element simulation.

The generalized singular value decomposition (GSVD) is a valuable tool that has many applications in computational science. However, computing the GSVD for large-scale problems is challenging. Motivated by applications in hyper-differential sensitivity analysis (HDSA), we propose new randomized algorithms for computing the GSVD which use randomized subspace iteration and weighted QR factorization. Detailed error analysis is given which provides insight into the accuracy of the algorithms and the choice of the algorithmic parameters. We demonstrate the performance of our algorithms on test matrices and a large-scale model problem where HDSA is used to study subsurface flow.

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

For a wide range of PDEs, the discontinuous Petrov–Galerkin (DPG) methodology of Demkowicz and Gopalakrishnan provides discrete stability starting from a coarse mesh and minimization of the residual in a user-controlled norm, among other appealing features. Research on DPG for transient problems has mainly focused on spacetime discretizations, which has theoretical advantages, but practical costs for computations and software implementations. The sole examination of time-stepping DPG formulations was performed by Führer, Heuer, and Gupta, who applied Rothe's method to an ultraweak formulation of the heat equation to develop an implicit time-stepping scheme; their work emphasized theoretical results, including error estimates in time and space. In the present work, we follow Führer, Heuer, and Gupta in examining the heat equation; our focus is on numerical experiments, examining the stability and accuracy of several formulations, including primal as well as ultraweak, and explicit as well as implicit and Crank–Nicolson time-stepping schemes. We are additionally interested in communication-avoiding algorithms, and we therefore include a highly experimental formulation that places all the trace terms on the right-hand side of the equation.

MTTKRP is the bottleneck operation in algorithms used to compute the CP tensor decomposition. For sparse tensors, utilizing the compressed sparse fibers (CSF) storage format and the CSF-oriented MTTKRP algorithms is important for both memory and computational efficiency on distributed-memory architectures. Existing intelligent tensor partitioning models assume the computational cost of MTTKRP to be proportional to the total number of nonzeros in the tensor. However, this is not the case for the CSF-oriented MTTKRP on distributed-memory architectures. We outline two deficiencies of nonzero-based intelligent partitioning models when CSF-oriented MTTKRP operations are performed locally: failure to encode processors' computational loads and increase in total computation due to fiber fragmentation. We focus on existing fine-grain hypergraph model and propose a novel vertex weighting scheme that enables this model encode correct computational loads of processors. We also propose to augment the fine-grain model by fiber nets for reducing the increase in total computational load via minimizing fiber fragmentation. In this way, the proposed model encodes minimizing the load of the bottleneck processor. Parallel experiments with real-world sparse tensors on up to 1024 processors prove the validity of the outlined deficiencies and demonstrate the merit of our proposed improvements in terms of parallel runtimes.

This article proposes a method of predicting the influence of random load behavior on the dynamics of dc microgrids and distribution systems. This is accomplished by combining stochastic load models and deterministic microgrid models. Together, these elements constitute a stochastic hybrid system. The resulting model enables straightforward calculation of dynamic state moments, which are used to assess the probability of desirable operating conditions. Specific consideration is given to systems based on the dual active bridge (DAB) topology. Bounds are derived for the probability of zero voltage switching (ZVS) in DAB converters. A simple example is presented to demonstrate how these bounds may be used to improve ZVS performance as an optimization problem. Predictions of state moment dynamics and ZVS probability assessments are verified through comparisons to Monte Carlo simulations.

Abstract not provided.

Numerous MELCOR modeling improvements and analyses have been performed in the time since the severe accidents at Fukushima Daiichi Nuclear Power Station that occurred in March 2011. This report briefly summarizes the related accident reconstruction and uncertainty analysis efforts. It further discusses a number of potential pursuits to further advance MELCOR modeling and analysis of the severe accidents at Fukushima Daiichi and severe accident modeling in general. Proposed paths forward include further enhancements to identified MELCOR models primarily impacting core degradation calculations, and continued application of uncertainty analysis methods to improve model performance and a develop deeper understanding of severe accident progression.

Abstract not provided.

Abstract not provided.

Abstract not provided.