Publications

Numerical simulations of perforation in steel plates involve the treatment of material failure during the perforation process. One way to model physical material separation is to delete failed elements from the analysis based on an appropriate failure criterion. Different algorithms were used in different transient finite element codes to delete failed elements. This investigation compares the results of PRONTO 2D and LS-DYNA2D codes for a specific steel plate perforation problem. Influences of the deletion algorithms on material parameters are discussed.

A framework for coupled elastoplastic and damage theories is developed, following a rigorous thermodynamic procedure. This framework is sufficiently general to include anisotropic plasticity and damage formulations. Both the plastic yield and damage functions are constructed using homogeneous functions of degree one. The principle of maximum dissipation or maximum entropy production is used to derive the evolution relations together with the loading and unloading conditions. In addition, the convexity of the undamaging elastic domain is shown. For plasticity the resulting evolution of the plastic strains corresponds to an associative flow. This general framework is shown to be sufficiently general to describe several popular theories for both plasticity and damage. Limitations of some existing damage theories are discussed.

Continuum damage theories describe the progressive reduction in stiffness and strength of brittle materials resulting from the initiation and growth of microcracks and microvoids. When brittle materials are loaded into the nonlinear regime, they often exhibit localized zones of intense deformation and the eventual formation of macrocracks. Criteria for diffuse and discontinuous bifurcations have previously been developed and used to study the initiation of necking and localization in elastic-plastic materials. In this investigation, the same bifurcation criteria are applied to continuum damage theories. Since the bifurcation criteria depend on the fourth-order tangent modulus tensor, the first step in this investigation is the derivation of the tangent modulus tensor for a general continuum damage theory. An eigenanalysis of the symmetric part of the tangent modulus tensor is then shown to fully characterize the potential diffuse and discontinuous bifurcations associated with a given continuum damage theory.

The motivation behind exchange I.2 of the USSR/US exchange program of cooperation for magnetic confinement fusion is to more closely simulate tokamak disruptions with a variety of plasma devices within the Soviet Union and the United States and to characterize the effect these simulated disruptions have on candidate PFC materials. Earlier work conducted in the Soviet Union by a team of Soviet and American researchers showed ablation of graphites exposed to a disruption like heat flux from a plasma flow was significantly less than that previously expected [J.M. Gahl et al., Proc. ICFRM-5, J. Nucl. Mater. 191-194 (1992) 454]. Work has continued and results from recent work at the University of New Mexico are in general agreement with earlier results from the Soviet Union. New results from work in the United States and the Soviet Union will be presented. © 1992 Elsevier Science Publishers B.V. All rights reserved.

This paper describes a current research program at Sandia National Laboratories whereby magnetic stripes are produced through the use of a new particle rotation technology. This new process allows the stripes to be produced in bulk and then held in a latent state so that they may be encoded at a later date. Since particle rotation is less dependent on the type of magnetic particle used, very high coercivity particles could provide a way to increase both magnetic tamper-resistance and accidental erasure protection. This research was initially funded by the Department of Energy, Office of Safeguard and Security as a portion of their Science and Technology Base Development, Advanced Security Concepts program. Current program funding is being provided by Sandia National Laboratories as part of their Laboratory Directed Research and Development program.

Given a graph G = (V, E) where each vertex v {element_of} V is assigned a weight w(v) and each edge e {element_of} E is assigned a cost c(e), the quotient of a cut partitioning the vertices of V into sets S and {bar S} is c(S, {bar S})/min{l_brace}w(S), w(S){r_brace}, where c(S, {bar S}) is the sum of the costs of the edges crossing the cut and w(S) and w({bar S}) are the sum of the weights of the vertices in S and {bar S}, respectively. The problem of finding a cut whose quotient is minimum for a graph has in recent years attracted considerable attention, due in large part to the work of Rao and Leighton and Rao. They have shown that an algorithm (exact or approximation) for the minimum-quotient-cut problem can be used to obtain an approximation algorithm for the more famous minimumb-balanced-cut problem, which requires finding a cut (S,{bar S}) minimizing c(S,{bar S}) subject to the constraint bW {le} w(S) {le} (1 {minus} b)W, where W is the total vertex weight and b is some fixed balance in the range 0 < b {le} {1/2}. Unfortunately, the minimum-quotient-cut problem is strongly NP-hard for general graphs, and the best polynomial-time approximation algorithm known for the general problem guarantees only a cut whose quotient is at mostO(lg n) times optimal, where n is the size of the graph. However, for planar graphs, the minimum-quotient-cut problem appears more tractable, as Rao has developed several efficient approximation algorithms for the planar version of the problem capable of finding a cut whose quotient is at most some constant times optimal. In this paper, we improve Rao`s algorithms, both in terms of accuracy and speed. As our first result, we present two pseudopolynomial-time exact algorithms for the planar minimum-quotient-cut problem. As Rao`s most accurate approximation algorithm for the problem -- also a pseudopolynomial-time algorithm -- guarantees only a 1.5-times-optimal cut, our algorithms represent a significant advance.

Advection-dominated flows occur widely in the transport of groundwater contaminants, the movements of fluids in enhanced oil recovery projects, and many other contexts. In numerical models of such flows, adaptive local grid refinement is a conceptually attractive approach for resolving the sharp fronts or layers that tend to characterize the solutions. However, this approach can be difficult to implement in practice. A domain decomposition method developed by Bramble, Ewing, Pasciak, and Schatz, known as the BEPS method, overcomes many of the difficulties. We demonstrate the applicability of the iterative BEPS ideas to finite-element collocation on trial spaces of piecewise Hermite bicubics. The resulting scheme allows one to refine selected parts of a spatial grid without destroying algebraic efficiencies associated with the original coarse grid. We apply the method to two dimensional time-dependent advection-diffusion problems.

The ability to design gating systems that reliably feed and support investment castings is often the result of ``cut-and-try`` methodology. Factors such as hot tearing, porosity, cold shuts, misruns, and shrink are defects often corrected by several empirical gating design iterations. Sandia National Laboratories is developing rules that aid in removing the uncertainty involved in the design of gating systems for investment castings. In this work, gating geometries used for filling of thin walled investment cast 17-4PH stainless steel flat plates were investigated. A full factorial experiment evaluating the influence of metal pour temperature, mold preheat temperature, and mold channel thickness were conducted for orientations that filled a horizontal flat plate from the edge. A single wedge gate geometry was used for the edge-gated configuration. Thermocouples placed along the top of the mold recorded metal front temperatures, and a real-time x-ray imaging system tracked the fluid flow behavior during filling of the casting. Data from these experiments were used to determine the terminal fill volumes and terminal fill times for each gate design.

This paper studied the behavior of retained system poles and transmission zeros in a control design model when the model is truncated. The sensitivity of the transmission zeros due to the tuncation of system dynamics was analytically obtained. The sensitivity of system poles to the truncation of system dynamics was shown to be zero as expected. The effects of actuator-sensor type and location was also studied. The results were illustrated with two example problems. The effect of transmission zero shifts in control design models and the controllers designed from them was illustrated with an example.

A telerobotic system demonstration was developed for the Department of Energy`s Accident Response group to highlight the applications of telerobotic vehicles to accident site inspection. The proof-of- principle system employs two mobile robots, Dixie and RAYBOT, to inspect a simulated accident site. Both robots are controlled serially from a single driving station, allowing an operator to take advantage of having multiple robots at the scene. The telerobotic system is described and some of the advantages of having more than one robot present are discussed. Future plans for the system are also presented.

This report will present a brief overview of the transient dynamics capabilities at Sandia National Laboratories, with an emphasis on recent new developments and current research. In addition, the Sandia National Laboratories (SNL) Engineering Analysis Code Access System (SEACAS), which is a collection of structural and thermal codes and utilities used by analysts at SNL, will be described. The SEACAS system includes pre- and post-processing codes, analysis codes, database translation codes, support libraries, Unix shell scripts for execution, and an installation system. SEACAS is used at SNL on a daily basis as a production, research, and development system for the engineering analysts and code developers. Over the past year, approximately 190 days of CPU time have been used by SEACAS codes on jobs running from a few seconds up to two and one-half days of CPU time. SEACAS is running on several different systems at SNL including Cray Unicos, Hewlett Packard HP-UX, Digital Equipment Ultrix, and Sun SunOS. An overview of SEACAS, including a short description of the codes in the system, will be presented. Abstracts and references for the codes are listed at the end of the report.

Computational mechanics simulation capability via the finite element method is being integrated into the FASTCAST project to allow realistic analyses of investment casting problems. Commercial and in-house software is being coupled to new, solid model based mesh generation capabilities to provide improved access to fluid, thermal and structural simulations. These simulations are being used for the validation of complex gating designs and the study of fundamental problems in casting.

A simple, approximate model of parachute inflation is described. The model is based on the traditional, practical treatment of the fluid resistance of rigid bodies in nonsteady flow, with appropriate extensions to accommodate the change in canopy inflated shape. Correlations for the steady drag and steady radial force as functions of the inflated radius are required as input to the dynamic model. In a novel approach, the radial force is expressed in terms of easily obtainable drag and reefing fine tension measurements. A series of wind tunnel experiments provides the needed correlations. Coefficients associated with the added mass of fluid are evaluated by calibrating the model against an extensive and reliable set of flight data. A parameter is introduced which appears to universally govern the strong dependence of the axial added mass coefficient on motion history. Through comparisons with flight data, the model is shown to realistically predict inflation forces for ribbon and ringslot canopies over a wide range of sizes and deployment conditions.

This document presents the quality assurance (QA) procedures for Parameter Selection and Expert Judgment Panels used by the performance Assessment Department of Sandia National Laboratories, which directly supports the Waste Isolation Pilot Plant (WIPP). Parameter Selection QA procedures described herein will be incorporated into the general Performance Assessment Quality Assurance Procedures, QAP 2-3; the Expert Judgment Panel procedures will be incorporated into QAP 2-6. Both sets of procedures will apply to all Sandia and Sandia contractor activities related to performance assessment (except where the contractor has its own approved QA procedures). This report presents the philosophy behind the QA procedures, provides the standards adopted for performance assessment Parameter Selection and Expert Judgment Panels, and discusses the implementation of these standards.

This 1991 report contains monitoring data from routine radiological and nonradiological environmental surveillance activities. Summaries of significant environmental compliance programs in progress such as National Environmental Policy Act (NEPA) documentation, environmental permits, environmental restoration (ER), and various waste management programs for Sandia National Laboratories in Albuquerque (SNL, Albuquerque) are included. The maximum offsite dose impact was calculated to be 1.3 {times} 10{sup {minus}3} mrem. The total population within a 50-mile radius of SNL, Albuquerque, received a collective dose of 0.53 person-rem during 1991 from SNL, Albuquerque, operations. As in the previous year, the 1991 operations at SNL, Albuquerque, had no discernible impact on the general public or on the environment.

This preliminary study analyzes the atmospheric entry of a solid core nuclear thermal rocket (NTR) engine under three accidental entry scenarios. Depending on the scenario, results of the analysis showed that, without external thermal protection, an aluminum pressure vessel will fail at altitudes ranging 25 to 73 km. subsequent release the core materials occurs. The graphitic based core materials will undergo partial ablation, with the percent mass loss depending on the geometry of the fuel elements. A carbon-phenolic thermal protection system was sized to prevent pressure vessel aerothermal failure. It was found to increase the mass of the NTR by approximately 15 percent.

This report examines containment requirements for spent-fuel transport containers that are transported under normal and hypothetical accident conditions. A methodology is described that estimates the probability of rod failure and the quantity of radioactive material released from breached rods. This methodology characterizes the dynamic environment of the cask and its contents and deterministically models the peak stresses that are induced in spent-fuel cladding by the mechanical and thermal dynamic environments. The peak stresses are evaluated in relation to probabilistic failure criteria for generated or preexisting ductile tearing and material fractures at cracks partially through the wall in fuel rods. Activity concentrations in the cask cavity are predicted from estimates of the fraction of gases, volatiles, and fuel fines that are released when the rod cladding is breached. Containment requirements based on the source term are calculated in terms of maximum permissible volumetric leak rates from the cask. Calculations are included for representative cask designs.

This report describes a mechanism for compiling the functional language SISAL for Sandia`s Epsilon-2 hybrid dataflow machine. The strategy couples the front-end of the standard SISAL compiler (which generates a data dependence graph intermediate form called IF1) with an optimizing code-generator for Epsilon-2. The Epsilon-2 code-generator is the back-end of a compiler for the functional language Id. It translates a data dependence graph intermediate form called Program Graphs into Epsilon-2 machine code. This report describes a translation path from IF1 graphs to Program Graphs. This report also comments on the relative merits of the IF1 and Program Graph representations.

This memorandum is a synopsis of the description and operation of the equipment used and the events accuring during the calibration of gauges on the vacuum station over the range of 0.0001 to 650 torr.

This report summarizes the environmental surveillance activities conducted by the US Environmental Protection Agency (EPA) and Reynolds Electrical and Engineering Company (REECO) for the Tonopah Test Range (TTR) operated by Sandia National Laboratories (SNL). Other environmental compliance programs such as the National Environmental Policy Act of 1969 (NEPA), environmental permits, environmental restoration, and waste management programs are also included. The 1991 SNL, TTR, operations had no discernible impact on the general public or the environment. This report 3-s prepared for the US Department of Energy (DOE) in compliance with DOE Order 5400.1.

Through a program sponsored by the US Department of Energy (DOE), Cummins Power Generation, Inc. (CPG) and Sandia National Laboratories (SNL) have entered into a joint venture to develop and commercialize economically competitive dish-Stirling systems for remote power applications. Sixteen systems, representing three generations of technology, will be developed, fielded, and tested in the Dish-Stirling Joint Venture Program (JVP). The JVP is funded equally by a consortium led by CPG and by the DOE. After completion of the program, CPG`s commercialization effort will continue with limited production expected to start in 1996. In this paper, the program plan and the technology used in the JVP are outlined. ne current status of the key system components, and the initial results of a system optimization study including current cost and performance estimates, are also provided.

Research in recent years has demonstrated the efficient use of solar thermal energy for driving endothermic chemical reforming reactions in which hydrocarbons are reacted to form syngas. Closed-loop reforming/methanation systems can be used for storage and transport of process heat and for short-term storage for peaking power generation. Open-loop reforming and gasification systems can be used for direct fuel production; for production of syngas feedstock for further processing to bulk ammonia, hydrogen, and liquid fuels; and for destruction of hazardous organic materials. To help identify the most promising areas for future development of this technology, we discuss in this paper the market potential of these applications.

FASTCAST is a Sandia National Laboratories program to produce investment cast prototypical hardware faster by integrating experimental and computational technologies into the casting process. FASTCAST uses the finite element method to characterize the metal flow and solidification processes to reduce uncertainty in the mold design. For the casting process to benefit from finite element analysis, analysis results must be available in a very short time frame. By focusing on the bottleneck of finite element model creation, automated mesh generation can drastically reduce the time span between geometry definition (design) and accurate analysis results. The increased availability of analysis results will diminish the need for trial and error approaches to acquiring production worthy mold and gating systems for investment casting. The CUBIT meshing tool kit is being developed to address the need for rapid mesh generation. CUBIT is being designed to effectively automate the generation of quadrilateral and hexahedral elements. It is a solid-modeler based, two- and three-dimensional preprocessor that prepares solid models for finite element analysis. CUBIT contains several meshing algorithms including two- and three-dimensional mapping, two- and three-dimensional paving (patented), and a general two and one-half dimensional sweeper based upon the plastering algorithm. This paper describes progress in the development of the CUBIT meshing toolkit.

The purpose of this NUREG is to provide technical information on the major components of entry control systems: identity verifiers, weapons detectors, explosives detectors, and special nuclear material (SNM) detectors. For each type of device, information is presented on principles of operation, hardware features, recommended installation, testing methods, and operational procedures. Applications to personnel, handcarried packages, bulk items, and vehicles are addressed.

We have investigated the chemistry of Cu(hfac){sub 2}, (hfac)Cu(VTMS), (hfac)Cu(2-butyne), and hfach on a Pt(111) surface. In contrast to what is observed on copper surfaces. Cu(hfac)2 and hfach lead to the formation of distinctly different adsorbed hfac species on Pt(111). This shows the importance of the copper atoms themselves in determining the surface chemistry of copper {beta}-diketonate CVD precursors. The hfac species on Pt(111) are considerably less stable than hfac on copper, suggesting that unimolecular decomposition may lead to impurity incorporation in the interfacial region when copper is deposited onto a more reactive substrate. In situ CVD studies with Cu(I) {beta}-diketonates show that the bimolecular disproportionation reaction leading to copper CVD is favored over unimolecular precursor decomposition at pressures above approximately 10{sup {minus}5} torr.