Sintering is one of the oldest processes used by man to manufacture materials dating as far back as 12,000 BC. While it is an ancient process, it is also necessary for many modern technologies such a multilayered ceramic packages, wireless communication devices, and many others. The process consists of thermally treating a powder or compact at a temperature below the melting point of the main constituent, for the purpose of increasing its strength by bonding together of the particles. During sintering, the individual particles bond, the pore space between particles is eliminated, the resulting component can shrinks by as much as 30 to 50% by volume, and it can distort its shape tremendously. Being able to control and predict the shrinkage and shape distortions during sintering has been the goal of much research in material science. And it has been achieved to varying degrees by many. The object of this project was to develop models that could simulate sintering at the mesoscale and at the macroscale to more accurately predict the overall shrinkage and shape distortions in engineering components. The mesoscale model simulates microstructural evolution during sintering by modeling grain growth, pore migration and coarsening, and vacancy formation, diffusion and annihilation. In addition to studying microstructure, these simulation can be used to generate the constitutive equations describing shrinkage and deformation during sintering. These constitutive equations are used by continuum finite element simulations to predict the overall shrinkage and shape distortions of a sintering crystalline powder compact. Both models will be presented. Application of these models to study sintering will be demonstrated and discussed. Finally, the limitations of these models will be reviewed.
Blastwalls are often assumed to be the answer for facility protection from malevolent explosive assault, particularly from large vehicle bombs (LVB's). The assumption is made that the blastwall, if it is built strong enough to survive, will provide substantial protection to facilities and people on the side opposite the LVB. This paper will demonstrate through computer simulations and experimental data the behavior of explosively induced air blasts during interaction with blastwalls. It will be shown that air blasts can effectively wrap around and over blastwalls. Significant pressure reduction can be expected on the downstream side of the blastwall but substantial pressure will continue to propagate. The effectiveness of the blastwall to reduce blast overpressure depends on the geometry of the blastwall and the location of the explosive relative to the blastwall.
The purpose of this study was to investigate the impact of instructions on aircraft visual inspection performance and strategy. Forty-two inspectors from industry were asked to perform inspections of six areas of a Boeing 737. Six different instruction versions were developed for each inspection task, varying in the number and type of directed inspections. The amount of time spent inspecting, the number of calls made, and the number of the feedback calls detected all varied widely across the inspectors. However, inspectors who used instructions with a higher number of directed inspections referred to the instructions more often during and after the task, and found a higher percentage of a selected set of feedback cracks than inspectors using other instruction versions. This suggests that specific instructions can help overall inspection performance, not just performance on the defects specified. Further, instructions were shown to change the way an inspector approaches a task.
This paper proposes a method for predicting the complexity of meshing computer aided design (CAD) geometries with unstructured, hexahedral, finite elements. Meshing complexity refers to the relative level of effort required to generate a valid finite element mesh on a given CAD geometry. A function is proposed to approximate the meshing complexity for single part CAD models. The function is dependent on a user defined element size as well as on data extracted from the geometry and topology of the CAD part. Several geometry and topology measures are proposed, which both characterize the shape of the CAD part and detect configurations that complicate mesh generation. Based on a test suite of CAD models, the function is demonstrated to be accurate within a certain range of error. The solution proposed here is intended to provide managers and users of meshing software a method of predicting the difficulty in meshing a CAD model. This will enable them to make decisions about model simplification and analysis approaches prior to mesh generation.
Memory may be the only system component that is more commoditized than a microprocessor. To simultaneously exploit this and address the impending memory wall, processing in memory (PIM) research efforts are considering ways to move processing into memory without significantly increasing the cost of the memory. As such, PIM devices may become the basis for future commodity clusters. Although these PIM devices may leverage new computational paradigms such as hardware support for multi-threading and traveling threads, they must provide support for legacy programming models if they are to supplant commodity clusters. This paper presents a prototype implementation of MPI over a traveling thread mechanism called parcels. A performance analysis indicates that the direct hardware support of a traveling thread model can lead to an efficient, lightweight MPI implementation.
Mechanisms for enhanced low-dose-rate sensitivity are described. In these mechanisms, bimolecular reactions dominate the kinetics at high dose rates thereby causing a sub-linear dependence on total dose, and this leads to a dose-rate dependence. These bimolecular mechanisms include electron-hole recombination, hydrogen recapture at hydrogen source sites, and hydrogen dimerization to form hydrogen molecules. The essence of each of these mechanisms is the dominance of the bimolecular reactions over the radiolysis reaction at high dose rates. However, at low dose rates, the radiolysis reaction dominates leading to a maximum effect of the radiation.
Properties of relevance for the equation of state for a high-density glass are discussed. We review the effects of failure waves, comminuted phase, and compaction on the validity of the Mie-Grueneisen EOS. The specific heat and the Grueneisen parameter at standard conditions for a {rho}{sub 0} = 5.085 g/cm{sup 3} glass ('Glass A') is then estimated to be 522 mJ/g/K and 0.1-0.3, respectively. The latter value is substantially smaller than the value of 2.1751 given in the SESAME tables for a high-density glass with {rho}{sub 0} = 5.46 g/cm{sup 3}. The present unusual value of the Grueneisen parameter is confirmed from the volume dependence determined from fitting the Mie-Grueneisen EOS to shock data in Ref. [2].
Submovements are hypothesized building blocks of human movement, discrete ballistic movements of which more complex movements are composed. Using a novel algorithm, submovements were extracted from the point-to-point movements of 41 persons recovering from stroke. Analysis of the extracted submovements showed that, over the course of therapy, patients' submovements tended to increase in peak speed and duration. The number of submovements employed to produce a given movement decreased. The time between the peaks of adjacent submovements decreased for inpatients (those less than 1 month post-stroke), but not for outpatients (those greater than 12 months post-stroke) as a group. Submovements became more overlapped for all patients, but more markedly for inpatients. The strength and consistency with which it quantified patients' recovery indicates that analysis of submovement overlap might be a useful tool for measuring learning or other changes in motor behavior in future human movement studies.
Proposed for publication in Coordinated & Multiple Views in Exploratory Visualization, Special Issue of Information Visualization Journal, Vol 2 No. 4, Palgrave/Macmillan.
PIM (Processor in Memory) architectures are being proposed for future supercomputers, because they reduce the problems that SMP MMPs have with latency. However, they do not meet the SMP MPP balance factors. Being relatively processor rich and memory starved, it is unclear whether an ASCI application could run on them, either as-is or with recoding. The KBA (Koch-Baker-Alcouffe) algorithm (Koch, 1992) for particle transport (radiation transport) is shown not to fit on PIMs as written. When redesigned with a 3-D allocation of cells to PIMs, the resulting algorithm is projected to execute an order of magnitude faster and more cost-effectively than the KBA algorithm, albeit with high initial hardware costs.
We present a series of electronic structure calculations that demonstrate a mechanism for spontaneous ionization of hydrogen at the Si-SiO{sub 2} interface. Specifically, we show that an isolated neutral hydrogen atom will spontaneously give up its charge and bond to a threefold coordinated oxygen atom. We refer to this entity as a proton. We have calculated the potential surface and found it to be entirely attractive. In contrast, hydrogen molecules will not undergo an analogous reaction. We relate these calculations both to proton generation experiments and to hydrogen plasma experiments.
Density functional theory is used to predict workfunctions, {psi}. For relaxed clean W(1 0 0), the local density approximation (LDA) agrees with experiment better than the newer generalized gradient approximation, probably due to the surface electron self-energy. The large Ba metallic radius indicates it covers W(1 0 0) at about 0.5 monolayer (ML). However, Ba{sup 2+}, O{sup 2-}, and metallic W all have similar radii. Thus 1 ML of BaO (one BaO unit for each two W atoms) produces minimum strain, indicating commensurate interfaces. BaO (1 ML) and Ba (1/2 ML) have the same {psi} to within 0.02 V, so at these coverages reduction or oxidation is not important. Due to greater chemical activity of ScO vs. highly ionic BaO, when mixing the latter with this suboxide of scandia, the overlayer always has BaO as the top layer and ScO as the second layer. The BaO/ScO bilayer has a rocksalt structure, suggesting high stability. In the series BaO/ScO/, BaO/YO/, and BaO/LaO/W(1 0 0), the latter has a remarkably low {psi} of 1.3 V (LDA), but 2 ML of rocksalt BaO also has {psi} at 1.3 V. We suggest BaO (1 ML) does not exist and that it is worthwhile to attempt the direct synthesis and study of BaO (2 ML) and BaO/LaO.
Several methods for preparing well equilibrated melts of long chains polymers are studied. We show that the standard method in which one starts with an ensemble of chains with the correct end-to-end distance arranged randomly in the simulation cell and introduces the excluded volume rapidly, leads to deformation on short length scales. This deformation is strongest for long chains and relaxes only after the chains have moved their own size. Two methods are shown to overcome this local deformation of the chains. One method is to first pre-pack the Gaussian chains, which reduces the density fluctuations in the system, followed by a gradual introduction of the excluded volume. The second method is a double-bridging algorithm in which new bonds are formed across a pair of chains, creating two new chains each substantially different from the original. We demonstrate the effectiveness of these methods for a linear bead spring polymer model with both zero and nonzero bending stiffness, however the methods are applicable to more complex architectures such as branched and star polymer.
This document provides a user manual for the SGOPT software library. SGOPT is a C++ class library for nonlinear optimization. This library uses an object-oriented design that allows the software to be extended to a new problem domains. Furthermore, this library was designed to that the interface is straightforward while providing flexibility to allow new algorithms to be easily added to this library. The SGOPT library has been used by several software projects at Sandia, and it is integrated into the DAKOTA design and analysis toolkit. This report provides a high-level description of the optimization algorithms provided by SGOPT and describes the C++ class hierarchy in which they are implemented. Finally, installation instructions are included.
An effort is underway at Sandia National Laboratories to develop a library of algorithms to search for potential interactions between surfaces represented by analytic and discretized topological entities. This effort is also developing algorithms to determine forces due to these interactions for transient dynamics applications. This document describes the Application Programming Interface (API) for the ACME (Algorithms for Contact in a Multiphysics Environment) library.
The ability for scientific simulation software to detect and recover from errors and failures of supporting hardware and software layers is becoming more important due to the pressure to shift from large, specialized multi-million dollar ASCI computing platforms to smaller, less expensive interconnected machines consisting of off-the-shelf hardware. As evidenced by the CPlant{trademark} experiences, fault tolerance can be necessary even on such a homogeneous system and may also prove useful in the next generation of ASCI platforms. This report describes a research effort intended to study, implement, and test the feasibility of various fault tolerance mechanisms controlled at the simulation code level. Errors and failures would be detected by underlying software layers, communicated to the application through a convenient interface, and then handled by the simulation code itself. Targeted faults included corrupt communication messages, processor node dropouts, and unacceptable slowdown of service from processing nodes. Recovery techniques such as re-sending communication messages and dynamic reallocation of failing processor nodes were considered. However, most fault tolerance mechanisms rely on underlying software layers which were discovered to be lacking to such a degree that mechanisms at the application level could not be implemented. This research effort has been postponed and shifted to these supporting layers.
