The cost of reliability: Iterative linear solvers and reactive fault tolerance
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Proceedings of the 6th International Workshop on Programming Models and Applications for Multicores and Manycores, PMAM 2015
The Bulk Synchronous Parallel programming model is showing performance limitations at high processor counts. We propose over-decomposition of the domain, operated on as tasks, to smooth out utilization of the computing resource, in particular the node interconnect and processing cores, and hide intra- and inter-node data movement. Our approach maintains the existing coding style commonly employed in computational science and engineering applications. Although we show improved performance on existing computers, up to 131,072 processor cores, the effectiveness of this approach on expected future architectures will require the continued evolution of capabilities throughout the codesign stack. Success then will not only result in decreased time to solution, but would also make better use of the hardware capabilities and reduce power and energy requirements, while fundamentally maintaining the current code configuration strategy.
Created at Sandia National Laboratories, the Portable File Format (PFF) allows binary data transfer across computer platforms. Although this capability is supported by many other formats, PFF files are still in use at Sandia, particularly in pulsed power research. This report provides detailed PFF specifications for accessing data without relying on legacy code.
Abstract not provided.
Abstract not provided.
The Sandia Data Archive (SDA) format is a specific implementation of the HDF5 (Hierarchal Data Format version 5) standard. The format was developed for storing data in a universally accessible manner. SDA files may contain one or more data records, each associated with a distinct text label. Primitive records provide basic data storage, while compound records support more elaborate grouping. External records allow text/binary files to be carried inside an archive and later recovered. This report documents version 1.0 of the SDA standard. The information provided here is sufficient for reading from and writing to an archive. Although the format was original designed for use in MATLAB, broader use is encouraged.
Abstract not provided.
Overview: A Pu shot scheduled for July 17 on the Z machine at SNL was cancelled this past summer. The LiF windows on the Pu targets were cracked during assembly because of configuration changes. Sandia management concluded that continuing with this experiment would present an unacceptable level of risk to the facility and possibly to the workers. In this report, we document the events that occurred which led to this decision and also present some lessons learned and plans and procedures put in place to reduce the likelihood of another such occurrence. The changes and this memorandum reflect the thinking of subject matter experts at both LANL and SNL. These changes represent significant improvements in both communication protocols and quality of the hardware assemblies.
Abstract not provided.
Abstract not provided.
Team Sandia California (Team H) used the Sandia code SIERRA Solid Mechanics: Implicit (SIERRA SM) to model the SFC2 challenge problem. SIERRA SM is a Lagrangian, three-dimensional, implicit code for the analysis of solids and structures. It contains a versatile library of continuum and structural elements, and an extensive library of material models. For all SFC2 related simulations, our team used Q1P0, 8 node hexahedral elements with element side lengths on the order 0.175 mm in failure regions. To model crack initiation and failure, element death removed elements from the simulation according to a continuum damage model. SIERRA SM’s implicit dynamics, implemented with an HHT time integration scheme for numerical damping [1], was used to model the unstable failure modes of the models. We chose SIERRA SM’s isotropic Elasto Viscoplastic material model for our simulations because it contains most of the physics required to accurately model the SFC2 challenge problem such as the flexibility to include temperature and rate dependence for a material.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Journal of Computational Science
Incorrect computer hardware behavior may corrupt intermediate computations in numerical algorithms, possibly resulting in incorrect answers. Prior work models misbehaving hardware by randomly flipping bits in memory. We start by accepting this premise, and present an analytic model for the error introduced by a bit flip in an IEEE 754 floating-point number. We then relate this finding to the linear algebra concepts of normalization and matrix equilibration. In particular, we present a case study illustrating that normalizing both vector inputs of a dot product minimizes the probability of a single bit flip causing a large error in the dot product's result. Moreover, the absolute error is either less than one or very large, which allows detection of large errors. Then, we apply this to the GMRES iterative solver. We count all possible errors that can be introduced through faults in arithmetic in the computationally intensive orthogonalization phase of GMRES, and show that when the matrix is equilibrated, the absolute error is bounded above by one.
Applied Optics
We find for infrared wavelengths that there are broad ranges of particle sizes and refractive indices that represent fog and rain, where circular polarization can persist to longer ranges than linear polarization. Using polarization tracking Monte Carlo simulations for varying particle size, wavelength, and refractive index, we show that, for specific scene parameters, circular polarization outperforms linear polarization in maintaining the illuminating polarization state for large optical depths. This enhancement with circular polarization can be exploited to improve range and target detection in obscurant environments that are important in many critical sensing applications. Initially, researchers employed polarizationdiscriminating schemes, often using linearly polarized active illumination, to further distinguish target signals from the background noise. More recently, researchers have investigated circular polarization as a means to separate signal from noise even more. Specifically, we quantify both linearly and circularly polarized active illumination and show here that circular polarization persists better than linear for radiation fog in the short-wave infrared, for advection fog in the short-wave and long-wave infrared, and large particle sizes of Sahara dust around the 4 μmwavelength. Conversely, we quantify where linear polarization persists better than circular polarization for some limited particle sizes of radiation fog in the long-wave infrared, small particle sizes of Sahara dust for wavelengths of 9-10.5 μm, and large particle sizes of Sahara dust through the 8-11 μm wavelength range in the long-wave infrared.
Proceedings of SPIE - The International Society for Optical Engineering
We present simulation results that show circularly polarized light persists through scattering environments better than linearly polarized light. Specifically, we show persistence is enhanced through many scattering events in an environment with a size parameter representative of advection fog at infrared wavelengths. Utilizing polarization tracking Monte Carlo simulations we show a larger persistence benefit for circular polarization versus linear polarization for both forward and backscattered photons. We show the evolution of the incident polarization states after various scattering events which highlight the mechanism leading to circular polarization's superior persistence.
NSTI: Advanced Materials - TechConnect Briefs 2015
SNL has developed a series of ionic-liquid electrolytes with accompanying non- Aqueous compatible membranes and flow cell designs for improved energy density redox flow batteries targeted to support increasing demands for stationary energy storage. The new electrolytes yield a higher energy density by chemically incorporating an electro- Active transition metal element into the solvent's molecular formula. Although ionic liquids have higher viscosities than conventional non- Aqueous electrolytes, they are promising for higher energy densities due to higher metal concentrations and wider voltage windows. We have addressed high viscosity by developing new materials through careful ligand and anion selection. We have also developed tunable membranes for non- Aqueous compatibility and rapid laboratory-scale prototyping to quickly screen materials and cell designs. We are projecting a four-fold improvement in energy density over the next two years.