Publications

This report describes Welch's method for computing Power Spectral Densities (PSDs). We first describe the bandpass filter method which uses filtering, squaring, and averaging operations to estimate a PSD. Second, we delineate the relationship of Welch's method to the bandpass filter method. Third, the frequency domain signal-to-noise ratio for a sine wave in white noise is derived. This derivation includes the computation of the noise floor due to quantization noise. The signal-to-noise ratio and noise flood depend on the FFT length and window. Fourth, the variance the Welch's PSD is discussed via chi-square random variables and degrees of freedom. This report contains many examples, figures and tables to illustrate the concepts. 26 refs.

Four expert-judgment teams have developed analyses delineating possible future societies in the next 10,000 years in the vicinity of the Waste Isolation Pilot Plant (WIPP). Expert-judgment analysis was used to address the question of future societies because neither experimentation, observation, nor modeling can resolve such uncertainties. Each of the four, four-member teams, comprised of individuals with expertise in the physical, social, or political sciences, developed detailed qualitative assessments of possible future societies. These assessments include detailed discussions of the underlying physical and societal factors that would influence society and the likely modes of human-intrusion at the WIPP, as well as the probabilities of intrusion. Technological development, population growth, economic development, conservation of information, persistence of government control, and mitigation of danger from nuclear waste were the factors the teams believed to be most important. Likely modes of human-intrusion were categorized as excavation, disposal/storage, tunneling, drilling, and offsite activities. Each team also developed quantitative assessments by providing probabilities of various alternative futures, of inadvertent human intrusion, and in some cases, of particular modes of intrusion. The information created throughout this study will be used in conjunction with other types of information, including experimental data, calculations from physical principles and computer models, and perhaps other judgments, as input to performance assessment.'' The more qualitative results of this study will be used as input to another expert panel considering markers to deter inadvertent human intrusion at the WIPP.

This paper summarizes the results of aging, condition monitoring, and accident testing of Class 1E cables used in nuclear power generating stations. Three sets of cables were aged for up to 9 months under simultaneous thermal ({approx_equal} 100{degrees}C) and radiation ({approx_equal}0.10 kGy/hr) conditions. After the aging, the cables were exposed to a simulated accident consisting of high dose rate irradiation ({approx_equal}6 kGy/hr) followed by a high temperature steam exposure. A fourth set of cables, which were unaged, were also exposed to the accident conditions. The cables that were aged for 3 months and then accident tested were subsequently exposed to a high temperature steam fragility test (up to 400{degrees}C), while the cables that were aged for 6 months and then accident tested were subsequently exposed to a 1000-hour submergence test in a chemical solution. The results of the tests indicate that the feasibility of life extension of many popular nuclear power plant cable products is promising and that mechanical measurements (primarily elongation, modulus, and density) were more effective than electrical measurements for monitoring age-related degradation. In the high temperature steam test, ethylene propylene rubber (EPR) cable materials generally survived to higher temperatures than crosslinked polyolefin (XLPO) cable materials. In dielectric testing after the submergence testing, the XLPO materials performed better than the EPR materials. This paper presents some recent experimental data that are not yet available elsewhere and a summary of findings from the entire experimental program.

Combustion of energetic materials involves processes in both gas and condensed phases and is governed by coupled thermal, physical, and chemical phenomena. Development of reliable models for design, performance, stability, and hazard analyses requires detailed understanding of three general chemical reaction regimes: (1) initial condensed-phase decomposition, (2) subsequent interaction of decomposition products with the remaining condensed phase, and (3) gas-phase reaction of decomposition products to form the ultimate combustion products. The first two regimes are the least understood and most difficult to study, particularly the initial condensed-phase decomposition. The basic difficulty in studying condensed phase phenomena has been the inability to probe directly chemistry in the condensed phase under isothermal condition and with the spatial and temporal resolution needed at higher temperatures and reaction rates. Thin-film samples provide a means to study condensed-phase chemistry at isothermal conditions and with microsecond temporal resolution. We are developing an experiment system that employs rapidly heated thin- film samples and multiple diagnostics to examine condensed-phase chemistry and monitor evolved gas species. Results from our initial work have been encouraging. Thin-film samples of several energetic materials have been prepared and appear to be representative of bulk materials. Furthermore, preliminary experiments indicate that all the use of these samples with two chemical diagnostic techniques, time-of- flight mass spectrometry (TOFMS) and time-resolved infrared spectral photography (TRISP), is viable. 5 refs., 8 figs.

The polarimetry problem (the measurement of the radar-cross-section polarization scattering matrix) is described. Two methods of calibrating a polarimetric radar are outlined. The first is a general multiple-calibration-target (MCT) method applicable to almost any radar system. The second is a simple, single-calibration-target (SCT) method applicable to systems which use a single antenna for both transmit/receive and a reciprocal RF network. The performance of the MCT method is examined through the use of Monte Carlo simulations. Finally, the SCT method is applied to measurements from the SCATTER facility, demonstrating about 40 dB isolation between polarization components in the frequency domain and in excess of 50 dB in the range domain.

Sandia National Laboratories, Albuquerque, manages the Utility Battery Exploratory Technology Development Program, which is sponsored by the US Department of Energy`s Office of Energy Management. In this capacity, Sandia is responsible for the engineering analyses and development of rechargeable batteries for utility-energy-storage applications. This report details the technical achievements realized during fiscal year 1991. Sodium/Sulfur, Zinc/Bromine, Nickel/Hydrogen, Aluminium/Air and Lead/Acid batteries are evaluated.

This report presents the results of single-well hydraulic tests performed in seven wells in the vicinity of the Waste Isolation Pilot Plant (WIPP) site in southeastern New Mexico in 1988 and 1989. The tests discussed in this report were performed in four of the five members of the Rustler Formation. The tests include: a slug-withdrawal test of the unnamed lower member of the Rustler ate well H-16; slug-withdrawal and slug-injection tests of the Culebra Dolomite Member at well AEC-7; slug-injection tests of the Culebra at well D-268; a pumping test of the Culebra at well H-18; pulse-injection, slug-injection, and slug-withdrawal tests of Magenta Dolomite Member at well H-2b1; pulse-withdrawal, slug-withdrawal, and slug-injection tests of the Magenta at well H-3b1; and pulse-withdrawal and slug-withdrawal tests of the Forty-niner Member at well H-3d. The tests were intended to provide data on the transmissivities of the Rustler members for use in regional-scale modeling of groundwater flow through the Rustler.

The ES&H Training Catalog is a tool to assist managers in determining which training courses they require their employees to complete. The narrative description under ``Who Shall Attend`` describes the characteristics of the employees and contractors under the direction of Sandia who are required by law, regulation, DOE Order, or SNL Directive to complete the training in order to be in compliance. The narrative is ``Who Should Attend`` describes the individuals for which the course is `highly recommended,`` although they are not mandated to attend.

SANET is a computer program intended for use in constructing, evaluating, and printing event trees for safety and reliability studies. SANET allows the user to graphically construct event trees, assign probabilities to the branches on the tree and include a variety of labels. Fast, publication quality printed output can be obtained. SANET requires an IBM compatible PC with a 80286, 80386 or 80486 processor, VGA graphics, a mouse and an HP Laserjet printer.

The technological innovation process comprises a range of stages, steps, and activities extending fro generation of new ideas through successful practical application of those ideas. This process constitutes the larger context within which technology transfer programs must necessarily operate if the goal is to transform new knowledge and technology into products which are competitive in the emerging global marketplace. A basic grasp of the principles and issues involved in the total innovation process is essential for developing and improving programs, prioritizing activities, and making strategic and operational decisions which will be appropriate and effective. This report attempts to provide a relatively brief overview of the total innovation process and related issues. This focus follows from the intent of the federal technology transfer initiative, which is essentially to facilitate the rapid advance of technological progress and to enhance national economic competitiveness. It is important to recognize that the transfer of technology between organizations or individuals is only one part of the complete process, although possibly a critical part in some situations. From an economic standpoint, technology transfer without resulting successful commercialization is meaningless. This report should be useful primer for people from any sector of the economy, even though it is intended to address the context for the transfer of federal laboratory technology in particular.

Two laboratory tests were designed to study the behavior of SnTe and CsOH in steam at {approximately}1230 K with the reactor materials Inconel 600, 304 stainless steel, silver and nickel, a nonoxidizing constituent of Inconels and steels in reactor accident environments. Thermochemical calculations examined the sensitivity of species in the H-O-Cs-Te-Sn system to temperature, to hydrogen and SnTe concentrations and to total system pressure. Test results indicate that should SnTe be formed from fission product tellurium and the tin in zircaloy cladding, it may not remain stable in steam in the presence of unoxidized (or lightly oxidized) metals. Calculations show a small amount of SnTe, in equilibrium with steam, decomposes to primarily Te and SnO. It is felt that these decomposition producets react with the reactor materials since Sn and Te are seen to be deposited separately and not as SnTe. The deposition velocity for SnO vapor in the system was estimated to be 0.57 m/s. The response of CsOH in the system was similar to behavior observed previously: some cesium combined with silicon found in Inconel and stainless steel oxides. At lower temperatures ({le}940 K) CsOH corroded Inconel, stainless steel and nickel providing enhanced surface area for additional trapping of species. There was no experimental evidence for the formation of cesium telluride; vapor equilibrium calculations did not predict formation either. 33 refs., 29 figs, 7 tabs.

An analytical expression and an integral representation are presented for the contribution, Y{sup (n)} = {Sigma}{sup (n)}{delta}{Omega}{sub n}, of n-event multiple scattering chains to the observed backscattering spectrum in heavy-ion backscattering (HIBS) measurements. The approximations introduced in deriving the results are chosen such that an upper limit is placed on Y{sup (n)} by the expressions. The Rutherford elastic scattering cross section is used to describe individual collisions between incident projectiles and target atoms. Screening of the Rutherford scattering cross section is included in an approximate fashion which maintains the upper limit estimate. Inelastic energy loss between collision events is assumed proportional to the projectile velocity. Specific application of theses expressions is made to HIBS detection of trace amounts of heavy atom impurities on a Si by 200--400 keV C{sup +} beams. The predicted multiple scattering background for this applications is compared with the predicted single scattering signal for 10{sup 10} atoms/cm{sup 2} of Fe, Cu, Zr, Sn, or Au, as surface impurity. The comparison shows that the multiple scattering background poses no barrier to extending the sensitivity of HIBS detection of impurities in this mass range to levels as low as 10{sup 8} atoms/cm{sup 2} for the upper part of the energy range considered. Comparison of calculations with and without screening included show that the screening of the Rutherford cross section by atomic electrons is a significant factor in preventing multiple scattering effects from interfering with HIBS spectrometry at impurity levels in the 10{sup 10} atoms/cm{sup 2} range.

The MELCOR code has been used to simulate the FLECHT SEASET natural circulation experiments done in a scale-model Westinghouse-PWR test facility, with code results compared to experimental data. Sensitivity studies have been done, for both single-phase and two-phase natural circulation conditions, on time step effects and machine dependencies; nodalization studies and studies on several code modelling options were also done. Good agreement is found between prediction and observation for steady-state, single-phase liquid natural circulation. The code could reproduce the major thermal/hydraulic response characteristics in two-phase natural circulation, but only through a number of nonstandard input modelling modifications; MELCOR cannot reproduce the requisite physical phenomena with ``normal`` input models. Because the same response is observed in similar tests at other facilities over a range of scales and is expected to occur in full-scale plants as well, the ability of the user to ``match`` the observed behavior through a small set of nonstandard input modelling changes allows MELCOR to be used in PRA studies in which such physics are expected to be encountered, while awaiting corrections to the code models involved. The time step control algorithm in MELCOR does not run this problem efficiently; a substantial reduction in time step results in significantly less oscillation predicted at only a small increase run time.

The Waste Isolation Pilot Plant (WIPP) is planned as a mined geologic repository for the disposal of transuranic (TRU) radioactive wastes generated by defense programs of the United States Department of Energy. One of the criteria for evaluating the suitability of the WIPP for disposal of TRU wastes is compliance with the United States Environmental Protection Agency`s (EPA) standards for such facilities. The Containment Requirements of those standards require calculating cumulative releases of radionuclides to the accessible environment by all combinations of events and processes (scenarios) that may affect the escape and transport of radionuclides from the repository for 10, 000 years after decommissioning of the facility. Because the release limits established by the EPA are probabilistic, scenario probabilities are also required. A panel of experts was convened to estimate the probabilities of occurrence of the events used in scenario development and to identify additional human-intrusion events for inclusion in a safety assessment of the WIPP. This report documents the background presentations that were made to the panel about the WIPP program, regulatory guidelines, and performance-assessment program, and site-specific and regional geologic and hydrologic characteristics that may affect the WIPP disposal system.

This report is a revision and update of the original geologic site characterization report that was published in 1980. Many of the topics addressed in the earlier report were predictive in nature and it is now possible to reexamine them some 12 years later, using the data from 17 new caverns and more than ten years of SPR storage experience. Revised maps of the salt configuration show on overhand and faults on the north side of the dome, defining more clearly the edge relationships with respect to the SPR caverns. Caprock faults may locally influence the pattern of subsidence, which is occurring primarily as a result of cavern creep closure. The greater subsidence rate occurring at West Hackberry will likely require mitigative action within a few years. Seismicity of low intensity recurs infrequently at West Hackberry, but a small earthquake in 1983 caused dish rattling in the immediate vicinity.

Excavation stability in an underground nuclear waste repository is required during construction, emplacement, retrieval (if required), and closure phases to ensure worker health and safety, and to prevent development of potential pathways for radionuclide migration in the post-closure period. Stable excavations are developed by appropriate excavation procedures, design of the room shape, design and installation of rock support reinforcement systems, and implementation of appropriate monitoring and maintenance programs. In addition to the loads imposed by the in situ stress field, the repository drifts will be impacted by thermal loads developed after waste emplacement and, periodically, by seismic loads from naturally occurring earthquakes and underground nuclear events. A priori evaluation of stability is required for design of the ground support system, to confirm that the thermal loads are reasonable, and to support the license application process. In this report, a design methodology for assessing drift stability is presented. This is based on site conditions, together with empirical and analytical methods. Analytical numerical methods are emphasized at this time because empirical data are unavailable for excavations in welded tuff either at elevated temperatures or under seismic loads. The analytical methodology incorporates analysis of rock masses that are systematically jointed, randomly jointed, and sparsely jointed. In situ thermal and seismic loads are considered. Methods of evaluating the analytical results and estimating ground support requirements for all the full range of expected ground conditions are outlines. The results of a preliminary application of the methodology using the limited available data are presented. 26 figs., 55 tabs.

The Yucca Mountain Site Characterization Project is studying Yucca Mountain in southwestern Nevada as a potential site for a high-level nuclear waste repository. Site characterization includes surface-based and underground testing. Analyses have been performed to design site characterization activities with minimal impact on the ability of the site to isolate waste, and on tests performed as part of the characterization process. One activity of site characterization is the construction of an Exploratory Studies Facility, for which many design options are being considered, including shafts, drifts, and ramps. The information in this report pertains to: (1) engineering calculations of the potential distribution of residual water from constructing the shafts and drifts; (2) numerical calculations predicting the movement of residual construction water from the shaft and drift walls into the rock; and (3) numerical calculations of the movement of residual water and how the movement is affected by ventilation. This document contains information that has been used in preparing Appendix 1 of the Exploratory Studies Facility Design Requirements document for the Yucca Mountain Project.

Numerical results are presented for the Performance Assessment Calculational Exercise (PACE-90). One- and two-dimensional water and solute transport are presented for steady infiltration into Yucca Mountain. Evenly distributed infiltration rates of 0.01, 0.1, and 0.5 mm/yr were considered. The calculations of solute transport show that significant amounts of radionuclides can reach the water table over 100,000 yr at the 0.5 mm/yr rate. For time periods less than 10,000 yr or infiltrations less than 0.1 mm/yr very little solute reaches the water table. The numerical simulations clearly demonstrate that multi-dimensional effects can result in significant decreases in the travel time of solute through the modeled domain. Dual continuum effects are shown to be negligible for the low steady state fluxes considered. However, material heterogeneities may cause local amplification of the flux level in multi-dimensional flows. These higher flux levels may then require modeling of a dual continuum porous medium.

In support of the development of American National Standards Institute standards for the transport of radioactive materials, Sandia has a program to characterize the normal transport environment. This program includes both analytical modeling of package and trailer responses, and over-the-road tests to measure those responses. This paper presents the results of a series of over-the-road tests performed using Chem-Nuclear equipment in the Barnwell, SC, area. The test events included a variety of road types such as rough concrete, shock events such as railroad grade crossings, and driver responses such as sharp turns. The response of the package and trailer to these events was measured with accelerometers at various locations to determine the inertial loads. Either load cells or strain gages were used to measure tiedown response. These accelerations and loads were measured on systems with flexible and ``rigid`` tiedowns. The results indicated that while significant accelerations occur on the trailer bed, these do not translate into equivalent loads in either the package or the tiedown system. This indicates that trailer-bed response should not be used in determining the load factor for fatigue calculations of the package components or in determining design loads for tiedowns.

Currently, there are several Greater Confinement Disposal (GCD) boreholes at the Radioactive Waste Management Site (RWMS) for the Nevada Test Site. These are intermediate-depth boreholes used for the disposal of special case wastes, that is, radioactive waste within the Department of Energy complex that do not meet the criteria established for disposal of high-level waste, transuranic waste, or low-level waste. A performance assessment is needed to evaluate the safety of the GCD site, and to examine the feasibility of the GCD disposal concept as a disposal solution for special case wastes in general. This report documents the effort in defining all the waste inventory presently disposed of at the GCD site, and the inventory and release model to be used in a performance assessment for compliance with the Environmental Protection Agency`s 40 CFR 191.

Cementitious materials, together with other materials, are being considered to seal a potential repository at Yucca Mountain. A concern with cementitious materials is the chemical and mineralogic changes that may occur as these materials age while in contact with local ground waters. A combined theoretical and experimental approach was taken to determine the ability to theoretically predict mineralogic changes. The cementitious material selected for study has a relatively low Ca:Si ratio approaching that of the mineral tobermorite. Samples were treated hydrothermally at 200{degrees}C with water similar to that obtained from the J-13 well on the Nevada Test Site. Post-test solutions were analyzed for pH as well as dissolved K, Na, Ca, Al, and Si. Solid phases formed during these experiments were characterized by scanning electron microscopy and X- ray diffraction. These findings were compared with predictions made by the geochemical modeling code EQ3NR/E06. It was generally found that there was good agreement between predicted and experimental results.

Two years ago, researchers at Sandia National Laboratories showed that a massively parallel computer with 1024 processors could solve scientific problems more than 1000 times faster than a single processor. Since then, interest in massively parallel processing has increased dramatically. This review paper discusses some of the applications of this emerging technology to important problems at Sandia. Particular attention is given here to the impact of massively parallel systems on applications related to national defense. New concepts in heterogenous programming and load balancing for MIMD computers are drastically increasing synthetic aperture radar (SAR) and SDI modeling capabilities. Also, researchers are showing that the current generation of massively parallel MIMD and SIMD computers are highly competitive with a CRAY on hydrodynamic and structural mechanics codes that are optimized for vector processors.

Analysis of an intermittent failure to write the "1" state to a particular memory location at low temperature (-55° C) in a 16K x 1 CMOS SRAM is presented. The failure was found to be due to an open metallization at a metal-to-silicon contact. The root cause of the failure was poor step coverage of the metallization over an oxide step. A variety of failure analysis techniques including dynamic electron beam analysis at low temperature using a Peltier cold stage were employed to study the intermittently failing SRAM. The failure site was located by using capacitive coupling voltage contrast analysis. PSPICE simulation, light emission microscopy, scanning electron microscopy, and focused-ion beam techniques were used to confirm the failure mechanism and location. The write cycle time of the failed IC was abnormally long, but within the allowable tester limit. The vulnerability of other ICs to failure by open metallization in metal-to silicon contacts is reviewed.

We present calculations which show the radial dependence of the KVV and L₂₃VV Auger matrix elements of silicon. We find greatly differing dependences, converging within ~1 a.u. of the nucleus in the case of the KVV, but not until ~4 a.u. in the case of the L₂₃VV, well beyond the bond midpoint of ~2.2 a.u. We also find quite different dependences for the various elements within a particular CVV transition. Because the local density of states (LDOS) is dependent on the radius of the sphere of integration, our results suggest that different CVV Auger processes on the same atom in fact probe different LDOSs, as do even different contributions within the same transition. (This effect is separate from the well-known matrix element property which weights angular-momentum components differently.) These results call into question both the single-site LDOS approximation when used in the interpretation of low-energy (<100 eV) Auger spectra, and the application to high-energy spectra of local densities of states obtained by integration over muffin-tin or Wigner–Seitz spheres which have a large radius compared to the region probed by the Auger process.

Two-dimensional Acousto-Optic (AO) correlators differ from the frequency plane correlators in that multiplying, shifting, and adding, rather than Fourier transforming are used to obtain the correlations. Thus, many of the available composite filter design techniques are not aimed at designing filters for use in AO correlators since they yield frequency-domain functions. In this paper, a method is introduced for designing filter impulse responses of arbitrary extent for implementation on AO correlators. These filters are designed to yield sharp correlation peaks. Simulation results are included to illustrate the viability of the proposed approach. Also included are some initial results from the first successful use of grey-scale composite filters on an AO correlator.