Publications

Sub-optimal design practices can reduce the radiation hardness of a circuit even though it is fabricated in a radiation hardened process. This is especially true for a nonvolatile memory, as compared to a standard digital circuit, where high voltages and unusual bias conditions are required. This paper will discuss the design techniques used in the development of a 64K EEPROM (Electrically Erasable Programmable Read Only Memory) to maximize radiation hardness. The circuit radiation test results will be reviewed in order to provide validation of the techniques. © 1993 IEEE

The U.S. Department of Energy (DOE) is developing the Waste Isolation Pilot Plant (WIPP) in southeastern New Mexico as a facility for the long-term disposal of defense-related transuranic (TRU) wastes. Use of the WIPP for waste disposal is contingent on demonstrations of compliance with applicable regulations of the U.S. Environmental Protection Agency (EPA). This paper addresses issues related to modeling gas and brine migration at the WIPP for compliance with both EPA 40 CFR 191 (the Standard) and 40 CFR 268.6 (the RCRA). At the request of the WIPP Project Integration Office (WPIO) of the DOE, the WIPP Performance Assessment (PA) Department of Sandia National Laboratories (SNL) has completed preliminary uncertainty and sensitivity analyses of gas and brine migration away from the undisturbed repository. This paper contains descriptions of the numerical model and simulations, including model geometries and parameter values, and a summary of major conclusions from sensitivity analyses. Because significant transport of contaminants can only occur in a fluid (gas or brine) medium, two-phase flow modeling can provide an estimate of the distance to which contaminants can migrate. Migration of gas or brine beyond the RCRA 'disposal-unit boundary' or the Standard's accessible environment constitutes a potential, but not certain, violation and may require additional evaluations of contaminant concentrations.

This paper presents an infinite impulse response (IIR) filtering technique for reducing structural vibration in remotely operated robotic systems. The technique uses a discrete filter between the operator's joy stick and the robot controller to alter the inputs of the system so that residual vibration and swing are reduced. A linearized plant model of the system is analyzed in the discrete time domain, and the filter is designed using pole-zero placement in the z-plane. This technique has been successfully applied to a two link flexible arm and a gantry crane with a suspended payload.

Before disposing of transuranic radioactive waste at the Waste Isolation Pilot Plant (WIPP), the United States Department of Energy (DOE) must evaluate compliance with long-term regulations of the United States Environmental Protection Agency (EPA). Sandia National Laboratories (SNL) is conducting iterative performance assessments (PAs) of the WIPP for the DOE to provide interim guidance while preparing for final compliance evaluations. This paper describes the 1992 preliminary comparison with Subpart B of the Environmental Standards for the Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Radioactive Wastes (40 CFR 191), which regulates long-term releases of radioactive waste. Results of the 1992 PA are preliminary, and cannot be used to determine compliance or noncompliance with EPA regulations because portions of the modeling system and data base are incomplete. Results are consistent, however, with those of previous iterations of PA, and the SNL WIPP PA Department has high confidence that compliance with 40 CFR 191B can be demonstrated. Comparison of predicted radiation doses from the disposal system also gives high confidence that the disposal system is safe for long-term isolation.

Before disposing of transuranic radioactive waste at the Waste Isolation Pilot Plant (WIPP), the United States Department of Energy (DOE) must evaluate compliance with long-term regulations of the United States Environmental Protection Agency (EPA), specifically the Environmental Standards for the Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Radioactive Wastes (40 CFR 191), and the Land Disposal Restrictions (40 CFR 268) of the Hazardous and Solid Waste Amendments to the Resource Conservation and Recovery Act (RCRA). Sandia National Laboratories (SNL) is conducting iterative performance assessments (PAs) of the WIPP for the DOE to provide interim guidance while preparing for final compliance evaluations. This paper provides background information on the regulations, describes the SNL WIPP PA Department's approach to developing a defensible technical basis for consistent compliance evaluations, and summarizes the major observations and conclusions drawn from the 1991 and 1992 PAs.

Two space applications are considered for (electrically) contactless coilguns: launch of small satellites into low-earth orbit, and launch of lunar liquid oxygen (LLOX) from the moon to the stationary Lagrangian point L2. For the earth-to-orbit (ETO) application, the baseline conceptual design consists of a 960-m long gun sited in a tunnel at 25° inclination. The gun launches an 1820-kg package that includes a 100-kg satellite and a 650-kg boost rocket for orbital insertion. For the lunar application, the launcher is 200 m long. A 100-kg load of LLOX is packaged in a 10-kg fiber-wrapped tank, accelerated at 2 kgees in an aluminum bucket (armature), and launched at 2.33 km/s at 30-minute intervals. The canisters arrive at L2 2.97 days later and are captured by robotic tugs that deliver them to a fuel depot. The total mass of LLOX delivered per year is 867 Mg (metric tons). © 1993 IEEE

We postulate that there is a required velocity of magnetic flux into a plasma-armature railgun to keep the armature from expanding toward the breech and thereby changing the nature of the forward-going force. Testing this idea with data from several guns we find confirmation, provided the required flux velocity increases with the square of the linear current density. This leads to the conclusion that these guns have a momentum limit while they are operating as true electromagnetic launchers. That limit is derived. High-performance railguns reach this limit during current rise, and most of their momentum comes from a hybrid phase. The hybrid phase has a pressure limitation at high velocity similar to electrothermal, light-gas, and powder guns. The 6 km/s “velocity limit” is explained by this pressure limitation. © 1993 IEEE

Different hardness-assurance tests are often required for advanced bipolar devices than for CMOS devices. In this work, the dose-rate dependence of bipolar current-gain degradation is mapped over a wide range of dose rates for the first time, and it is very different from analogous MOSFET curves. Annealing experiments following irradiation show negligible change in base current at room temperature, but significant recovery at temperatures of 100°C and above. In contrast to what is observed in MOSFET’s, irradiation and annealing tests cannot be used to predict the low-dose-rate response of bipolar devices. A comparison of x-ray-induced and 60Co gamma-ray-induced gain degradation is reported for the first time for bipolar transistors. The role of the emitter bias during irradiation is also examined. Implications for hardening and hardness assurance are discussed. © 1993 IEEE

Vadose-zone moisture transport near an impermeable barrier has been under study at a field site near Albuquerque, NM since 1990. Moisture content and temperature have been monitored in the subsurface on a regular basis; both undergo a seasonal variation about average values. Even though the slab introduces two-dimensional effects on the scale of the slab, moisture and heat transport is predominantly vertical. Numerical simulations, based on the models developed by Philip and de Vries (1957) and de Vries (1958), indicate that the heat flow is conduction-dominated while the moisture movement is dominated by diffusive vapor distillation. Model predictions of the magnitude and extent of changes in moisture content underneath the slab are in reasonable agreement with observation.

A Sandia National Laboratories/AT&T Bell Laboratories Team is developing a soft x-ray projection lithography tool that uses a compact laser plasma as a source of 14 nm x-rays. Optimization of the 14 nm x-rays source brightness is a key issue in this research. This paper describes our understanding of the source as it has been obtained through the use of computer simulations utilizing the LASNEX radiation-hydrodynamics code.

Ion-beam-induced charge-collection imaging (IBICC) has been used to study the SEU mechanisms of the Sandia TA670 16K-bit SRAM. Quantitative charge-collection spectra from known regions of the memory cell have been derived with this technique. For 2.4-MeV He ions at normal incidence, charge collection depth for a reverse-biased p+ drain strike is estimated to be 4.8±0.4 μm. Heavy-ion strikes to the reverse-biased p-well result in nearly complete collection of deposited charge to a depth of 5.5±0.5 μm. A charge amplification effect in the n-on drain is identified and is due to either bipolar amplification or a shunt effect in the parasitic vertical npn bipolar transistor associated with the n+/n substrate, p-well, and n+ drain. This effect is present only when the n+ drain is at 0V bias. When coupled with previous SEU-imaging, these results strongly suggest that the dominant SEU mechanism in this SRAM is a heavy-ion strike to the n-on transistor drain. © 1993 IEEE

The Vapnik-Chervonenkis (V-C) dimension is an important combinatorial tool in the analysis of learning problems in the PAC framework. For polynomial learnability, we seek upper bounds on the V-C dimension that are polynomial in the syntactic complexity of concepts. Such upper bounds are automatic for discrete concept classes, but hitherto little has been known about what general conditions guarantee polynomial bounds on V-C dimension for classes in which concepts and examples are represented by tuples of real numbers. In this paper, we show that for two general kinds of concept class the V-C dimension is polynomially bounded as a function of the syntactic complexity of concepts. One is classes where the criterion for membership of an instance in a concept can be expressed as a formula (in the first-order theory of the reals) with fixed quantification depth and exponentially-bounded length, whose atomic predicates are polynomial inequalities of exponentially-bounded degree. The other is classes where containment of an instance in a concept is testable in polynomial time, assuming we may compute standard arithmetic operations on reals exactly in constant time. Our results show that in the continuous case, as in the discrete, the real barrier to efficient learning in the Occam sense is complexity-theoretic and not information-theoretic. We present examples to show how these results apply to concept classes defined by geometrical figures and neural nets, and derive polynomial bounds on the V-C dimension for these classes.

During hydrocarbon reservoir stimulations, such as hydraulic fracturing, the cracking and slippage of the formation results in the emission of seismic energy. The objective of this study was to determine the properties of these induced micro-seisms. A hydraulic fracture experiment was performed in the Piceance Basin of Western Colorado to induce and record micro-seismic events. The formation was subjected to four processes; breakdown/ballout, step-rate test, KCL mini-fracture, and linear-gel mini-fracture. Micro-seisms were acquired with an advanced three-component wall-locked seismic accelerometer package, placed in an observation well 211 ft offset from the fracture well. During the two hours of formation treatment, more than 1200 micro-seisms with signal-to-noise ratios in excess of 20 dB were observed. The observed micro- seisms had a nominally flat frequency spectrum from 100 Hz to 1500 Hz and lack the spurious tool-resonance effects evident in previous attempts to measure micro-seisms. Both p-wave and s-wave arrivals are clearly evident in the data set, and hodogram analysis yielded coherent estimates of the event locations. This paper describes the characteristics of the observed micro- seismic events (event occurrence, signal-to-noise ratios, and bandwidth) and illustrates that the new acquisition approach results in enhanced detectability and event location resolution.

Many designs use EPLDs (Erasable Programmable Logic Devices) to implement control logic and state machines. If the design is slow, timing through the EPLD is not crucial so designers often treat the device as a black box. In high speed designs, timing through the EPLD is critical. In these cases a thorough understanding of the device architecture is necessary. Lessons learned in the implementation of a high-speed design using the Altera EPM5130 are discussed.

A new Assembly Test Chip, ATC04, designed to measure mechanical stresses at the die surface has been built and tested. This CMOS chip 0.25 in. on a side, has an array of 25 piezoresistive stress sensing cells, four resistive heaters and two ring oscillators. The ATCO4 chip facilitates making stress measurements with relatively simple test equipment and data analysis. The design, use, and accuracy of the chip are discussed and initial results are presented from three types of stress measurement experiments: four-point bending calibration, single point bending of a substrate with an ATC04 attached by epoxy, and stress produced by a liquid epoxy encapsulant.

The feasibility of utilizing a groundbased laser without an orbital mirror for space debris removal is examined. Technical issues include atmospheric transmission losses, adaptive-optics corrections of wavefront distortions, laser field-of-view limitations, and laser-induced impulse generation. The physical constraints require a laser with megawatt output, long run-time capability, and wavelength with good atmospheric transmission characteristics. It is found that a 5-MW reactor-pumped laser can deorbit debris having masses of the order of one kilogram from orbital altitudes to be used by Space Station Freedom. Debris under one kilogram can be deorbited after one pass over the laser site, while larger debris can be deorbited or transferred to alternate orbits after multiple passes over the site.

Flux-flow based devices such as the superconducting flux flow transistor and magnetically controlled long junctions have been made from thin films of TICaBaCuO and YBaCuO. The devices are based on the magnetic control of flux flow in their respective structures: a long junction or an array of weak links. The equivalent circuits of the two devices are similar: a low impedance input control line, an output impedance of 3-20 Q and an active current-controlled element. The long junctions have tended to be slower, have lower gain and be somewhat less noisy than their counterparts. The performance of circuits such as narrowband and distributed amplifiers (SO GHz bandwidths, noise figures < 3 dB), phase shifters (continuous with < 2 dB loss 4-40 GHz), logic gates (2-3 ps gate delays) and memories made using these devices will be compared and analyzed. © 1993 IEEE

Automatic motion planning of a spray cleaning robot with collision avoidance is presented in this paper. In manufacturing environments, electronic and mechanical components are traditionally cleaned by spraying or dipping them using chlorofluorocarbon (CFC) solvents. As new scientific data show that such solvents are major causes for stratospheric ozone depletion, an alternate cleaning method is needed. Part cleaning with aqueous solvents is environmentally safe, but can require precision spraying at high pressures for extended time periods. Operator fatigue during manual spraying can decrease the quality of the cleaning process. By spraying with a robotic manipulator, the necessary spray accuracy and consistency to manufacture high-reliability components can be obtained. Our motion planner was developed to automatically generate motions for spraying robots based on the part geometry and cleaning process parameters. For spraying paint and other coatings a geometric description of the parts and robot may be sufficient for motion planning, since coatings are usually done over the visible surfaces. For spray cleaning, the requirement to reach hidden surfaces necessitates the addition of a rule-based method to the geometric motion planning.

The geochemical properties of a porous sand and several tracers (Ni, Br, and Li) have been characterized for use in a caisson experiment designed to validate sorption models used in models of reactive transport. The surfaces of the sand grains have been examined by a combination of techniques including potentiometric titration, acid leaching, optical microscopy, and scanning electron microscopy with energy-dispersive spectroscopy. The surface studies indicate the presence of small amounts of carbonate, kaolinite and iron-oxyhydroxides. Adsorption of nickel, lithium and bromide by the sand was measured using batch techniques. Bromide was not sorbed by the sand. A linear (Kd) or an isotherm sorption model may adequately describe transport of Li; however, a model describing the changes of pH and the concentrations of other solution species as a function of time and position within the caisson and the concomitant effects on Ni sorption may be required for accurate predictions of nickel transport.

For problems where media properties are measured at one scale and applied at another, scaling laws or models must be used in order to define effective properties at the scale of interest. The accuracy of such models will play a critical role in predicting flow and transport through the Yucca Mountain Test Site given the sensitivity of these calculations to the input property fields. Therefore, a research program has been established to gain a fundamental understanding of how properties scale with the aim of developing and testing models that describe scaling behavior in a quantitative manner. Scaling of constitutive rock properties is investigated through physical experimentation involving the collection of suites of gas permeability data measured over a range of discrete scales. Also, various physical characteristics of property heterogeneity and the means by which the heterogeneity is measured and described and systematically investigated to evaluate their influence on scaling behavior. This paper summarizes the approach that is being taken toward this goal and presents the results of a scoping study that was conducted to evaluate the feasibility of the proposed research.

Experimental results exploring gravity-driven wetting front instability in a pre-wetted, rough-walled analog fracture are presented. Initial conditions considered include a uniform moisture field wetted to field capacity of the analog fracture and the structured moisture field created by unstable infiltration into an initially dry fracture. As in previous studies performed under dry initial conditions, instability was found to result both at the cessation of stable infiltration and at flux lower than the fracture capacity under gravitational driving force. Individual fingers were faster, narrower, longer, and more numerous than observed under dry initial conditions. Wetting fronts were found to follow existing wetted structure, providing a mechanism for rapid recharge and transport.

In an attempt to achieve completeness and consistency, the performance-assessment analyses developed by the Yucca Mountain Project are tied to scenarios described in event trees. Development of scenarios requires describing the constituent features, events, and processes in detail. Several features and processes occurring at the waste packages and the rock immediately surrounding the packages (i.e., the near field) have been identified: the effects of radiation on fluids in the near-field rock, the path-dependency of rock-water interactions, and the partitioning of contaminant transport between colloids and solutes. This paper discusses some questions regarding these processes that the near-field performance-assessment modelers will need to have answered to specify those portions of scenarios dealing with the near field.

Experiments investigating the behavior of individual, gravity-driven fingers in an initially dry, rough-walled analog fracture are presented. Fingers were initiated from constant flow to a point source. Finger structure is described in detail; specific phenomena observed include: desaturation behind the finger-tip, variation in finger path, intermittent flow structures, finger-tip bifurcation, and formation of dendritic sub-fingers. Measurements were made of finger-tip velocity, finger width, and finger-tip length. Non-dimensional forms of the measured variables are analyzed relative to the independent parameters, flow rate and gravitational gradient.

The potential repository system is intended to isolate high-level radioactive waste at Yucca Mountain according to the performance objective - 10 CFR 60.112. One subsystem that may contribute to achieving this objective is the sealing subsystem. This subsystem is comprised of sealing components in the shafts, ramps, underground network of drifts, and the exploratory boreholes. Sealing components can be rigid, as in the case of a shaft seal, or can be more compressible, as in the case of drift fill comprised of mined rockfill. This paper presents the preliminary seismic response of discrete sealing components in welded and nonwelded tuff. Special consideration is given to evaluating the stress in the seal, and the behavior of the interface between the seal and the rock. The seismic responses are computed using both static and dynamic analyses. Also presented is an evaluation of the maximum seismic response encountered by a drift seal with respect to the angle of incidence of the seismic wave. Mitigation strategies and seismic design considerations are proposed which can potentially enhance the overall response of the sealing component and subsequently, the performance of the overall repository system.

This analytical study examines the effect of initial thermodynamic conditions on the loads generated by the combustion of homogeneous hydrogen-air-steam mixtures. The effect of initial temperature, pressure, hydrogen concentration, and steam concentration is evaluated for two cases, (1) constant volume and (2) constant initial pressure. For each case, the Adiabatic, Isochoric, Complete Combustion (AICC), Chapman-Jouguet (CJ), and normally reflected CJ pressures are calculated for a range of hydrogen and steam concentrations representative of the entire flammable regime. For detonation loads, pressure profiles and time-histories are also evaluated in one-dimensional Cartesian geometry. The results show that to a first approximation, the AICC and CJ pressures are directly proportional to the initial density. Increasing the hydrogen concentration up to stoichiometric concentrations significantly increases the AICC, CJ, and reflected CJ pressures. For the constant volume case, the AICC, CJ, and reflected CJ pressures increase with increasing hydrogen concentration on the rich side of stoichiometric concentrations. For the constant initial pressure case, the AICC, CJ, and reflected CJ pressures decrease with increasing hydrogen concentration on the rich side of stoichiometric values. The addition of steam decreases the AICC, CJ, and reflected CJ pressures for the constant initial pressure case, but increases them for the constant volume case. For detonations, the pressure time-histories can be normalized with the AICC pressure and the reverberation time for Cartesian geometry. © 1993.

We present a new methodology for predicting the fatigue life of solder joints for electronics applications. This approach involves the integration of experimental and computational techniques. The first stage involves correlating the manufacturing and processing parameters with the starting microstructure of the solder joint. The second stage involves a series of experiments that characterize the evolution of the microstructure during thermal cycling. The third stage consists of a computer modeling and simulation effort that utilizes the starting microstructure and experimental data to produce a reliability prediction of the solder joint. This approach is an improvement over current methodologies because it incorporates the microstructure and properties of the solder directly into the model and allows these properties to evolve as the microstructure changes during fatigue.

The results of experiments on the energy losses in switches are presented. All experiments were circuit code modeled after the development of a new switch loss version based on Braginskii's spark channel theory for gasses. The circuit code predicts switch energy loss and peak currents as a function of time. Slight constant offsets between the theory and data that depended on the switching dielectric were observed. Plasma conductivity for each tested dielectric was modified to lessen this offset. An excellent agreement between the theory and experiment was also observed. The new model successfully predicted the series resistance of Marx generators.

Yucca Mountain, Nevada is a potential site for a high-level radioactive-waste repository. Uncertainty and sensitivity analyses were performed to estimate critical factors in the performance of the site with respect to a criterion in terms of pre-waste-emplacement ground-water travel time. The degree of failure in the analytical model to meet the criterion is sensitive to the estimate of fracture porosity in the upper welded unit of the problem domain. Fracture porosity is derived from a number of more fundamental measurements including fracture frequency, fracture orientation, and the moisture-retention characteristic inferred for the fracture domain.

This report contains an outline of the Solar Thermal Design Assistance Center's (STDAC) major activities and accomplishments in Fiscal Year 1992 (FY92). The report describes the resources allocated to fund STDAC and the personnel needed to carry out STDAC activities and accomplishments. It also contains a comprehensive list of persons that called STDAC for consultation in FY92.

Images taken with a synthetic aperture radar (SAR) on an airplane were distorted with phase errors generated by a computer program that simulates the propagation of radar waves through the disturbed ionosphere. The simulation is for an orbiting SAR imaging a scene on the ground. Both the spatially-invariant (decorrelation length projected onto the ground much larger than the scene size) and spatially-variant (decorrelation length much smaller than the scene size) cases are described. The spatially-invariant phase errors can be removed using several different algorithms. Problems and strategies in restoring SAR images distorted with spatially-variant phase errors are discussed.

Future challenges facing the nonproliferation community will undoubtedly change the normal way of doing business'' in international safeguards. New technology will emerge in support of compliance concepts such as transparency and openness, regional security assurance, bilateral cooperation, and special. or non-routine inspections. Technologies address in remote unattended monitoring, integrated on-site monitoring, environmental monitoring, satellite and aerial over-flight systems, equipment for special inspectios, and sharable data information fusion and management, are just a few examples of potential technologics for new nonproliferation monitoring regimes.

We describe an algorithm for the static load balancing of scientific computations that generalizes and improves upon spectral bisection. Through a novel use of multiple eigenvectors, our new spectral algorithm can divide a computation into 4 or 8 pieces at once. This leads to balanced partitions that have lower communication overhead and are less expensive to compute than those of spectral bisection. In addition, our approach automatically works to minimize message contention on a hypercube or mesh architecture.

This paper describes a collaborative effort between Sandia National Laboratories and the Rocketdyne Division of Rockwell International Corporation to develop an automated braze paste dispensing system for rocket engine nozzle manufacturing. The motivation for automating this manufacturing process is to reduce the amount of labor and excess material required. A critical requirement for this system is the automatic location of key nozzle features using non-contact sensors. Sandia has demonstrated that the low-cost Multi-Axis Seam Tracking (MAST) capacitive sensor can be used to accurately locate the nozzle surface and tube gaps.

This paper examines some of the metrics that are commonly used to design correlation filter's for optical pattern recognition, including: the Fisher ratio, the signal-to-noise ratio, the equal correlation peak (ECP) constraint, and normalized correlation. Attention is given to the underlying assumptions that are required to move from Bayesian decision theory to a particular metric or design principle. Since a Bayes classifier is statistically optimum, this provides a means for assessing the merit of a particular approach. Although the authors only examine a few metrics in this paper, the approach is general and should be useful for assessing the merit and applicability of any of the numerous filter designs that have been proposed in the optical pattern recognition community.

This report discusses electronic isolators which are used to maintain electrical separation between safety and non-safety systems in nuclear power plants. The concern is that these devices may fail allowing unwanted signals or energy to act upon safety systems, or preventing desired signals from performing their intended function. While operational history shows many isolation device problems requiring adjustments and maintenance, we could not find incidents where there was a safety implication. Even hypothesizing multiple simultaneous failures did not lead to significant contributions to core damage frequency. Although the analyses performed in this study were not extensive or detailed, there seems to be no evidence to suspect that isolation device failure is an issue which should be studied further.

This report documents the fiscal year 1992 activities of the, Utility Battery Storage Systems Program (UBS) of the US Department of Energy (DOE), Office of Energy Management (OEM). The UBS program is conducted by Sandia National Laboratories (SNL). UBS is responsible for the engineering development of integrated battery systems for use in utility-energy-storage (UES) and other stationary applications. Development is accomplished primarily through cost-shared contracts with industrial organizations. An important part of the development process is the identification, analysis, and characterization of attractive UES applications. UBS is organized into five projects: Utility Battery Systems Analyses; Battery Systems Engineering; Zinc/Bromine; Sodium/Sulfur; Supplemental Evaluations and Field Tests. The results of the Utility Systems Analyses are used to identify several utility-based applications for which battery storage can effectively solve existing problems. The results will also specify the engineering requirements for widespread applications and motivate and define needed field evaluations of full-size battery systems.

The sixth experiment of the Integral Effects Test (IET-6) series was conducted to investigate the effects of high pressure melt ejection on direct containment heating. Scale models of the Zion reactor pressure vessel (RPV), cavity, instrument tunnel, and subcompartment structures were constructed in the Surtsey Test Facility at Sandia National Laboratories. The RPV was modeled with a melt generator that consisted of a steel pressure barrier, a cast MgO crucible, and a thin steel inner liner. The melt generator/crucible had a hemispherical bottom head containing a graphite limitor plate with a 4-cm exit hole to simulate the ablated hole in the RPV bottom head that would be formed by ejection of an instrument guide tube in a severe nuclear power plant accident. The cavity contained 3.48 kg of water, which corresponds to condensate levels in the Zion plant, and the containment basement floor was dry. A 43-kg initial charge of iron oxide/aluminum/chromium thermite was used to simulate corium debris on the bottom head of the RPV. Molten thermite was ejected by steam at an initial pressure of 6.3 MPa into the reactor cavity. The Surtsey vessel atmosphere contained pre-existing hydrogen to represent partial oxidation of the zirconium in the Zion core. The initial composition of the vessel atmosphere was 87.1 mol.% N[sub 2], 9.79 mol.% O[sub 2], and 2.59 mol.% H[sub 2], and the initial absolute pressure was 198 kPa. A partial hydrogen burn occurred in the Surtsey vessel. The peak vessel pressure increase was 279 kPa in IET-6, compared to 246 kPa in the IET-3 test. The total debris mass ejected into the Surtsey vessel in IET-6 was 42.5 kg. The gas grab sample analysis indicated that there were 180 g[center dot] moles of pre-existing hydrogen, and that 308[center dot]moles of hydrogen were produced by steam/metal reactions. About 335 g[center dot]moles of hydrogen burned, and 153 g[center dot]moles remained unreacted.

The Sandia National Laboratories (SNL) Engineering Analysis Code Access System (SEACAS) is a collection of structural and thermal codes and utilities used by analysts at SNL. The system includes pre- and post-processing codes, analysis codes, database translation codes, support libraries, UNIX{trademark} shell scripts, 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 180 days of Cray Y-MP{trademark} CPU time have been used at SNL by SEACAS codes. The job mix includes jobs using only a few seconds of CPU time, up to jobs using two and one-half days of CPU time. SEACAS is running on several different systems at SNL including Cray Unicos, Hewlett Packard HP-UX{trademark}, Digital Equipment Ultrix{trademark}, and Sun SunOS{trademark}. This document is a short description of the codes the SEACAS system.

This bulletin from Sandia Laboratories presents current research on testing technology. Fiber optics systems at the Nevada Test Site is replacing coaxial cables. The hypervelocity launcher is being used to test orbital debris impacts with space station shielding. A digital recorder makes testing of high-speed water entries possible. Automobile engine design is aided by an instrumented head gasket that detects the combustion zone. And composite-to-metal strength and fatigue tests provide new data on joint failures in wind turbine joint tests.

The CONTAIN computer code is a best-estimate, integrated analysis tool for predicting the physical, chemical, and radiological conditions inside a nuclear reactor containment building following the release of core material from the primary system. CONTAIN is supported primarily by the U. S. Nuclear Regulatory Commission (USNRC), and the official code versions produced with this support are intended primarily for the analysis of light water reactors (LWR). The present manual describes CONTAIN LMR/1B-Mod. 1, a code version designed for the analysis of reactors with liquid metal coolant. It is a variant of the official CONTAIN 1.11 LWR code version. Some of the features of CONTAIN-LMR for treating the behavior of liquid metal coolant are in fact present in the LWR code versions but are discussed here rather than in the User`s Manual for the LWR versions. These features include models for sodium pool and spray fires. In addition to these models, new or substantially improved models have been installed in CONTAIN-LMR. The latter include models for treating two condensables (sodium and water) simultaneously, sodium atmosphere and pool chemistry, sodium condensation on aerosols, heat transfer from core-debris beds and to sodium pools, and sodium-concrete interactions. A detailed description of each of the above models is given, along with the code input requirements.

Inelastic material constitutive relations for elastoplasticity coupled with continuum damage mechanics are investigated. For elastoplasticity, continuum damage mechanics, and the coupled formulations, rigorous thermodynamic frameworks are derived. The elastoplasticity framework is shown to be sufficiently general to encompass J{sub 2} plasticity theories including general isotropic and kinematic hardening relations. The concepts of an intermediate undamaged configuration and a fictitious deformation gradient are used to develop a damage representation theory. An empirically-based, damage evolution theory is proposed to overcome some observed deficiencies. Damage deactivation, which is the negation of the effects of damage under certain loading conditions, is investigated. An improved deactivation algorithm is developed for both damaged elasticity and coupled elastoplasticity formulations. The applicability of coupled formulations is validated by comparing theoretical predictions to experimental data for a spectrum of materials and loads paths. The pressure-dependent brittle-to-ductile transitional behavior of concrete is replicated. The deactivation algorithm is validated using tensile and compression data for concrete. For a ductile material, the behavior of an aluminum alloy is simulated including the temperature-dependent ductile-to-brittle behavior features. The direct application of a coupled model to fatigue is introduced. In addition, the deactivation algorithm in conjunction with an assumed initial damage and strain is introduced as a novel method of simulating the densification phenomenon in cellular solids.

Target recognition requires the ability to distinguish targets from non-targets, a capability called one-class generalization. Many neural network pattern classifiers fail as one-class classifiers because they use open decision boundaries. To function as one-class classifier, a neural network must have three types of generalization: within-class, between-class, and out-of-class. We discuss these three types of generalization and identify neural network architectures that meet these requirements. We have applied our one-class classifier ideas to the problem of automatic target recognition in synthetic aperture radar. We have compared three neural network algorithms: Carpenter and Grossberg`s algorithmic version of the Adaptive Resonance Theory (ART-2A), Kohonen`s Learning Vector Quantization (LVQ), and Reilly and Cooper`s Restricted Coulomb Energy network (RCE). The ART 2-A neural network gives the best results, with 100% within-class, between-class, and out-of-class generalization. Experiments show that the network`s performance is sensitive to vigilance and number of training set presentations.

This report discusses the seventh experiment of the Integral Effects Test (IET-7) series. The experiment was conducted to investigate the effects of preexisting hydrogen in the Surtsey vessel on direct containment heating. Scale models of the Zion reactor pressure vessel (RPV), cavity, instrument tunnel, and subcompartment structures were constructed in the Surtsey Test Facility at Sandia National Laboratories. The RPV was modeled with a melt generator that consisted of a steel pressure barrier, a cast MgO crucible, and a thin steel inner liner. The melt generator/crucible had a hemispherical bottom head containing a graphite limitor plate with a 4-cm exit hole to simulate the ablated hole in the RPV bottom head that would be formed by ejection of an instrument guide tube in a severe nuclear power plant accident. The cavity contained 3.48 kg of water, and the containment basement floor inside the cranewall contained 71 kg of water, which corresponds to scaled condensate levels in the Zion plant. A 43-kg initial charge of iron oxide/aluminum/chromium thermite was used to simulate corium debris on the bottom head of the RPV. Molten thermite was ejected by steam at an initial pressure of 5.9 MPa into the reactor cavity.

Charpy V-notch specimens (ASTM Type A) and 5.74-mm diameter tension test specimens of the Shippingport Reactor Neutron Shield Tank (NST) (outer wall material) were irradiated together with Charpy V-notch specimens of the Oak Ridge National Laboratory (ORNI), High,, Flux Isotope Reactor (HFIR) vessel (shell material), to 5.07 {times} 10{sup 17} n/cm{sup 2}, E > 1 MeV. The irradiation was performed in the Ford Nuclear Reactor (FNR), a test reactor, at a controlled temperature of 54{degrees}C (130{degrees}F) selected to approximate the prior service temperatures of the cited reactor structures. Radiation-induced elevations in the Charpy 41-J transition temperature and the ambient temperature yield strength were small and independent of specimen test orientation (ASTM LT vs. TL). The observations are consistent with prior findings for the two materials (A 212-B plate) and other like materials irradiated at low temperature (< 200{degrees}C) to low fluence. The high radiation embrittlement sensitivity observed in HFIR vessel surveillance program tests was not found in the present accelerated irradiation test. Response to 288{degrees}C-168 h postirradiation annealing was explored for the NST material. Notch ductility recovery was found independent of specimen test orientation but dependent on the temperature within the transition region at which the specimens were tested.

A method is presented for determining the nonlinear stability of undamped flexible structures spinning about a principal axis of inertia. Equations of motion are developed for structures that are free of applied forces and moments. The development makes use of a floating reference frame which follows the overall rigid body motion. Within this frame, elastic deformations are assumed to be given functions of n generalized coordinates. A transformation of variables is devised which shows the equivalence of the equations of motion to a Hamiltonian system with n + 1 degrees of freedom. Using this equivalence, stability criteria are developed based upon the normal form of the Hamiltonian. It is shown that a motion which is spin stable in the linear approximation may be unstable when nonlinear terms are included. A stability analysis of a simple flexible structure is provided to demonstrate the application of the stability criteria. Results from numerical integration of the equations of motion are shown to be consistent with the predictions of the stability analysis. A new method for modeling the dynamics of rotating flexible structures is developed and investigated. The method is similar to conventional assumed displacement (modal) approaches with the addition that quadratic terms are retained in the kinematics of deformation. Retention of these terms is shown to account for the geometric stiffening effects which occur in rotating structures. Computational techniques are developed for the practical implementation of the method. The techniques make use of finite element analysis results, and thus are applicable to a wide variety of structures. Motion studies of specific problems are provided to demonstrate the validity of the method. Excellent agreement is found both with simulations presented in the literature for different approaches and with results from a commercial finite element analysis code. The computational advantages of the method are demonstrated.

We describe an algorithm for the static load balancing of scientific computations that generalizes and improves upon spectral bisection. Through a novel use of multiple eigenvectors, our new spectral algorithm can divide a computation into 4 or 8 pieces at once. These multidimensional spectral partitioning algorithms generate balanced partitions that have lower communication overhead and are less expensive to compute than those produced by spectral bisection. In addition, they automatically work to minimize message contention on a hypercube or mesh architecture. These spectral partitions are further improved by a multidimensional generalization of the Kernighan-Lin graph partitioning algorithm. Results on several computational grids are given and compared with other popular methods.

A robotic rover vehicle designed for use in the exploration of the Lunar surface is described. The Robotic All-Terrain Lunar Exploration Rover (R-A.T.L.E.R-) is a four wheeled all-wheel-drive dual-body vehicle. A uniquely simple method of chassis articulation is employed which allows all four wheels to remain in contact with the ground, even while climbing over step-like obstacles as large as [approximately]1.3 wheel diameters. Skid steering and modular construction are used to produce a simple, rugged, highly agile mobility chassis with fewer parts required compared to other designs being considered for planetary exploration missions. The design configuration, mobility parameters, and performance of several existing R.A.T.L.E.R prototypes are discussed, with emphasis on an analysis of the configuration parameters which directly affect the designs mobility performance.

This report describes preliminary experiments to investigate the feasibility of using electron beam (e-beam) radiolysis to destroy the organic compounds in simulated Hanford tank waste. For these experiments a simulated Hanford Tank 101-SY waste mixture was radiolyzed in a {sup 60}Co facility to simulate radiolysis in the waste tank. This slurry was then exposed without dilution to dose levels up to 1600 Mrad at instantaneous dose rates of 2.5 {times} 10{sup 8} and 2. 7 {times} 10{sup 11} rad/s. The inferred dose to destroy all the organic material in the simulated waste, assuming destruction is linear with dose, is 1000 Mrads for the higher dose rate. The cost for organic destruction of Hanford waste at a treatment rate of 20 gpm is roughly estimated to be $10. 60 per gallon. Such a system would treat all the waste in a 1 million gallon Hanford tank in about 40 days. Estimates of capital costs are given in the body of this report. While ferrocyanide destruction was not experimentally investigated in this work, previous experiments by others suggest that ferrocyanide would also be destroyed in such a system.

This report describes the Training and Qualification Program at the Saturn Facility. The main energy source at Saturn is the Saturn accelerator which is used to test military hardware for vulnerability to X-rays, as well as to perform various types of plasma radiation source experiments. The facility is operated and maintained by a staff of twenty scientists and technicians. This program is designed to ensure these personnel are adequately trained and qualified to perform their jobs in a safe and efficient manner. Copies of actual documents used in the program are included as appendices. This program meets all the requirements for training and qualification in the DOE Orders on Conduct of Operations and Quality Assurance, and may be useful to other organizations desiring to come into compliance with these orders.

Experiments were run to determine if oxidized Kovar could be chemically cleaned so that copper would wet the Kovar in a wet hydrogen atmosphere at 1100{degrees}C. We found that a multi-stepped acid etch process cleaned the Kovar so that copper would wet it. We also found that the degree of copper cracking after melting and cool-down correlated well with the degree of wetting.

The relatively thin web of salt that separates Bayou Choctaw Caverns 15 and 17 was evaluated using the finite-element method. The stability calculations provided insight as to whether or not any operationrestrictions or recommendations are necessary. Because of the uncertainty in the exact dimensions of the salt web, various web thicknesses were examined under different operating scenarios that included individual cavern workovers and drawdowns. Cavern workovers were defined by a sudden drop in the oil side pressure at the wellhead to atmospheric. Workovers represent periods of low cavern pressure. Cavern drawdowns were simulated by enlargening the cavern diameters, thus decreasing the thickness of the web. The calculations predict that Cavern 15 dominates the behavior of the web because of its larger diameter. Thus, giventhe choice of caverns, Cavern 17 should be used for oil withdrawal in order to minimize the adverse impacts on web resulting from pressure drops or cavern enlargement. From a stability point of view, maintaining normal pressures in Cavern 15 was found to be more important than operating the caverns as a gallery where both caverns are maintained at the same pressure. However, during a workover, it may be prudent to operate the caverns under similar pressures to avoid the possibility of a sudden pressure surge at the wellhead should the web fail.

A feasibility study for developing an improved tool and improved models for performing event assessments is described. The study indicates that the IRRAS code should become the base tool for performing event assessments, but that modifications would be needed to make it more suitable for routine use. Alternative system modeling approaches are explored and an approach is recommended that is based on improved train-level models. These models are demonstrated for Grand Gulf and Sequoyah. The insights that can be gained from importance measures are also demonstrated. The feasibility of using Individual Plant Examination (IPE) submittals as the basis for train-level models for precursor studies was also examined. The level of reported detail was found to vary widely, but in general, the submittals did not provide sufficient information to fully define the model. The feasibility of developing an industry risk profile from precursor results and of trending precursor results for individual plants were considered. The data sparsity would need to be considered when using the results from these types of evaluations, and because of the extremely sparse data for individual plants we found that trending evaluations for groups of plants would be more meaningful than trending evaluations for individual plants.

The discrete Fourier transform (DFT) is used frequently used in the computation of the signal-to-noise ratio (SNR) and harmonic distortion. To estimate the SNR or harmonic distortion, a sine wave is applied to the digitizing system under test. When the data record contains an integer number of cycles of the sine wave, energy from the sine wave and its harmonics does not leak into adjacent DFT frequency bins. Each harmonic occupies one an only one DFT frequency bin. To find the root-mean-square (RMS) value of a harmonic from its DFT, one computes the magnitude of the DFT value at the single frequency of the harmonic. When the DFTs of the fundamental and its harmonics are single lines, the SNR and harmonic distortion are easy to compute. When the data record contains a non-integer number of cycles of the sine wave, energy leaks from the sine wave and its harmonics to adjacent frequencies. The literature contains several approaches to problem of determining which DFT components correspond to a sine wave harmonic. This paper describes how to estimate the RMS value of a sine wave from its DFT with special attention to the selection of the DFT window. The set of DFT frequencies which comprise a harmonic depends on the DFT window, the length of the DFT, and the number of bits of the digitizer. Criteria are developed for choosing the DFT frequencies that correspond to a sine wave. These criteria lead to better choices of DFT windows for SNR and harmonic distortion calculations.

Sandia has developed an advanced operational control system approach, caged Graphical Programming, to design and operate robotic waste cleanup and other hazardous duty robotic systems. The Graphical Programming approach produces robot systems that are faster to develop and use, safer in operation, and cheaper overall than altemative teleoperation or autonomous robot control systems. The Graphical Programming approach uses 3-D visualization and simulation software with intuitive operator interfaces for the programming and control of complex robotic systems. Graphical Programming Supervisor software modules allow an operator to command and simulate complex tasks in a graphic preview mode and, when acceptable, command the actual robots and monitor their motions with the graphic system. Graphical Progranuning Supervisors maintain registration with the real world and allow the robot to perform tasks that cannot be accurately represented with models alone by using a combination of model and sensor-based control. This paper describes the Graphical Programming approach, several example control systems that use Graphical Programming, and key features necessary for implementing successful Graphical Programming systems.

Hazardous operations which in the past have been completed by technicians are under increased scrutiny due to high costs and low productivity associated with providing protective clothing and environments. As a result, remote systems are needed to accomplish many hazardous materials handling tasks such as the clean up of waste sites in which the exposure of personnel to radiation, chemical, explosive, and other hazardous constituents is unacceptable. Traditional remote manual operations have proven to have very low productivity when compared with unencumbered humans. Computer models augmented by sensing and structured, modular computing environments are proving to be effective in automating many unstructured hazardous tasks.

Testing and analysis of shock wave characteristics such as produced by detonators and ground shock propagation frequently require a method of measuring velocity and displacement of the surface of interest. One method of measurement is doppler interferometry. The VISAR (Velocity Interferometer System for Any Reflector) uses doppler interferometry and has pined wide acceptance as the preferred tool for shock measurement. An important asset of VISAR is that it measures velocity and displacement non intrusively. The conventional VISAR is not well suited for portability because of its sensitive components, large power and cooling requirements, and hazardous laser beam. A new VISAR using the latest technology in solid state lasers and detectors has been developed and tested. To further enhance this system's versatility, the unit is fiber optic coupled which allows remote testing, permitting the VISAR to be placed over a kilometer away from the target being measured. Because the laser light is contained in the fiber optic, operation of the system around personnel is far less hazardous. A software package for data reduction has also been developed for use with a personal computer. These new advances have produced a very versatile system with full portability which can be totally powered by batteries or a small generator. This paper describes the solid state VISAR and its peripheral components, fiber optic coupling methods and the fiber optic coupled sensors used for sending and receiving laser radiation.

Thick film Au metallizations are commonly used as conductors in hybrid microelectronics that operate at high frequencies. Discrete components are attached to these conductors with 5OPb/5OIn solder. Intermetallic compounds form and grow in the solid state with time; AuIn[sub 2] is the primary compound formed in the Au-5OPb/50In system. A hybrid failed after being artificially aged to consume all of the Au and then subjected to normal thermal cycle and vibration testing. Postmortem analysis revealed that three capacitors had debonded. The failed parts were studied to determine the failure mechanism(s) and define a parametric study to characterize the mechanisms that bond the Au thick film to the ceramic substrate.

This is a peer review report of the Environmental and Molecule Sciences Laboratory. Although the Pretreatment, Treatment, and Waste Forms comments are focused specifically on pretreatment treatment, and waste forms, the group recognizes that the life cycle designation is a somewhat arbitrary breakdown of a series of activities that form a continuum in the environmental restoration and waste management program. Consequently, some of the comments made here are relevant in a broader context or even for EM as a whole. particular, characterization activities pervade all life cycles in environmental restoration/waste management. As we use the term in this section, characterization'' refers to the process monitoring and control that are required during pretreatment and treatment. Most of the technology presentations during the review identified links to support this area and delineated to varying degrees the specific ties to the Hanford Site cleanup requirements. Overall, the EMSL especially its planned facilities, are most impressive, and DOE and PNL are urged to proceed with all due haste toward its completion. Specific issues or concerns identified during the review are included in the following section.

Scientific visualization is playing an increasingly important role in the analysis and interpretation of massively parallel CFD simulations due to the enormous volume of data that can be generated on these machines. In this paper we will describe the development of a visualization technique based on a parallel analogue to the Marching Cubes algorithm. The algorithm has been developed for Multiple-Instruction, Multiple-Data (MIMD) massively parallel computers and is designed to take advantage of the heterogeneous programming capabilities of the MIMD architecture. We examine several different configurations and conclude that for producing animations the best one, in terms of both frame generation time and disk usage, is to run the two applications heterogeneously and send the resulting geometry description directly to a workstation for rendering, thereby totally eliminating the use of files from the animation process.

Sandia National Laboratories is prototyping a service for providing a distributed visualization capability between its Albuquerque, New Mexico and Livermore, California sites. The service will allow TCP/EP LAN users to process computationally intensive codes on the Cray Y-MP 8/864 and use the Application Visualization System, (AVS) running on a visualization server to display the results to desktops in Livermore (or Albuquerque). The long-haul environment presents several technical challenges such as high delays (up to 42 milliseconds), potential error rates, and security concerns which can affect the quality of the services. The network design requires low latency switch gear and high speed LAN connections (FDDI) to make the visualization/computational services useful to the customer. This paper will describe tie network design used and will report performance characteristics of the applications utilizing this service.

AC impedance and x-ray diffraction measurements versus temperature and pressure ternary intercalation compounds are reported. Results support our previously established empirical correlation between superconducting onset temperature and 300K fcc lattice constant. Compounds which do not follow this correlation either phase separate or are otherwise unstable at low T and/or high P.

We have demonstrated removal of H from a H-passivated Si(100) surface by low energy ([approximately]50 eV) He ion bombardment. The extent of the removal of H from the surface can be controlled by varying the duration of He bombardment and plasma parameters. This, in turn, means that the growth rate by this method can also be controlled. Si[sub 2]H[sub 6] was shown to adsorb in a self-limiting manner on the Si(100) surface.

The earth`s ionosphere consists of an ionized plasma which will interact with any electromagnetic wave propagating through it. The interaction is particularly strong at vhf and uhf frequencies but decreases for higher microwave frequencies. These interaction effects and their relationship to the operation of a wide-bandwidth, synthetic-aperture, space-based radar are examined. Emphasis is placed on the dispersion effects and the polarimetric effects. Results show that high-resolution (wide-bandwidth) and high-quality coherent polarimetrics will be very difficult to achieve below 1 GHz.

SAFSIM (System Analysis Flow SIMulator) is a FORTRAN computer program that provides engineering simulations of user-specified flow networks at the system level. It includes fluid mechanics, heat transfer, and reactor dynamics capabilities. SAFSIM provides sufficient versatility to allow the simulation of almost any flow system, from a backyard sprinkler system to a clustered nuclear reactor propulsion system. In addition to versatility, speed and robustness are primary goals of SAFSIM development. The current capabilities of SAFSIM are summarized and some sample applications are presented. It is applied here to a nuclear thermal propulsion system and nuclear rocket engine test facility.

Photovoltaic (PV) systems offer a cost-effective solution to provide electrical power for a wide variety of applications, with battery performance playing a major role in their success. This paper presents some of the results of an industry meeting regarding battery specifications and ratings that photovoltaic system designers require, but do not typically have available to them. Communications between the PV industry and the battery industry regarding appropriate specifications have been uncoordinated and poor in the past. This paper also discusses the effort under way involving the PV industry and battery manufacturers, and provides a working draft of specifications to develop and outline the information sorely needed on batteries. The development of this information is referred to as ``Application Notes for Batteries in Photovoltaic Systems.`` The content of these ``notes`` has been compiled from various sources, including the input from the results of a survey on battery use in the photovoltaic industry. Only lead-acid batteries are discussed

In one design of molten-salt central receivers, the molten salt flows in a serpentine path, down one panel of tubes then up the next and down again continuing in this fashion through the receiver. There have been concerns about this design because in the down flow sections, the heat flux incident on the tubes can cause flow instability since the flow is in direct opposition to the buoyant forces. In extreme cases the buoyant forces can cause flow stagnation or reversal. An analysis of flow stability within individual tubes and down flow sections of receiver panels is presented. When the partial derivative of the pressure drop with respect to mass flow rate is negative ({partial_derivative}{Delta}P/{partial_derivative}{sup {lg_bullet}} < 0), the flow is unstable and could cause serious damage to the receiver. Stability maps are developed that show safe operating regimes where inertial forces dominate over buoyant forces. The data is then normalized using the Grashof and Reynolds numbers.

United States Department of Energy has established its first Advanced Research Objective in the Solids Transport Program. The scientific, engineering, and management goals are discussed in some detail. Scientific progress to date is summarized. Comments are made on the technical direction of further Advanced Research Objectives.

We have developed a Li/SOCl{sub 2} ``D`` cell for applications requiring 10 to 15 years life at very low drain rates, typically less than 150 {mu}A. Maximizing cell safety and reliability, while delivering very good energy density, have been the goals of our study. We have achieved these goals by designing the cell to be application specific. The low-rate cell has been optimized to deliver up to 16 Ah at drain rates of less than 70 mA. By virtue of its low surface area, 145 cm{sub 2}, the cell has demonstrated excellent safety behavior. Safety testing has been performed on individual cells as well as on two-cell and four-cell batteries. Single cells did not vent when short-circuited. We were able to produce benign venting in a two cell string, but only when the string was partially discharged before shorting. The vent mechanism is a 300 psi rupture pressure burst disc manufactured by BS&B Safety Systems. We define benign venting as full opening of the 3/8 in. dia vent hole without deformation of the case. Material is expelled from the cell without flame, and the cell stack remains largely intact. We have not produced venting of the Sandia-designed low rate cell under any other abuse test conditions. The vent functions as an ultimate safety mechanism in the case of severe abuse, but resistance to venting under normal use and mild abuse conditions is key to the achievement of high reliability.

A new class of ion-exchange materials that can selectively separate low parts per million level concentrations of Cs{sup +} from 3--6 molar concentrations of Na{sup +} over a wide pH range has recently been developed as a result of a collaborative effort between Sandia National Laboratories and Texas A&M University. The materials, called crystalline silicotitanates, show potential for application in the treatment of aqueous nuclear waste solutions.

In order to develop a procedure for measuring cation diffusion coefficients below 1000{degrees}C, we have examined the suitability of several diffusion couple configurations involving single crystals of garnet. Initial experiments using an enriched {sup 25}MgCl{sup 2} proved ineffective in providing a uniform and coherent surface for analysis by ion microprobe. A technique was developed using thin film deposition. Thin films ({approximately} 1000 {Angstrom}) of MgO{sub x} (x < 1) can be applied to polished mineral surfaces by evaporating MgO powder under high vacuum with a thermal-resistance strip heater. Thermal resistance evaporation is efficient. Samples of single crystal grossular and pyrope garnets with thin films of MgO, as created by these techniques, were annealed for various times at 800, 900, and 1000{degrees}C, at several log fO{sub 2} values, and 1 atm. Optical, SEM, and ion microprobe analyses reveal no disruption of the interface. Profiles of elemental counts vs depth exhibit expected patterns going through the thin film into the garnet substrate. Our experimental matrix of garnet diffusion runs includes over 60 cut and polished crystals of pyrope composition that are being run at various oxygen fugacity conditions from 600 to 1000{degrees}C.

The Los Alamos National Laboratory (LANL) Mechanical Engineering and Electronics Division, in partnership with Sandia National Laboratories and Programmed Composites, is advancing the development of thin-walled, high modulus short-fiber compression-molded composite materials fabrication. In this paper, we investigate component uniformity, structural integrity, thermal conductivity, and radiation resistance; discuss the scanning-electron microscopic inspection of the graphite fiber distribution and orientation, and describe the process used in selecting the reinforcement fiber length and modulus and for choosing the hydrophobic, cyanate-ester resin.

During the last two decades there has been considerable interest in developing alternatives to conventional chemical propulsion for space missions. Laser propulsion has been identified as a serious contender for the task of inexpensively delivering small payloads to low-earth orbit. Recent advances in the development of lasers powered directly by nuclear reaction products offer the potential for new propulsion methods, namely, reactor-laser propulsion. Such systems would allow ``nuclear propulsion`` without placing nuclear systems in space.

Alternate organometallic Sb sources are being investigated to improve the characteristics of InSb grown by MOCVD. InSb grown using trimethylindium (TMIn) and trimethylantimony (TMSb) or triethylantimony (TESb) yielded similar quality materials under similar growth conditions. InSb grown using triethylindium (TEIn) and TESB under similar growth conditions yielded very poor quality n-type material. Three new organometallic Sb sources, triisopropyl-antimony (TIPSb), tris(dimethylamino)antimony (TDMASb), and tertiarybutyldimethylantimony (TBDMSb) are being investigated. Growth of InSb using TIPSb, TDMASb, or TBDMSb and TMIn was investigated over 350 to 475{degrees}C. InSb grown from TDMASb had similar properties to InSb grown from TMIn and TMSb when using a similar temperature and V/III ratio range. Growth rates of InSb using TMIn and either TIPSb or TBDMSb at temperatures {le} 425{degrees}C were proportional to both TMIn flow rate and temperature. Surface morphology of InSb grown using either TIPSb or TBDMSb was rough for growth temperatures {le} 425{degrees}C; this may be due to complex decomposition and methyl groups on surface. The InSb with the highest mobility was grown at 400{degrees}C and a V/III ratio of 3 using TIPSb. It was n-type with a carrier concentration of 2.5 {times} 10{sup 15} cm{sup {minus}3} and a mobility of 78,160 cm{sup 2}/Vs at 77 K. Both n- and p-type InSb were grown using TBDMSb with mobilities up to 67,530 and 7773 cm{sup 2}/Vs, respectively at 77 K. Mobility for InSb using either TIPSb or TBDMSb was optimized by going to lower temperatures, pressures, V/III ratios; however, surface morphology improved with higher temperature, pressure, V/III ratio. High mobility InSb with smooth surfaces at T {le} 425{degrees}C was not obtained with TIPSb or TBDMSb and TMIn.

We present the use of an STM to make quantitative observations of time-dependent mass flow associated with the decay of two-dimensional clusters on the Au(lll) surface. When formed and observed in air, layered islands with well-defined edges located on larger terraces are generally found to decay in such a way that their areas decrease linearly in time over periods ranging from minutes to several hours depending on the island size. This is in contrast to the behavior of similar features formed and observed under ultra high vacuum conditions, which do not appear to decay over experimental periods of several days. The linear decay is consistent with models that have been used previously to describe growth of 2-dimensional clusters on surfaces. We discuss possible decay mechanisms`, and the role that adsorbates may play in influencing the decay.

In support of the UST, WeDID, VOC/Non-Arid, and VOC/Arid, and VOC/Arid Integrated Demonstrations, organic contaminants and nitrates in Hanford Underground Storage Tank (UST) wastes, polymeric organics in weapon components, and chlorocarbon soil contaminants have been destroyed by exposure to high-temperature steam during bench tests with a quartz reactor and full-scale tests that used the Synthetica Detoxifier, a commercial one-ton-per-day steam reforming waste destruction system. Reactivation of Granular Activated Carbon (GAC) in the Detoxifier and Thermal Gravimetric Analyses (TGA) of the decomposition of sodium nitrate have also been performed.

We have successfully grown ZnMnTe alloys by molecular beam epitaxy using GaAs as a substrate. Bulk MnTe has the wurtzite crystal structure but the structural phase of the material was confirmed to be zinc-blende by standard {theta}-2{theta} x-ray diffraction techniques. The composition was also determined using x-ray diffraction techniques. Manganese concentration was also estimated from magnetization measurements taken as a function of temperature. Magnetoluminescence studies were performed at 1.4K on the acceptor-bound exciton in the semimagnetic semiconductor ZnMnTe alloys. As expected, the photoluminescence peak energy decreased with increasing magnetic field.

At some landfills, hazardous wastes were placed into disposal pits with other waste oils. Volatile organic wastes, such as cleaning solvents, when combined with oil are much more difficult to remove from the soils because the oil makes the volatile organic chemical evaporate much slower. The typical vacuum extraction remediation method could become a lengthy campaign. Since all chemicals evaporate faster when they are heated, if the contaminated soil could be heated, the chemicals would be easier to remove from the soil. By using heating techniques developed by the Oil and Gas Industry to enhance the removal of oil and gas resources from the soil, the problem of removing contamination from the soils could be solved. The Thermal Enhanced Vapor Extraction System (TEVES) demonstration will combine vacuum vapor extraction technology with powerline frequency soil heating and radiofrequency soil heating to accelerate the soil decontamination process. The premise that soil heating technology can actually reduce the costs of soil decontamination is complicated by the high capital costs of the soil heating equipment and the cost of electrode installation. By performing this field demonstration, Sandia will be able to collect the information needed to see if this new technology will improve the decontamination of soils.

Sandia National Laboratories is a large multiprogram Research and Development laboratory which is operated by a contractor for the US Department of Energy. In the Laboratories, lasers are both the subject of research and the tools that are used in other research, development, and testing activities. Since 1979, laser safety training has been the primary focus of the Laboratories` laser safety program. Approximately 1 100 personnel have been trained in formal courses during that time period. The formal course, presented on site by a contractor, consists of two full days of instruction. The course contents include the following topics: Laser technology and safety overview; Federal and ANSI laser standards summaries; Biological effects of laser radiation; Classification of lasers; Laser hazard analysis; Review of ANSI Z136.1 control measures; Laser eye protection. Recent emphasis on regulatory requirements, conduct of operations, and quality management has revealed a need to change the laser safety training curriculum. A new course for users of low power lasers (Class 2 and 3a) is being developed. A refresher course, a management awareness (self-study) course, and major changes in the current course are planned.

Short communication.

A solar photocatalytic process has been under development at both Sandia National Laboratories and the National Renewable Energy Laboratory (formerly the Solar Energy Research Institute). This process uses solar ultraviolet light to activate a titanium dioxide catalyst which oxidizes organic contaminants in water. In the summer of 1991, a solar photocatalytic detoxification of water system was installed and tested at a California Superfund Site located at Lawrence Livermore National Laboratory. The site was designated a Superfund Site because of widespread groundwater contamination which resulted from the release of chlorinated solvents, principally trichloroethylene, when the site was a Naval Air Station in the early 1940s. The objectives of these experiments were to measure the effects of process variables and the process efficiency in an actual remediation setting, to collect experimental data and operating experience in photocatalytic oxidation of organic contaminants, to develop accurate models of the system operation and to develop control strategies.

The Materials Interface Interaction Tests (MIIT) program involves the comparative performance-evaluation testing of multiple US and foreign nuclear waste glasses (nonradioactive), potential canister and overpack metals, brine, and geologic materials in the rock salt repository environment at the Waste Isolation Pilot Plant (WIPP) facility. We emplaced about 2000 materials specimens onto fiftn, separate test assemblies and exposed them to a heated, salt-brine environment at the WIPP for multi-year periods. We successfully terminated the in situ conduct of the MIIT in July 1991, after five years of testing, and retrieved all samples for posttest laboratory analyses. These 5-year glass and metal samples, along with samples previously retrieved after 0.5, 1, and 2 years, are being analyzed in multiple international laboratories, in a cooperative testing effort. Individual test participants will present available laboratory results, and interpretations, on MIIT specimens in this workshop. Our focus in this paper is to summarize technical details and repository-relevant observations on the in situ conduct, sampling, and termination operations of the MIIT experimental program. Such information should be useful for the interpretation of the laboratory-based analyses. This information also will be relevant and instructive for other organizations contemplating, planning, or conducting additional materials-related, in situ tests.

A comparison between numerical calculations with use of commercial thermal analysis software packages and experimental data simulating a horizontally oriented spent fuel rod array was performed. Twelve cases were analyzed using air and helium for the fill gas, with three different heat dissipation levels. The numerically predicted temperatures are higher than the experimental data for all levels of heat dissipation with air as the fill gas. The temperature differences are 4{degree}C and 23{degree}C for the low heat dissipation and high heat dissipation, respectively. The temperature predictions using helium as a fill gas are lower than the experimental data for the low and medium heat dissipation levels. The temperature predictions are 1{degree}C and 6{degree}C lower than the experimental data for the low and medium heat dissipation, respectively. For the high heat dissipation level, the temperature predictions are 16{degree}C higher than the experimental data. Differences between the predicted and experimental temperatures can be attributed to several factors. These factors include experimental uncertainty in the temperature and heat dissipation measurements, actual convection effects not included in the model, and axial heat flow in the experimental data. This work demonstrates that horizontally oriented spent fuel rod surface temperature predictions can be made using existing commercial software packages. This work also shows that end effects, such as axial heat transfer through the spent fuel rods, will be increasingly important as the amount of dissipated heat increases.

Magnetic force microscopy (MFM) has been applied to image currents in internal IC conductors. We present a model for the MFM imaging of IC currents, describe MFM signal generation, and demonstrate the ability to analyze current direction and magnitude with a sensitivity of {approximately} 1 mA dc and {approximately} 1 {mu}A ac. Our experimental results are a significant improvement on the 100 mA ac resolution previously reported using an electron beam to detect IC currents [1].

The Waste Isolation Pilot Plant: (WIPP) in southeastern New Mexico is being developed by the US Department of Energy as a disposal facility for transuranic waste. In support of this project, Sandia National Laboratories is conducting an ongoing performance assessment (PA) for the WIPP. The ordered triple- representation for risk proposed by Kaplan and Garrick is used to provide a clear conceptual structure for this PA. This presentation describes how the preceding representation provides a basis in the WIPP PA for (1) the definition of scenarios and the calculation of scenario probabilities and consequences, (2) the separation of subjective and stochastic uncertainties, (3) the construction of the complementary cumulative distribution functions required in comparisons with the US Environmental Protection Agency`s standard for the geologic disposal of radioactive waste (i.e., 40 CFR Part 191, Subpart B), and (4) the performance of uncertainty and sensitivity studies. Results obtained in a preliminary PA for the WIPP completed in December of 1991 are used for illustration.

Short communication.

As applications for hybrid circuits and multichip modules create demand for higher density circuits and higher power components, new substrate materials are required to deal with the heat generated on the circuit. Sandia National Laboratories is developing diamond substrate technology to meet the requirements of high thermal conductivity. Thin film processes were developed and characterized to delineate conductor-resistor networks on free standing diamond substrates having fine line gold conductors and low and high sheet resistivity resistors. Thin film hybrid circuit technology was developed on CVD-processed, polycrystalline diamond substrates having as-deposited surface finishes as well as those with polished surfaces. Conductors were defined by pattern plating gold and resistors were processed from sputtered tantalum nitride films which were deposited to sheet resistivities of 5 and/or 100 ohms per square. Resistor films on diamond substrates were evaluated for Temperature Coefficient of Resistance (TCR), stability with time and temperature, and trimmability using YAG laser processing. Plated gold conductors were patterned on diamond to feature sizes of 25 microns and successfully tested for adhesion and bondability. Advanced YAG laser trimming techniques were developed to allow resistor trims on both low and high value resistors to within 1% of desip value while maintaining required resistor stability, new trim techniques were needed to offset the carbonization of diamond in the laser trim area. Reliability studies were carried out on the diamond thin film networks which showed them to compare favorably with the same thin film technology on alumina substrates.

Photocreation mechanisms and properties of nitrogen dangling bonds in amorphous hydrogenated silicon nitride (a-SiN{sub x}:H) thin films are investigated. We find that the creation kinetics are strongly dependent on the post-deposition anneal; this thermal process can be described by a simple exponential function which yields an activation energy of 0.8 eV. The compositional dependence of the nitrogen dangling bond center suggests that its energy level lies close to the valence band edge, in agreement with theoretical calculations. This energy level position can explain why a-SiN{sub x}:H films often become conducting following a high post-deposition anneal.

Path planning needs to be fast to facilitate real-time robot programming. Unfortunately, current planning techniques are still too slow to be effective, as they often require several minutes, if not hours of computation. To overcome this difficulty, we present an adaptive algorithm that uses previous experience to speed up future performance. It is a learning algorithm suitable for incrementally-changing environments such as those encountered in manufacturing of evolving products and waste-site remediation. The algorithm extends our previous work for stationary environments in two directions: For minor environmental change, an object-attached experience abstraction scheme is introduced to increase the flexibility of the learned experience; for major environmental change, an on-demand experience repair scheme is also introduced to retain those experiences that remain valid and useful. In addition to presenting this algorithm, we identify three other variants with different repair strategies. To analyze the respective performance of these algorithms, we develop an analytic model that quantifies and relates training effort, experience value and utility, and environmental change through intuitive terms of energy and work. It is a general and simple model that should be very useful in characterizing other types of learning processes as well. Using this model, we formalize the concept of incremental change, and prove the optimality of our proposedalgorithm under such change. Empirically, we also characterize the performance curve of each variant, confirm our theoretical optimality results, and demonstrate the practicality of our algorithm.

A method of solving the two-phase fluid flow equations using a genetic algorithm on a NCUBE multiprocessor computer is presented. The topics discussed are the two-phase flow equations, the genetic representation of the unknowns, the fitness function, the genetic operators, and the implementation of the algorithm on the NCUBE computer. The efficiency of the implementation is investigated using a pipe blowdown problem. Effects of varying the genetic parameters and the number of processors are presented.

A method is described for parallelizing molecular dynamics (MD) simulations by block-decomposing the matrix of bonded and non-bonded force computations. It is particularly useful for organic simulations (polymers, proteins) because unlike spatial-decomposition methods, it requires no geometric information about the simulation domain. Because its communication cost scales as N/{radical}P. rather than N as in the all-to-all broadcast or ring-exchange techniques commonly used in this type of MD simulation, larger numbers of processors can be used effectively, yielding greater parallel speed-ups.

Two major initiatives are underway in the US that are creating a significant financial impact on both the US taxpayer and on users of electric power. First, the US Department of Energy (DOE) has been tasked with cleaning-up the defense complex. This task is managed under the direction of the Office of Environmental Restoration and Waste Management (EM) of the DOE. The waste that EM must address includes radioactive, hazardous, and mixed that consists of both radioactive and hazardous constituents. Second, the DOE is required by the Nuclear Waste Policy Act (NWPA) to take title to commercial nuclear spent fuel assemblies starting in 1998. The DOE Office of Civilian Radioactive Waste Management (OCRWM) was established to carry out this charter. Since a final repository is not scheduled for opening until 2010 at the earliest, the DOE is planning on providing a Monitored Retrievable Storage (MRS) facility for centralized storage to bridge the time gap between 1998 and 2010. The NWPA requires that nuclear utilities pay a fee into a specific fund that Congress uses to pay the DOE for the development of the MRS, the transportation system, and the repository. This fund, along with the EM budget, constitutes a multi-billion dollar effort to manage DOE nuclear waste and to store and dispose of commercial spent nuclear fuel. These two seemingly unrelated problems have aspects of commonality that can be considered for the benefit of both programs, the US taxpayer, and the utility rate payer. Both programs are the responsibility of the DOE, and both will require engineered packages for storage, transportation, and disposal of the EM waste and commercial spent fuel. Rather than using specialized systems for each step (storage, transport, and disposal), a concept for a Universal Container System has been developed that could potentially simplify the overall waste management system, reduce expensive handling operations, and reduce total system cost.

Inadvertent alloying of Cu braze metal can compromise metal/ceramic seals. Electron microprobe analyses have quantified alloying of Cu brazes in metal/ceramic feedthroughs. Pin material and processing parameters above 1084C both affect alloying levels. Using either Kovar or Ni-plated 316L stainless steel pins limits alloying compared to Palco pins. Minimizing the time during which the braze is molten also avoids excessive alloying. The original thickness of the Ni plating on the Mo-Mn metallization of the ceramic also influences the alloying content of these brazes. Metal/ceramic brazes made with long brazing cycles, Mo-Mn metallization, and Kovar components grow a layer of Mo{sub 6}(Fe{sub 3.5}CO{sub 3.5}){sub 7} on the metallization. Layer thicknesses observed do not appear to compromise joint integrity. Ni additions of approximately 10 and 20 wt.% to Cu apparently increases the stress required for stress relaxation during cooldown. to maintain creep rates required for stress relaxation during cooldown. Relative to unalloyed Cu, this strengthening effect tends to increase as temperature is decreased.

In this paper we investigate the applicability of the feature extraction mechanisms found in the neurophysiology of mammals to the problem of object recognition in synthetic aperture radar imagery. Our approach is to present multiple views of objects to be recognized to a two-stage self-organizing neural network architecture. The first stage, a two-layer Neocognitron, performs feature extraction in each layer The resulting feature vectors are presented to the second stage, an ART-2A classifier self-organizing neural network which clusters the features into multiple object categories. The feature extraction operators resulting from the self-organization process are compared to the feature extraction mechanisms found in the neurophysiology of vision. In a previous paper, the Neocognitron was trained on raw SAR imagery. The architecture was able to recognize a simulated vehicle at arbitrary azimuthal orientations at a single depression angle while rejecting clutter as well as other vehicles. Feature extraction on raw imagery yielded features that were robust but very difficult to interpret. In this paper we report the results of some new experiments in which the self-organization process is applied separately to shadow and bright returns from objects to be recognized. Feature extraction on shadow returns yield oriented contrast edge operators suggestive of bipartite simple cells observed in the striate cortex of mammals. Feature extraction on the specularity patterns in bright returns yield a collection of operators resembling a twodimensional Haar basis set. We compare the performance of the earlier two-stage neural network trained on raw imagery with a modified network using the new feature set.

The Mixed Waste Landfill Integrated Demonstration (MWLID) is testing noninvasive site characterization methods at several locations, including the Chemical Waste Landfill (CWL) at Sandia National Laboratories. The CWL comprises shallow, unlined pits that were used for the disposal of acids, oils, solvents, and inorganic compounds from 1962 until the CWL was closed in 1985. The soils of the landfills are alluvial, predominantly sand, gravels and cobbles with small quantities of silts and clays. The focus of this study is an unlined chromic acid pit (UCAP). The UCAP pit is rectangular (approximately 4 {times} 10 m), which reportedly received unknown volumes of chromium in the form of chromic acid (liquid) and other hazardous materials. At this location, we have demonstrated a continuous waveform (CW) system for site characterization. During this year, we will also utilize a crossborehole pulsed radar system. Both methods are sensitive to variations in either electrical conductivity or dielectric constant in the soils or host rock at a waste site. These earth properties are some of the most responsive geophysical indicators of metallic, acidic and water-based subsurface contaminants.

Using Saturn as a driver, we are pursuing both photoresonantly pumped and photoionization/recombination lasers. Our lasing targets are gas cells with thin windows that are pumped by a z pinch 2 cm away radiating 10 TW. In both schemes the lasant and gas fill is neon. To increase our chances of measuring the resonantly photopumped lasing transition we have introduced potassium into a sodium z pinch and have eliminated oxygen from the gas cell windows. We have measured the spatial dependence of ionization balance across the gas cell, and this measurement is consistent with propagation of a shock front across the gas cell target. We have measured blue-shifted satellites to several Li-like neon transitions that may indicate return-current driven jetting a high 1.5e8 cm/sec velocity. Using a gold z-pinch we have shown that key radiation is necessary to excite the He-like lines of neon. An attempt at a single shot gain measurement also indicates that radiation is not the only source of gas cell heating.

This paper presents a full-field dye concentration measurement technique that extends our experimental capabilities to the measurement of transient dye concentration fields within steady state flow fields under unsaturated or saturated conditions. Simple light absorption theory provides a basis for translating images into high resolution dye concentration fields. A series of dye pulse experiments that demonstrate the combined use of the full-field saturation and dye concentration techniques was conducted at four different degrees of saturation. Each of these experimental sequences was evaluated with respect to mass balance, the results being within 5% of the known dye mass input. An image windowing technique allowed us to see increased dispersion due to decreasing moisture content, tailing of concentration at the rear of the dye pulse and slight velocity changes of the dispersive front due to changes in moisture content. The exceptional resolution of dye concentration in space and time provided by this laboratory technique allows systematic experimentation for examining basic processes affecting solute transport within saturated/unsaturated porous media. Future challenges for this work will be to use these techniques to analyze more complex systems involving heterogeneities, scaling laws, and detailed investigations of the relationship between transverse and longitudinal dispersion in unsaturated media.

The Beam Characterization System is being employed at the Sandia`s National Solar Thermal Test Facility to characterize the optical performance of heliostats, point-focus solar collectors, and their optical sub-elements as part of the on-going task to develop solar thermal technologies. With this measurement system, images of concentrated solar flux are acquired using digital imaging and processed to obtain such measures of the collector`s optical performance as beam power, flux distribution, and beam diameter. Key system elements are a diffusely reflective target for imaging collector beams, meteorological instrumentation including a flux gauge to measure flux at one point in the beam, and a calibration technique to establish a pixel-intensity-to-flux-density conversion factor for the image. The system is employed in a variety of collector tests such as beam quality, tracking error, and wind effects. The paper describes the Beam Characterization System and its components and presents, for illustration, sample test results. An analysis of the Beam Characterization System`s sources of measurement error is presented. Lastly, measurement techniques that employ the BCS to align heliostats and to measure or estimate collector surface slope errors are described.

We have measured the laser emission spectra of several vertical cavity surface emitting lasers following pulsed laser excitation, with a time resolution of < 1 ps. Correlations between the observed pulse widths and cavity lifetimes were observed.

A series of constant strain rate, unconfined compression experiments was performed on saturated welded tuff specimens collected from Busted Butte near Yucca Mountain, Nevada. Twenty specimens were loaded to failure at strain rates ranging from 10{sup {minus}9}s{sup {minus}1} to 10{sup {minus}3}s{sup {minus}1}, under ambient pressure and temperature conditions. The strength of the specimens showed a continuous decrease with decreasing strain rate between 10{sup {minus}9} s{sup {minus}1} and 10{sup {minus}5} s{sup {minus}1}. At the highest strain rate, 10{sup {minus}3} s{sup {minus}1}, strengths were less than those observed at 10{sup {minus}5} s{sup {minus}1}, likely due to hydrofracturing within the specimen at rapid loading rates. Reduction in strength, corresponding to the decrease in strain rate, is explained in terms of stress corrosion cracking. A detailed examination of six specimens tested at a strain rate of 10{sup {minus}9} s{sup {minus}1}, using acoustic wave velocities and CT scans, shows a correlation between the nature of the microstructure of the specimens and the observed strengths and elastic moduli.

Pore scale invasion percolation theory is modified for imbibition of.wetting fluids into fractures. The effects of gravity, local aperture field geometry, and local in-plane air/water interfacial curvatureare included in the calculation of aperture filling potential which controls wetted structure growth within the fracture. The inclusion of gravity yields fingers oriented in the direction of the gravitational gradient. These fingers widen and tend to meander and branch more as the gravitational gradient decreases. In-plane interfacial curvature also greatly affects the wetted structure in both horizontal and nonhorizontal fractures causing the formation of macroscopic wetting fronts. The modified percolation model is used to simulate imbibition into an analogue rough-walled fracture where both fingering and horizontal imbibition experiments were previously conducted. Comparison of numerical and experimental results showed reasonably good agreement. This process oriented physical and numerical modeling is-a necessary step toward including gravity-driven fingering in models of flow and transport through unsaturated, fractured rock.