Publications

SrCo{sub 0.5}FeO{sub x} (SCF) is an attractive material for oxygen separation membranes and for use in catalytic membrane reactors. While tubes of this material have been prepared by extrusion, further improvements in oxygen transport performance may be gained by preparing thinner membranes on porous supports. In this paper, we will discuss the deposition of thick films by spray deposition and centrifugal casting, and thin films by pyrolysis of chemical precursors. For the chemically derived thin films, porous MgO supports were used as membrane supports. Three types of precursor solutions were employed for dipcoating: a Pechini type solution, a nitrate-based solution, and a citrate-based solution. To prevent the infiltration of the precursor into the support, the support was backfilled with a material that decomposed at higher temperatures than the precursors. Cracking due to the volume changes during drying and pyrolysis of the precursors is discussed. Thick films were prepared by spray coating and centrifugal casting. Spray deposition of thick film membranes was accomplished by air brushing SCF from a water-based suspension onto the surface of a porous MgO support. Films on the interior surface of the supports were prepared by centrifugal casting using a xylene/butanol-based SCF suspension. Unlike extruded tubes, thick films undergo constrained sintering due to the presence of the support, which greatly reduces the densification rate. For membranes prepared by both approaches, we will discuss the effects of heating schedules on membrane microstructure, densification behavior, and cracking.

Plasma-induced-damage often degrades the electrical and optical properties of compound semiconductor devices. Despite the fact that the binding energy of GaN is larger than that for more conventional III--V compounds, etch damage is still a concern. Photoluminescence measurements and atomic force microscopy have been used to determine the damage induced in GaN by exposure to both electron cyclotron resonance (ECR) and inductively coupled plasmas (ICP) generated Ar plasmas.

This paper presents a portion of the work on specification, design, and implementation of safety-critical systems such as reactor control systems. A natural approach to this problem, once all the requirements are captured, would be to state the requirements formally and then either to prove (preferably via automated tools) that the system conforms to spec (program verification), or to try to simultaneously generate the system and a mathematical proof that the requirements are being met (program derivation). An obstacle to this is frequent presence of partially defined operations within the software and its specifications. Indeed, the usual proofs via first order logic presuppose everywhere defined operations. Recognizing this problem, David Gries, in ``The Science of Programming,`` 1981, introduced the concept of partial functions into the mainstream of program correctness and gave hints how his treatment of partial functions could be formalized. Still, however, existing theorem provers and software verifiers have difficulties in checking software with partial functions, because of absence of uniform first order treatment of partial functions within classical 2-valued logic. Several rigorous mechanisms that took partiality into account were introduced [Wirsing 1990, Breu 1991, VDM 1986, 1990, etc.]. However, they either did not discuss correctness proofs or departed from first order logic. To fill this gap, the authors provide a semantics for software correctness proofs with partial functions within classical 2-valued 1st order logic. They formalize the Gries treatment of partial functions and also cover computations of functions whose argument lists may be only partially available. An example is nuclear reactor control relying on sensors which may fail to deliver sense data. This approach is sufficiently general to cover correctness proofs in various implementation languages.

The US Department of Energy (DOE) uses sensitive or classified parts and material that must be protected and accounted for. The authors believe there is a need for an automated system that can help protect and monitor these parts and material. In response to this need Sandia National Laboratories (SNL) has developed a real-time personnel and material tracking system called PAMTRAK that has been installed at selected DOE facilities. PAMTRAK safeguards sensitive parts and material by tracking tags worn by personnel and by monitoring sensors attached to the parts or material. This paper describes the goals when designing PAMTRAK, the PAMTRAK system components, the current installations, and the benefits a site can expect when using PAMTRAK.

Planar, surface micromachined pressure sensors have been fabricated by an extension of the chemical-mechanical polishing (CMP) process. CMP eliminates many of the fabrication problems associated with the photolithography, dry etch, and metallization of non-planar devices. Furthermore, CMP adds additional design flexibility. The sensors are based upon deformable, silicon nitride diaphragms with polysilicon piezoresistors. Absolute pressure is detected by virtue of reference pressure cavities underneath the diaphragms. Process details are discussed and characteristics from many devices are presented.

Phosphorus diffusions are used in the fabrication process for nearly all crystalline-silicon (c-Si) photovoltaic solar cells to form the emitter of the solar cell. These phosphorous diffusions are also well known to have beneficial gettering benefits, i.e., deleterious metallic impurities are gettered from the bulk of the c-Si substrate into the phosphorous doped layer. In this study, we examined the effect of oxidations performed after the phosphorus diffusion. We were particularly interested in using the oxidation to passivate the surface of the phosphorus diffusion. Post-diffusion oxidations or moderate temperature steps in oxidizing ambients are also commonly found in commercial fabrication sequences of c-Si solar cells. we found that the bulk lifetime was degraded in Czochralski (Cz) silicon due to the post-diffusion oxidation unless there was a gettering agent present during the oxidation. Possible explanations for these results are presented at the end of the paper.

The dehydrogenation of propane and isobutane was studied in hydrogen permselective packed bed membrane reactors and conventional packed bed reactors. Two different types of developmental membranes were investigated: sol-gel derived silica-based membranes and a pure palladium thin film supported by a porous ceramic substrate. The palladium membranes deactivated and eventually failed when exposed to both isobutane and propane dehydrogenation temperatures above 773 K. Moderate improvements in propylene and isobutylene yields were obtained with the silica-based membrane reactors. An isobutylene yield of 48 mole percent was obtained at a liquid hourly space velocity (LHSV) of 1.8 and temperature of 798 K compared to a yield of 39 percent in a conventional reactor operated with the same flow rate. Similar improvements in propylene yield were obtained when the silica-based membranes were tested in propane dehydrogenation experiments. There was no significant difference in the reaction selectivities for the desired olefin products when the membrane and conventional reactors were operated with the scone LHSV However, for a constant value of the olefin yield, the membrane reactors had a higher reaction selectivity since the desired yield was achieved at a higher LHSV where there was less time for side products to form. Catalyst deactivation rates were generally greater in the membrane reactors, especially when the reactors were operated with high hydrogen removal rates at temperatures of 773 K and above.

GaAs Junction Field Effect Transistors (JFETs) have attracted renewed attention for low-power, low-voltage electronics. JFETs have a significant advantage over MESFETs for low-power operation due to their higher gate barrier to current flow resulting from p/n junction gate. This paper reports recent advances in an all ion implanted self-aligned GaAs JFET with a gate length down to 0.3 {mu}m. By employing shallopw SiF implants next to the gate, dielectric sidewall spacers, and 50 keV source and drain implants, JFETs with a f{sub t} up to 49 GHz with good pinchoff and subthreshold characteristics have been realized. In addition, the JFET benefits from the use of shallow Zn or Cd implantation to form abrupt p{sup +}/n gate profiles.

Deep etching of GaAs is a critical process step required for many device applications including fabrication of through-substrate via holes for monolithic microwave integrated circuits (MMICs). Use of high-density plasmas, including inductively coupled plasmas (ICP), offers an alternative approach to etching vias as compared to more conventional parallel plate reactive ion etch systems. This paper reports ICP etching of GaAs vias at etch rates of about 5.3 {mu}m/min with via profiles ranging from highly anistropic to conical.

In this study, 33 different conversion coatings were applied to 5 different Al alloy substrates. Salt spray exposure testing and EIS (electrochemical impedance spectroscopy) were conducted for comparison. A relation was developed.

With the demonstration of implant doping of GaN and the resulting need to perform the activation anneal at 1100 C, details of thermal stability of the GaN surface needs to be understood. This work reports on the use of a sputtered AlN encapsulant to preserve the surface of GaN during such annealing. The surface was characterized by formation of Pt/Au Schottky contacts and by AES. Schottky contacts deposited an GaN annealed wtih the AlN encapsulant displayed good rectification properties while those formed on GaN annealed uncapped approached ohmic behavior. AES analysis supports the hypothesis that the uncapped sample has lost N from the very near surface which creates N-vacancies that act as donors and thereby form an n{sup +}-surface layer.

This paper presents design, analysis, and first results of the high brightness electron beam experiments currently under investigation at Sandia. Anticipated beam parameters are: energy 12 MeV, current 35-40 kA, rms radius 0.5 mm, pulse duration 40 ns FWHM. The accelerator is SABRE, a pulsed LIVA modified to higher impedance, and the electron source is a magnetically immersed foilless electron diode. 20 to 30 Tesla solenoidal magnets are required to insulate the diode and contain the beam to its extremely small sized (1 mm) envelope. These experiments are designed to push the technology to produce the highest possible electron current in a submillimeter radius beam. Design, numercial simulations, and first experimental results are presented.

Development of a complementary heterostructure field effect transistor (CHFET) technology for low-power, mixed-mode digital-microwave applications is presented. Digital CHFET technology with independently optimizable transistors has been shown to operate with 319 ps loaded gate delays at 8.9 fJ. Power consumption is dominated by leakage currents of the p-channel FET, while performance is determined by the characteristics of 0.7 {mu}m gate length devices. As a microwave technology, the nJFET forms the basis of low-power cirucitry without any modification to the digital process. Narrow band amplification with a 0.7x100 {mu}m nJFET has been demonstrated at 2.1-2.4 GHz with gains of 8-10 dB at 1 mW power. These amplifiers showed a minimum noise figure of 2.5 dB. Next generation CHFET transistors with sub 0.5 {mu}m gate lengths have also been developed. Cutoff frequencies of 49 and 11.5 GHz were achieved for n- and p-channel FETs with 0.3 and 0.4 {mu}m gates, respectively. These FETs will enable enhancements in both digital and microwave circuits.

Extracting information from unstructured text has become an important research area in recent years due to the large amount of text now electronically available. This status report describes the findings and work done during the second year of a two-year Laboratory Directed Research and Development Project. Building on the first-year`s work of identifying important entities, this report details techniques used to group words into semantic categories and to output templates containing selective document content. Using word profiles and category clustering derived during a training run, the time-consuming knowledge-building task can be avoided. Though the output still lacks in completeness when compared to systems with domain-specific knowledge bases, the results do look promising. The two approaches are compatible and could complement each other within the same system. Domain-independent approaches retain appeal as a system that adapts and learns will soon outpace a system with any amount of a priori knowledge.

This paper provides an overview of the history and process of establishing a cooperative research and development agreement (CRADA) between Sandia National Laboratories and Magnavox Electronic Systems Company for the design, development, and testing of a 360-degree scanning, imaging, intrusion detection sensor. The subject of the CRADA is the Advanced Exterior Sensor (AES). It is intended for exterior use at ranges from 50 to 1,500 meters and uses a combination of three sensing technologies (infrared, visible, and radar) and a new data processing method to provide low false-alarm intrusion detection and tracking combined with immediate visual assessment. The establishment of this CRADA represents a new paradigm in the cooperation between the Department of Defense, the Department of Energy, the National Laboratories and Private Industry. Although a formal document has now been executed, a CRADA is, nonetheless, primarily an agreement to work with each other to achieve goals that might otherwise be unattainable.

In this paper, the authors explore the feasibility of using the distributed Bragg reflector, grown on the substrate for a VCSEL (Vertical Cavity Surface Emitting Laser), to provide waveguiding within the substrate. This waveguiding could serve as an interconnection among VCSELs in an array. Before determining the feasibility of waveguide interconnected VCSELs, two analysis methods are presented and evaluated for their applicability to this problem. The implementations in Mathematica of both these methods are included. Results of the analysis show that waveguiding in VCSEL structures is feasible. Some of the many possible uses of waveguide interconnected VCSELs are also briefly discussed. The tools and analysis presented in this report can be used to evaluate such system concepts and to do detailed design calculations.

This is the final report of a one-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). If a beam-chopping system could be developed for the Los Alamos Meson Physics Facility low-energy beam line, there would be potential to operate the Los Alamos Neutron Scattering Center (LANSCE) at much higher power and duty factor and enable such operation with a radio-frequency quadrapole (RFQ) injector. This would greatly extend the capability of the facility. To accommodate LANSCE operation in the new configuration, a chopped beam must be created in the low-energy transport line before the RFQ. Chopping in this region has never been demonstrated and constitutes the major uncertainty of the proposal and determines the critical path for project completion. This study produces a better understanding of the physics involved in chopping an H-beam in a dilute plasma background, and in transporting a chopped H-beam through a neutralized or partially neutralized plasma channel, as well as an estimate for the optimum neutralization strategy for the beam chopping and transport between the ion source and the RFQ.

This 1995 report contains data from routine radiological and non-radiological environmental monitoring activities. Summaries of significant environmental compliance programs in progress, such as National Environmental Policy Act documentation, environmental permits, environmental restoration and various waste management programs at Sandia National Laboratories in Albuquerque, New Mexico, are included.

DOE installations possessing sufficient quantities of fissile material to potentially constitute a critical mass, such that the excessive exposure of personnel to radiation from a nuclear accident is possible, are required to provide nuclear accident dosimetry services. This document describes the personal nuclear accident dosimeter (PNAD) used by SNL and prescribes methodologies to initially screen, and to process PNAD results. In addition, this report describes PNAD dosimetry results obtained during the Nuclear Accident Dosimeter Intercomparison Study (NAD23), held during 12-16 June 1995, at Los Alamos National Laboratories. Biases for reported neutron doses ranged from -6% to +36% with an average bias of +12%.

Ten commercial MgO powders were evaluated for their suitability to act as a binder in the separator of thermal batteries to immobilize the electrolyte when it is molten. One brand in particular, Maglite S from Calgon, outperformed all the others. This report describes the results of a characterization study of this MgO as well as similar materials from other commercial vendors. The study objective was to define the critical properties of Maglite S MgO that are responsible for its superior performance in thermal-battery separators. Separator mixes were prepared with the various MgO powders and the resulting powders and pellets were characterized, to correlate key physical properties of these materials to select physical and chemical properties of the MgO powders used in their preparation. The MgO pore-size distribution was the only parameter that could be related to the deformation and electrolyte-leakage behavior of separator pellets. A potential replacement for the Maglite S is currently being qualified, since Maglite S MgO is no longer available.

This report summarizes the environmental surveillance activities conducted by Sandia National Laboratories. Sandia National Laboratories` responsibility for environmental surveillance results extends to those activities performed by Sandia National Laboratories or under its direction. Results from other organizations environmental surveillance activities are included to provide a measure of completeness. Other environmental compliance programs such as the National Environmental Policy Act of 1969, environmental permits, and environmental restoration and waste management programs are also included in this report, prepared for the U.S. Department of Energy (DOE) in compliance with DOE Order 5400.1.

Capillary barriers consisting of tilted fine-over-coarse layers have been suggested as landfill covers as a means to divert water infiltration away from sensitive underground regions under unsaturated flow conditions, especially for arid and semi-arid regions. Typically, the HELP code is used to evaluate landfill cover performance and design. Unfortunately, due to its simplified treatment of unsaturated flow and its essentially one-dimensional nature, HELP is not adequate to treat the complex multidimensional unsaturated flow processes occurring in a tilted capillary barrier. In order to develop the necessary mechanistic code for the performance evaluation of tilted capillary barriers, an efficient and comprehensive unsaturated flow code needs to be selected for further use and modification. The present study evaluates a number of candidate mechanistic unsaturated flow codes for application to tilted capillary barriers. Factors considered included unsaturated flow modeling, inclusion of evapotranspiration, nodalization flexibility, ease of modification, and numerical efficiency. A number of unsaturated flow codes are available for use with different features and assumptions. The codes chosen for this evaluation are TOUGH2, FEHM, and SWMS{_}2D. All three codes chosen for this evaluation successfully simulated the capillary barrier problem chosen for the code comparison, although FEHM used a reduced grid. The numerical results are a strong function of the numerical weighting scheme. For the same weighting scheme, similar results were obtained from the various codes. Based on the CPU time of the various codes and the code capabilities, the TOUGH2 code has been selected as the appropriate code for tilted capillary barrier performance evaluation, possibly in conjunction with the infiltration, runoff, and evapotranspiration models of HELP. 44 refs.

An oblique penetration modeling procedure is evaluated by correlation with onboard acceleration data from a series of six penetration tests into Antelope Dry Lake soil at Tonopah Test Range, Nevada. The modeling represents both the loading which is coupled to the penetrator bending and the penetrator structure including connections between the major subsections. Model results show reasonable agreement with the data which validates the modeling procedure within a modest uncertainty related to accelerometer clipping and rattling of the telemetry package. The experimental and analytical results provide design guidance for the location and lateral restraint of components to reduce their shock environment.

CHEM{_}MACCS has been developed from the radiological accident consequence code, MACCS, to perform probabilistic calculations of potential off-site consequences of the accidental atmospheric release of hazardous chemicals. The principal phenomena considered in CHEM{_}MACCS are atmospheric transport, mitigative actions based on dose projection, dose accumulation by a number of pathways, and early and latent health effects. CHEM{_}MACCS provides the following capabilities: (1) statistical weather sampling data (8,760 hourly data points per year), (2) population dose and health effect risk calculations based on site-specific population data, (3) health effects calculations including the consideration of potential site specific mitigative actions (evacuation and shielding), and (4) modeling of multiple release segments. Three different sample problems are contained in this report to show how to use CHEM{_}MACCS. Three test problems are run to compare CHEM{_}MACCS and D2PC. The doses versus the downwind centerline distances from the source for the given doses are in very close agreement.

An autonomous gas chromatograph system was designed and built to support the Thermal Enhanced Vapor Extraction System (TEVES) demonstration. TEVES is a remediation demonstration that seeks to enhance an existing technology (vacuum extraction) by adding a new technology (soil heating). A pilot scale unit was set up at one of the organic waste disposal pits at the Sandia National Laboratories Chemical Waste Landfill (CWL) in Tech Area 3. The responsibility for engineering a major part of the process instrumentation for TEVES belonged to the Manufacturing Control Subsystems Department. The primary mission of the one-of-a-kind hardware/software system is to perform on-site gas sampling and analysis to quantify a variety of volatile organic compounds (VOCs) from various sources during TEVES operations. The secondary mission is to monitor a variety of TEVES process physical parameters such as extraction manifold temperature, pressure, humidity, and flow rate, and various subsurface pressures. The system began operation in September 1994 and was still in use on follow-on projects when this report was published.

In support of the Motorola CRADA, the capabilities of the computational fluid dynamics code FIDAP (Fluid Dynamics International) for simulating problems involving fluid flow, heat transport, and chemical reactions have been assessed and enhanced as needed for semiconductor-processing applications (e.g. chemical vapor deposition). A novel method of treating surface chemical species that uses only pre-existing FIDAP commands is described and illustrated with test problems. A full-Jacobian treatment of the chemical reaction rate expressions during formation of the stiffness matrix has been implemented in FIDAP for both the Arrhenius-parameter and user-subroutine methods of specifying chemical reactions, where the Jacobian terms can be calculated analytically or numerically. This formulation is needed to obtain convergence when reaction rates become large compared to transport rates (stiff chemistry). Several test problems are analyzed, and in all cases this approach yields good convergence behavior, even for extremely stiff fluid-phase and surface reactions. A stiff segregated algorithm has been developed and implemented in FIDAP. Analysis of test problems indicates that this algorithm yields improved convergence behavior compared with the original segregated algorithm. This improved behavior enables segregated techniques to be applied to problems with stiff chemistry, as required for large three-dimensional multi-species problems.

The current that flows between the electrical test equipment and the nuclear explosive must be limited to safe levels during electrical tests conducted on nuclear explosives at the DOE Pantex facility. The safest way to limit the current is to use batteries that can provide only acceptably low current into a short circuit; unfortunately this is not always possible. When it is not possible, current limiters, along with other design features, are used to limit the current. Three types of current limiters, the fuse blower, the resistor limiter, and the MOSFET-pass-transistor limiters, are used extensively in Pantex test equipment. Detailed failure mode and effects analyses were conducted on these limiters. Two other types of limiters were also analyzed. It was found that there is no best type of limiter that should be used in all applications. The fuse blower has advantages when many circuits must be monitored, a low insertion voltage drop is important, and size and weight must be kept low. However, this limiter has many failure modes that can lead to the loss of over current protection. The resistor limiter is simple and inexpensive, but is normally usable only on circuits for which the nominal current is less than a few tens of milliamperes. The MOSFET limiter can be used on high current circuits, but it has a number of single point failure modes that can lead to a loss of protective action. Because bad component placement or poor wire routing can defeat any limiter, placement and routing must be designed carefully and documented thoroughly.

Transient dynamics simulations are commonly used to model phenomena such as car crashes, underwater explosions, and the response of shipping containers to high-speed impacts. Physical objects in such a simulation are typically represented by Lagrangian meshes because the meshes can move and deform with the objects as they undergo stress. Fluids (gasoline, water) or fluid-like materials (earth) in the simulation can be modeled using the techniques of smoothed particle hydrodynamics. Implementing a hybrid mesh/particle model on a massively parallel computer poses several difficult challenges. One challenge is to simultaneously parallelize and load-balance both the mesh and particle portions of the computation. A second challenge is to efficiently detect the contacts that occur within the deforming mesh and between mesh elements and particles as the simulation proceeds. These contacts impart forces to the mesh elements and particles which must be computed at each timestep to accurately capture the physics of interest. In this paper we describe new parallel algorithms for smoothed particle hydrodynamics and contact detection which turn out to have several key features in common. Additionally, we describe how to join the new algorithms with traditional parallel finite element techniques to create an integrated particle/mesh transient dynamics simulation. Our approach to this problem differs from previous work in that we use three different parallel decompositions, a static one for the finite element analysis and dynamic ones for particles and for contact detection. We have implemented our ideas in a parallel version of the transient dynamics code PRONTO-3D and present results for the code running on a large Intel Paragon.

Sandia National Laboratories and Oceania, Inc. entered into a Cooperative Research and Development Agreement (CRADA) in November 1993 to provide ``Information Integrity and Privacy for Computerized Medical Patient Records`` (CRADA No. SC93/01183). The main objective of the project was to develop information protection methods that are appropriate for databases of patient records in health information systems. This document describes the findings and alternative solutions that resulted from this CRADA.

PHASER (Probabilistic Hybrid Analytical System Evaluation Routine) is a software tool that has the capability of incorporating subjective expert judgment into probabilistic safety analysis (PSA) along with conventional data inputs. An earlier report described the PHASER methodology, but only gave a cursory explanation about how dependence was incorporated in Version 1.10 and about how ``Importance`` and ``Sensitivity`` measures were to be incorporated in Version 2.00. A more detailed description is given in this report. The basic concepts involve scale factors and confidence factors that are associated with the stochastic variability and subjective uncertainty (which are common adjuncts used in PSA), and the safety risk extremes that are crucial to safety assessment. These are all utilized to illustrate methodology for incorporating dependence among analysis variables in generating PSA results, and for Importance and Sensitivity measures associated with the results that help point out where any major sources of safety concern arise and where any major sources of uncertainty reside, respectively.

This work addresses a new kind of silicon based chemical sensor that combines the reliability and stability of silicon microelectronic field effect devices with the highly selective and sensitive immunoassay. The sensor works on the principle that thin SiN layers on lightly doped Si can detect pH changes rapidly and reversibly. The pH changes affect the surface potential, and that can be quickly determined by pulsed photovoltage measurements. To detect other species, chemically sensitive films were deposited on the SiN where the presence of the chosen analyte results in pH changes through chemical reactions. A invention of a cell sorting device based on these principles is also described. A new method of immobilizing enzymes using Sandia`s sol-gel glasses is documented and biosensors based on the silicon wafer and an amperometric technique are detailed.

Microsensors do not have the selectivity to chemical species available in large laboratory instruments. This project employed arrays of catalytically gated silicon microsensors with different catalysts to create data streams which can be analyzed by pattern recognition programs. One of the most significant accomplishments of the program was the demonstration of that mixtures of H{sub 2} with the oxidants NO{sub x} and O{sub 2} could distinguished from one another by the use of different catalytic metals on the Sandia Robust Hydrogen (SRH) sensors and the newly developed pattern recognition algorithm. This sensor system could be used to identify explosive gas mixtures and analyze exhaust streams for pollution control.

This manual describes the use of MPSalsa, an unstructured finite element (FE) code for solving chemically reacting flow problems on massively parallel computers. MPSalsa has been written to enable the rigorous modeling of the complex geometry and physics found in engineering systems that exhibit coupled fluid flow, heat transfer, mass transfer, and detailed reactions. In addition, considerable effort has been made to ensure that the code makes efficient use of the computational resources of massively parallel (MP), distributed memory architectures in a way that is nearly transparent to the user. The result is the ability to simultaneously model both three-dimensional geometries and flow as well as detailed reaction chemistry in a timely manner on MT computers, an ability we believe to be unique. MPSalsa has been designed to allow the experienced researcher considerable flexibility in modeling a system. Any combination of the momentum equations, energy balance, and an arbitrary number of species mass balances can be solved. The physical and transport properties can be specified as constants, as functions, or taken from the Chemkin library and associated database. Any of the standard set of boundary conditions and source terms can be adapted by writing user functions, for which templates and examples exist.

The USW SD-7 drill hole is one of several holes drilled under Site Characterization Plan Study 8.3.1.4.3.1, also known as the Systematic Drilling Program, as part of the U.S. Department of Energy characterization program at Yucca Mountain, Nevada. The Yucca Mountain site has been proposed as the potential location of a repository for high-level nuclear waste. The SD-7 drill hole is located near the southern end of the potential repository area and immediately to the west of the Main Test Level drift of the Exploratory Studies Facility. The hole is not far from the junction of the Main Test Level drift and the proposed South Ramp decline. Drill hole USW SD-7 is 2675.1 ft (815.3 m) deep, and the core recovered nearly complete sections of ash-flow tuffs belonging to the lower half of the Tiva Canyon Tuff, the Pah Canyon Tuff, and the Topopah Spring Tuff, all of which are part of the Miocene Paintbrush Group. Core was recovered from much of the underlying Calico Hills Formation, and core was virtually continuous in the Prow Pass Tuff and the Bullfrog Tuff. The SD-7 drill hole penetrated the top several tens of feet into the Tram Tuff, which underlies the Prow Pass and Bullfrog Tuffs. These latter three units are all formations of the Crater Flat Group, The drill hole was collared in welded materials assigned to the crystal-poor middle nonlithophysal zone of the Tiva Canyon Tuff; approximately 280 ft (85 m) of this ash-flow sheet was penetrated by the hole. The Yucca Mountain Tuff appears to be missing from the section at the USW SD-7 location, and the Pah Canyon Tuff is only 14.5 ft thick. The Pah Canyon Tuff was not recovered in core because of drilling difficulties, suggesting that the unit is entirely nonwelded. The presence of this unit is inferred through interpretation of down-hole geophysical logs.

This LDRD program No. 3505.230 explored a new approach to monolithic integration of active waveguides and rare-earth solid state lasers directly onto III-V substrates. It involved selectively incorporating rare-earth ions into spin-on glasses (SOGs) that could be solvent cast and then patterned with conventional microelectronic processing. The patterned, rare-earth spin-on glasses (RESOGs) were to be photopumped by laser diodes prefabricated on the wafer and would serve as directly integrated active waveguides and/or rare-earth solid state lasers.

Under the auspices of the {open_quotes}Government/Industry Wind Technology Applications Project{close_quotes} [{open_quotes}Letter of Interest{close_quotes} (LOI) Number RC-1-11101], Flo Wind Corp. has successfully developed, tested, and delivered a high-energy rotor upgrade candidate for their 19-meter Vertical Axis Wind Turbine. The project included the demonstration of the innovative extended height-to-diameter ratio concept, the development of a continuous span single-piece composite blade, the demonstration of a continuous blade manufacturing technique, the utilization of the Sandia National Laboratories developed SNLA 2150 natural laminar flow airfoil and the reuse of existing wind turbine and wind power plant infrastructure.

PHASER (Probabilistic Hybrid Analytical System Evaluation Routine) is a computer code for solving the top event probability of a system fault tree. It has the capability for easy migration of individual basic event probabilities from a zero-{open_quotes}scale{close_quotes}-factor (completely subjective) state to one in which the analyst has total knowledge (completely stochastic) about each basic event. The code implements a fuzzy algebra solution for subjective data, a probabilistic solution for stochastic data, and a hybrid mathematics solution for data that are partly subjective and partly stochastic. Events that are not completely subjective or completely stochastic are hybrid events and are internally handled as such. The stochastic and fuzzy ranges of uncertainty in the top event probability are also computed for the analyst. These are provided in the form of a fuzzy function for the subjective uncertainty, a probability density function (PDF) for the stochastic variability, and the overall {open_quotes}confidence{close_quotes} factors for the two constituents of uncertainty, giving a complete hybrid result. PHASER interfaces with other Sandia codes (SABLE, LHS and LHSPOST) to assist the user in determining cutsets, and to compute probability density functions.

In 1983, high-level radioactive waste repository performance requirements related to groundwater travel time were defined by NRC subsystem regulation 10 CFR 60.113. Although DOE is not presently attempting to demonstrate compliance with that regulation, understanding of the prevalence of fast paths in the groundwater flow system remains a critical element of any safety analyses for a potential repository system at Yucca Mountain, Nevada. Therefore, this analysis was performed to allow comparison of fast-path flow against the criteria set forth in the regulation. Models developed to describe the conditions for initiation, propagation, and sustainability of rapid groundwater movement in both the unsaturated and saturated zones will form part of the technical basis for total- system analyses to assess site viability and site licensability. One of the most significant findings is that the fastest travel times in both unsaturated and saturated zones are in the southern portion of the potential repository, so it is recommended that site characterization studies concentrate on this area. Results support the assumptions regarding the importance of an appropriate conceptual model of groundwater flow and the incorporation of heterogeneous material properties into the analyses. Groundwater travel times are sensitive to variation/uncertainty in hydrologic parameters and in infiltration flux at upper boundary of the problem domain. Simulated travel times are also sensitive to poorly constrained parameters of the interaction between flow in fractures and in the matrix.

This report provides a users` guide for parallel processing ITS on a UNIX workstation network, a shared-memory multiprocessor or a massively-parallel processor. The parallelized version of ITS is based on a master/slave model with message passing. Parallel issues such as random number generation, load balancing, and communication software are briefly discussed. Timing results for example problems are presented for demonstration purposes.

The work is part of the rock mechanics effort for the Yucca Mountain Site Characterization Program. The laboratory-scale experiments are intended to provide high quality data on the mechanical behavior of jointed structures that can be used to validate complex numerical models for rock-mass behavior. Frictional sliding between simulated rock joints was studied using phase shifting moire interferometry. A model, constructed from stacks of machined and sandblasted granite plates, contained a central hole bore normal to the place so that frictional slip would be induced between the plates near the hole under compressive loading. Results show a clear evolution of slip with increasing load. Since the rock was not cycled through loading- unloading, the quantitative differences between the three data sets are probably due to a ``wearing-in`` effect. The highly variable spatial frequency of the data is probably due to the large grain size of the granite and the stochastic frictional processes. An unusual feature of the evolution of slip with increasing load is that as the load gets larger, some plates seem to return to a null position. Figs, 6 refs.

The Department of Energy is investigating Yucca Mountain, Nevada as a potential site for commercial radioactive waste disposal in a mined geologic repository. One critical aspect of site suitability is the tectonic stability of the repository site. The levels of risk from both actual fault displacements in the repository block and ground shaking from nearby earthquakes are being examined. In particular, it is necessary to determine the expected level of ground shaking at the repository depth for large seismic sources such as nearby large earthquakes or underground nuclear explosions (UNEs). Earthquakes are expected to cause the largest ground motions at the site, however, only underground nuclear explosion data have been obtained at the repository depth level (about 350m below the ground level) to date. In this study we investigate ground motion from Nevada Test Site underground nuclear explosions recorded at Yucca Mountain to establish a compressional velocity model for the uppermost 350m of the mountain. This model is useful for prediction of repository-level ground motions for potential large nearby earthquakes.

The Do-It-Now (DIN) building maintenance system is proposed to reduce the cost of routine building maintenance and repairs and to improve customer satisfaction with maintenance services. DIN uses a team approach to periodically inspect buildings and provide maintenance services on the spot. It emphasizes communications between the customers and the craftspeople performing the work. The system was designed using a reengineering approach that characterized the existing maintenance work control system, analyzed comparable systems in other DOE laboratories, envisioned an ideal system, and proposed a workable, testable system for initial implementation. At each stage, input was solicited from customer representatives and Facilities management to ensure meeting customer requirements with an implementable system.

This report describes the analysis and conclusion of an investigation of the carbon monoxide emissions resulting from Sandia National Laboratories and Department of Energy (DOE) commuter and on-base traffic for the Clean Air Act (CAA) Conformity Determination. Albuquerque/Bernalillo County was classified as a nonattainment area by the Environmental Protection Agency. Nonattainment area is an area which is shown by monitored data or which is calculated by air quality modeling to exceed any National Ambient Air Quality Standard (NAAQS) for the pollutant. Albuquerque/Bernalillo County exceeds the NAAQS for carbon monoxide and ozone. The Conformity Determination was needed to complete the CAA Title V Permitting process for SNL and the DOE. The analysis used the EPA approved MOBILE5a Carbon Monoxide (CO) emissions modeling program. This analysis will provide a baseline for mobile sources to allow Sandia to estimate any future activity and how that activity will impact CO emissions. The General Conformity Rule (AQCR 43) requires that operations which will increase CO emissions in nonattaimnent or maintenance areas such as Bernalillo County undergo conformity analyses to determine whether or not they will impact ambient air quality in the area.

The US DOE has amassed over 555,000 metric tons of depleted uranium from its uranium enrichment operations. Rather than dispose of this depleted uranium as waste, this study explores a beneficial use of depleted uranium as metal shielding in casks designed to contain canisters of vitrified high-level waste. Two high-level waste storage, transport, and disposal shielded cask systems are analyzed. The first system employs a shielded storage and disposal cask having a separate reusable transportation overpack. The second system employs a shielded combined storage, transport, and disposal cask. Conceptual cask designs that hold 1, 3, 4 and 7 high-level waste canisters are described for both systems. In all cases, cask design feasibility was established and analyses indicate that these casks meet applicable thermal, structural, shielding, and contact-handled requirements. Depleted uranium metal casting, fabrication, environmental, and radiation compatibility considerations are discussed and found to pose no serious implementation problems. About one-fourth of the depleted uranium inventory would be used to produce the casks required to store and dispose of the nearly 15,400 high-level waste canisters that would be produced. This study estimates the total-system cost for the preferred 7-canister storage and disposal configuration having a separate transportation overpack would be $6.3 billion. When credits are taken for depleted uranium disposal cost, a cost that would be avoided if depleted uranium were used as cask shielding material rather than disposed of as waste, total system net costs are between $3.8 billion and $5.5 billion.

An integral part of the licensing procedure for the potential nuclear waste repository at Yucca Mountain, Nevada, involves characterization of the in situ rheology for the design and construction of the facility and the emplacement of canisters containing radioactive waste. The data used to model the thermal and mechanical behavior of the repository and surrounding lithologies include dry and saturated bulk densities, average grain density, porosity, compressional and shear wave velocities, elastic moduli, and compressional and tensional fracture strengths. In this study, a suite of experiments was performed on cores recovered from boreholes UE25 NRG-2, 2A, 2B, and 3 drilled in support of the Exploratory Studies Facility (ESF) at Yucca Mountain. The holes penetrated the Timber Mountain tuff and two thermal/mechanical units of the Paintbrush tuff. The thermal/mechanical stratigraphy was defined by Ortiz to group rock horizons of similar properties for the purpose of simplifying modeling efforts. The relationship between the geologic stratigraphy and the thermal/mechanical stratigraphy for each borehole is presented. The tuff samples in this study have a wide range of welding characteristics (usually reflected in sample porosity), and a smaller range of mineralogy and petrology characteristics. Generally, the samples are silicic, ash-fall tuffs that exhibit large variability in their elastic and strength properties.

Experimental results are presented for bulk and mechanical properties measurements on specimens of the Paintbrush tuff recovered from boreholes UE25 NRG-4 and -5, at Yucca Mountain, Nevada. Measurements have been performed on three thermal/mechanical units, PTn, TSwl, and TSw2. On each specimen the following bulk properties have been reported: dry bulk density, saturated bulk density, average grain density, and porosity. Unconfined compression to failure, confined compression to failure, and indirect tensile strength tests were performed on selected specimens recovered from the boreholes. In addition, compressional and shear wave velocities were measured on specimens designated for unconfined compression and confined compression experiments. Measurements were conducted at room temperature on nominally water-saturated specimens. The nominal rate for the fracture experiments was 10{sup -5}s{sup -1}.

Sandia National Laboratories, New Mexico (SNL/NM) is managed and operated by Sandia Corporation, a Lockheed Martin Company. SNL/NM is located on land owned by the U.S. Department of Energy (DOE) within the boundaries of the Kirtland Air Force Base (KAFB) in Albuquerque, New Mexico. The major responsibilities of SNL/NM are the support of national security and energy projects. Low-level radioactive waste (LLW) is generated by some of the activities performed at SNL/NM in support of the DOE. This report describes potential environmental effects of the shipments of low-level radioactive wastes to other sites.

Ambient air samples were taken at two locations in the East Mountain Area in conjunction with thermal testing at the Lurance Canyon Burn Site (LCBS). The samples were taken to provide measurements of particulate matter with a diameter less than or equal to 10 micrometers (PM{sub 10}) and volatile organic compounds (VOCs). This report summarizes the results of the sampling performed in 1995. The results from small-scale testing performed to determine the potentially produced air pollutants in the thermal tests are included in this report. Analytical results indicate few samples produced measurable concentrations of pollutants believed to be produced by thermal testing. Recommendations for future air sampling in the East Mountain Area are also noted.

One of the most common oil-field treatments is hot oiling to remove paraffin from wells. Even though the practice is common, the thermal effectiveness of the process is not commonly understood. In order for producers to easily understand the thermodynamics of hot oiling, a simple tool is needed for estimating downhole temperatures. Such a tool has been developed that was distributed as a compiled, public-domain-software spreadsheet. That spreadsheet has evolved into an interactive from on the World Wide Web and has been adapted into a Windows{trademark} program by Petrolite, St. Louis MO. The development of such a tools was facilitated by expressing downhole temperatures in terms of analytic formulas. Considerable algebraic work is required to develop such formulas. Also, the data describing hot oiling is customarily a mixture of practical units that must be converted to a consistent set of units. To facilitate the algebraic manipulations and to assure unit conversions are correct, during development parallel calculations were made using the spreadsheet and a symbolic mathematics program. Derivation of the formulas considered falling film flow in the annulus and started from the transient differential equations so that the effects of the heat capacity of the tubing and casing could be included. While this approach to developing a software product does not have the power and sophistication of a finite element or difference code, it produces a user friendly product that implements the equations solved with a minimum potential for bugs. This allows emphasis in development of the product to be placed on the physics.

This paper considers the protein structure prediction problem for lattice and off-lattice protein folding models that explicitly represent side chains. Lattice models of proteins have proven extremely useful tools for reasoning about protein folding in unrestricted continuous space through analogy. This paper provides the first illustration of how rigorous algorithmic analyses of lattice models can lead to rigorous algorithmic analyses of off-lattice models. The authors consider two side chain models: a lattice model that generalizes the HP model (Dill 85) to explicitly represent side chains on the cubic lattice, and a new off-lattice model, the HP Tangent Spheres Side Chain model (HP-TSSC), that generalizes this model further by representing the backbone and side chains of proteins with tangent spheres. They describe algorithms for both of these models with mathematically guaranteed error bounds. In particular, the authors describe a linear time performance guaranteed approximation algorithm for the HP side chain model that constructs conformations whose energy is better than 865 of optimal in a face centered cubic lattice, and they demonstrate how this provides a 70% performance guarantee for the HP-TSSC model. This is the first algorithm in the literature for off-lattice protein structure prediction that has a rigorous performance guarantee. The analysis of the HP-TSSC model builds off of the work of Dancik and Hannenhalli who have developed a 16/30 approximation algorithm for the HP model on the hexagonal close packed lattice. Further, the analysis provides a mathematical methodology for transferring performance guarantees on lattices to off-lattice models. These results partially answer the open question of Karplus et al. concerning the complexity of protein folding models that include side chains.