Low residual stress silicon oxynitride thin films are investigated for use as a replacement for silicon dioxide (SiO2) as sacrificial layer in surface micromachined microelectrical-mechanical systems (MEMS). It is observed that the level of residual stress in oxynitrides is a function of the nitrogen content in the film. MEMS film stacks are prepared using both SiO2 and oxynitride sacrificial layers. Wafer bow measurements indicate that wafers processed with oxynitride release layers are significantly flatter. Polycrystalline Si (poly-Si) cantilevers fabricated under the same conditions are observed to be flatter when processed with oxynitride rather than SiO2 sacrificial layers. These results are attributed to the lower post-processing residual stress of oxynitride compared to SiO2.
Two major problems associated with Si-based MEMS (MicroElectroMechanical Systems) devices are stiction and wear. Surface modifications are needed to reduce both adhesion and friction in micromechanical structures to solve these problems. In this paper, we will present a CVD (Chemical Vapor Deposition) process that selectively coats MEMS devices with tungsten and significantly enhances device durability. Tungsten CVD is used in the integrated-circuit industry, which makes this approach manufacturable. This selective deposition process results in a very conformal coating and can potentially address both stiction and wear problems confronting MEMS processing. The selective deposition of tungsten is accomplished through the silicon reduction of WF6. The self-limiting nature of this selective. We deposition process ensures the consistency necessary for process control. The tungsten is deposited after the removal of the sacrificial oxides to minimize stress and process integration problems. Tungsten coating adheres well and is hard and conducting, requirements for device performance. Furthermore, since the deposited tungsten infiltrates under adhered silicon parts and the volume of W deposited is less than the amount of Si consumed, it appears to be possible to release stuck parts that are contacted over small areas such as dimples. The wear resistance of selectively coated W parts has been shown to be significantly improved on microengine test structures.
A new type of surface micromachined ratcheting actuation system has been developed at the Microelectronics Development Laboratory at Sandia National Laboratories. The actuator uses a torsional electrostatic comb drive that is coupled to an external ring gear through a ratcheting scheme. The actuator can be operated with a single square wave, has minimal rubbing surfaces, maximizes comb finger density, and can be used for open-loop position control. The prototypes function as intended with a minimum demonstrated operating voltage of 18V. The equations of motion are developed for the torsional electrostatic comb drive. The resonant frequency, voltage vs. displacement and force delivery characteristics are predicted and compared with the fabricated device's performance.
An increasing number of applications are requiring fiber transmission of high-intensity laser pulses. Our particular interests have led us to examine carefully the fiber transmission of Q-switched pulses from multimode Nd:YAG lasers at their fundamental wavelength. The maximum pulse energy that can be transmitted through a particular fiber is limited by the onset of laser-induced breakdown and damage mechanisms. Laser breakdown at the fiber entrance face is often the first limiting process to be encountered, but other mechanisms can result in catastrophic damage at either fiber face, within the initial `entry' segment of the fiber, and at other internal sites along the fiber path. In the course of our studies we have examined a number of factors that govern the relative importance of different mechanisms, including laser characteristics, the design and alignment of injection optics, fiber end-face preparation, and fiber routing. The present study emphasizes the important criteria for injection optics in high-intensity fiber transmission, and illustrates the opportunities that now exist for innovative designs of optics to meet these criteria. Our consideration of diffractive optics to achieve desired injection criteria began in 1993, and we have evaluated a progression of designs since that time. In the present study, two recent designs for injection optics are compared by testing a sufficient number of fibers with each design to establish statistics for the onset of laser-induced breakdown and damage. In this testing we attempted to hold constant other factors that can influence damage statistics. Both designs performed well, although one was less successful in meeting all injection criteria and consequently showed a susceptibility to a particular damage process.
A transient, multidimensional model has been developed to simulate proton exchange membrane fuel cells. The model accounts simultaneously for electrochemical kinetics, current distribution, hydrodynamics, and multicomponent transport. A single set of conservation equations valid for flow channels, gas-diffusion electrodes, catalyst layers, and the membrane region are developed and numerically solved using a finite-volume-based computational fluid dynamics technique. The numerical model is validated against published experimental data with good agreement. Subsequently, the model is applied to explore hydrogen dilution effects in the anode feed. The predicted polarization curves under hydrogen dilution conditions are in qualitative agreement with recent experiments reported in the literature. The detailed two-dimensional electrochemical and flow/transport simulations further reveal that in the presence of hydrogen dilution in the fuel stream, hydrogen is depleted at the reaction surface, resulting in substantial anode mass transport polarization and hence a lower current density that is limited by hydrogen transport from the fuel stream to the reaction site. Finally, a transient simulation of the cell current density response to a step change in cell voltage is reported.
Correctly identifying the possible alteration products and accurately predicting their occurrence in a repository-relevant environment are the key for source-term calculations in a repository performance assessment. Uraninite in uranium deposits has long been used as a natural analog to spent fuel in a repository because of their chemical and structural similarity. In this paper, a SEM/AEM investigation has been conducted on a partially alterated uraninite sample from a uranium ore deposit of Shinkolobwe of Congo. The mineral formation sequences were identified: uraninite→uranyl hydrates→uranyl silicates→Ca-uranyl silicates or uraninite→uranyl silicates→Ca-uranyl silicates. Reaction-path calculations were conducted for the oxidative dissolution of spent fuel in a representative Yucca Mountain groundwater. The predicted sequence is in general consistent with the SEM observations. The calculations also show that uranium carbonate minerals are unlikely to become major solubility-controlling mineral phases in a Yucca Mountain environment. Some discrepancies between model predictions and field observations are observed. Those discrepancies may result from poorly constrained thermodynamic data for uranyl silicate minerals.
Crystalline phases of pyrochlore (e.g., CaPuTi2O7, CaUTi2O7) have been proposed as a durable ceramic waste form for disposal of high level radioactive wastes including surplus weapons-usable plutonium. In this paper, we use a linear free energy relationship to predict the Gibbs free energies of formation of pyrochlore phases (CaMTi2O7). The Pu-pyrochlore phase is predicted to be stable with respect to PuO2, CaTiO3, and TiO2 at room temperatures. Pu-pyrochlore is expected to be stable in a geologic repository where silica and carbonate components are absent or limited. We suggest that a repository in a salt formation be an ideal environment for disposal of high level, pyrochlore-based ceramic wastes. In such environment, adding CaO as a backfill will make pyrochlore minerals thermodynamically stable and therefore effectively prevents actinide release from these mineral phases.
Interfaces play an important role in determining the effect of electric fields on the mechanism of the formation of spinel by solid-state reaction. The reaction occurs by the movement of phase boundaries but the rate of this movement can be affected by grain boundaries in the reactants or in the reaction product. Only by understanding these relationships will it be possible to engineer their behavior. As a particular example of such a study, MgIn2O4 can be formed by the reaction between single-crystal MgO substrate and a thin film of In2O3 with or without an applied electric field. High-resolution backscattered electron (BSE) imaging and electron backscattered diffraction (EBSD) in a scanning electron microscope (SEM) has been used to obtain complementary chemical and crystallographic information.
A new class of semiconductor lasers that can potentially produce much more short pulse energy is presented. This new laser is not limited in volume or aspect ratio by the depth of a p-n junction and are created from current filaments in semi-insulating GaAs. A current filament semiconductor lasers (CFSL) that have produced 75 nJ of 890 nm radiation in 1.5 ns were tested. A filaments as long as 3.4 cm and several hundred microns in diameter in high gain GaAs photoconductive switches were observed. Their smallest dimension can be more than 100 times the carrier diffusion length in GaAs. The spectral narrowing, lasing thresholds, beam divergence, temporal narrowing and energies which imply lasing for several configurations of CFSL are reported.
Failure analysis (FA) tools have been applied to analyze tungsten coated polysilicon microengines. These devices were stressed under accelerated conditions at ambient temperatures and pressure. Preliminary results illustrating the failure modes of microengines operated under variable humidity and ultra-high drive frequency will also be shown. Analysis of tungsten coated microengines revealed the absence of wear debris in microengines operated under ambient conditions. Plan view imaging of these microengines using scanning electron microscopy (SEM) revealed no accumulation of wear debris on the surface of the gears or ground plane on microengines operated under standard laboratory conditions. Friction bearing surfaces were exposed and analyzed using the focused ion beam (FIB). These cross sections revealed no accumulation of debris along friction bearing surfaces. By using transmission electron microscopy (TEM) in conjunction with electron energy loss spectroscopy (EELS), we were able to identify the thickness, elemental analysis, and crystallographic properties of tungsten coated MEMS devices. Atomic force microscopy was also utilized to analyze the surface roughness of friction bearing surfaces.
One method of providing the mode selectivity necessary to insure single mode operation in a large diameter VCSEL is to independently control the size of the gain region and that of the optical mode. Numerical simulations quantity this approach by predicting lateral mode discrimination for different sized gain apertures. Calculations are experimentally confirmed by the fabrication and testing of 850 nm VCSELs employing hybrid ion implantation/selective oxidation that produce a single-mode output of more than 5 mW.
In this paper we concern ourselves with modified versions of the traditional brachistochrone and tautochrone problems. In the modified version of each problem the constant gravity model is replaced with an attractive inverse square law, consequently we name these the 1/r2 brachistochrone and 1/r2 tautochrone problems. With regard to the 1/r2 brachistochrone problem, we show that the shape of the minimizing curve is formally constructed from an infinite series of elliptic integrals, and we use a numerical optimal control technique to generate the trajectories. The 1/r2 tautochrone problem is solved using fractional calculus techniques and we show that the solution satisfies Lagrange's rule for tautochronous curves.
Electrostatic discharge (ESD) and electrical overstress (EOS) damage of Micro-Electro-Mechanical Systems (MEMS) has been identified as a new failure mode. This failure mode has not been previously recognized or addressed primarily due to the mechanical nature and functionality of these systems, as well as the physical failure signature that resembles stiction. Because many MEMS devices function by electrostatic actuation, the possibility of these devices not only being susceptible to ESD or EOS damage but also having a high probability of suffering catastrophic failure due to ESD or EOS is very real. Results from previous experiments have shown stationary comb fingers adhered to the ground plane on MEMS devices tested in shock, vibration, and benign environments. Using Sandia polysilicon microengines, we have conducted tests to establish and explain the ESD/EOS failure mechanism of MEMS devices. These devices were electronically and optically inspected prior to and after ESD and EOS testing. This paper will address the issues surrounding MEMS susceptibility to ESD and EOS damage as well as describe the experimental method and results found from ESD and EOS testing. The tests were conducted using conventional IC failure analysis and reliability assessment characterization tools. In this paper we will also present a thermal model to accurately depict the heat exchange between an electrostatic comb finger and the ground plane during an ESD event.
In 1993, the Government Performance and Results Act (GPRA, PL 103-62) was enacted. GPRA, which applies to all federal programs, has three components: strategic plans, annual performance plans, and metrics to show how well annual plans are being followed. As part of meeting the GRPA requirement in FY2000, a 14-member external peer review panel (the Garwin Committee) was convened on May 17-19, 2000 to review Sandia National Laboratories' Pulsed Power Programs as a component of the Performance Appraisal Process negotiated with the Department of Energy (DOE). The scope of the review included activities in inertial confinement fission (ICF), weapon physics, development of radiation sources for weapons effects simulation, x-ray radiography, basic research in high energy density physics (HEDP), and pulsed power technology research and development. In his charge to the committee, Jeffrey Quintenz, Director of Pulsed Power Sciences (1600) asked that the review be based on four criteria (1) quality of science, technology, and engineering, (2) programmatic performance, management, and planning, (3) relevance to national needs and agency missions, and (4) performance in the operation and construction of major research facilities. In addition, specific programmatic questions were posed by the director and by the DOE-Defense Programs (DP). The accompanying report, produced as a SAND document, is the report of the committee's findings.
This paper identifies active research topics concerning human machine interfaces for intelligent machine systems. The paper was compiled by performing a series of literature searches and organizing the information according to the author's interest in better directing his own Human Machine Interface (HMI) research. Introductory literature from outside the HMI communities is also referenced to provide context.
This report summarizes the results of work conducted by the Semiconductor Manufacturing Technologies Program at Sandia National Laboratories (Sandia) during 1999. This work was performed by one working group: the Semiconductor Equipment Technology Center (SETEC). The group's projects included Numerical/Experimental Characterization of the Growth of Single-Crystal Calcium Fluoride (CaF{sub 2}); The Use of High-Resolution Transmission Electron Microscopy (HRTEM) Imaging for Certifying Critical-Dimension Reference Materials Fabricated with Silicon Micromachining; Assembly Test Chip for Flip Chip on Board; Plasma Mechanism Validation: Modeling and Experimentation; and Model-Based Reduction of Contamination in Gate-Quality Nitride Reactor. During 1999, all projects focused on meeting customer needs in a timely manner and ensuring that projects were aligned with the goals of the National Technology Roadmap for Semiconductors sponsored by the Semiconductor Industry Association and with Sandia's defense mission. This report also provides a short history of the Sandia/SEMATECH relationship and a brief on all projects completed during the seven years of the program.
Thermally-Induced Voltage Alteration (TIVA) is a relatively new technique for locating electrical defects in integrated circuits [1,2]. This paper describes a novel application of TIVA, to locate design anomalies. A newly designed integrated circuit with high and inconsistent Quiescent Power Supply Current (IDDQ) was initially diagnosed with limited success using various failsite isolation techniques. The TIVA technique was successful in accurately locating design anomalies. Results from TIVA identified a spurious ring oscillator in the design. Design modifications carried out using a focussed ion beam (FIB), verified the accuracy of the results from TIVA. This study clearly extends the use of TIVA beyond that of locating electrical defects and anomalies into the realm of design debugging.
LM111 voltage comparators exhibit a wide range of total-dose-induced degradation. Simulations show this variability may be a natural consequence of the low base doping of the substrate PNP (SPNP) input transistors. Low base doping increases the SPNPs collector to base breakdown voltage, current gain, and densities. The build-up of oxide trapped charge (N OT) and interface traps (N IT) is shown to be a function of pre-irradiation bakes. Experimental data indicate that, despite its structural similarities to the LM111, irradiated input transistors of the LM124 operational amplifier do not exhibit the same sensitivity to variations in pre-irradiation thermal cycles. Further disparities in LM111 and LM124 responses may result from a difference in the oxide defect build-up in the two part types. Variations in processing, packaging, and circuit effects are suggested as potential explanations.
The worst case bias during total dose irradiation of partially depleted SOI transistors from two technologies is correlated to the device architecture. Experiments and simulations are used to analyze SOI back transistor threshold voltage shift and charge trapping in the buried oxide.
This report describes the features of Aspen-EE (Electricity Enhancement), a new model for simulating the interdependent effects of market decisions and disruptions in the electric power system on other critical infrastructures in the US economy. Aspen-EE extends and modifies the capabilities of Aspen, an agent-based model previously developed by Sandia National Laboratories. Aspen-EE was tested on a series of scenarios in which the rules governing electric power trades were changed. Analysis of the scenario results indicates that the power generation company agents will adjust the quantity of power bid into each market as a function of the market rules. Results indicate that when two power markets are faced with identical economic circumstances, the traditionally higher-priced market sees its market clearing price decline, while the traditionally lower-priced market sees a relative increase in market clearing price. These results indicate that Aspen-EE is predicting power market trends that are consistent with expected economic behavior.
A record high fundamental-mode power of 5.1 mW was achieved from coupled-resonator vertical-cavity lasers (CRVCLs). In conventional VCSELs, the extent to which the gain volume may be increased is limited by the onset of multi-mode operation. Results indicate that this limitation is circumvented in a coupled-resonator device allowing high power fundamental-mode operation.
This report discusses an automatic target recognition (ATR) algorithm for synthetic aperture radar (SAR) imagery that is based on coherent change detection techniques. The algorithm relies on templates created from training data to identify targets. Objects are identified or rejected as targets by comparing their SAR signatures with templates using the same complex correlation scheme developed for coherent change detection. Preliminary results are presented in addition to future recommendations.
Current approaches to reliability are not adequate to keep pace with the need for faster, better and cheaper products and systems. This is especially true in high consequence of failure applications. The original proposal for the LDRD was to look at this challenge and see if there was a new paradigm that could make reliability predictions, along with a quantitative estimate of the risk in that prediction, in a way that was faster, better and cheaper. Such an approach would be based on the underlying science models that are the backbone of reliability predictions. The new paradigm would be implemented in two software tools: the Virtual Reliability Realization System (VRRS) and the Reliability Expert System (REX). The three-year LDRD was funded at a reduced level for the first year ($120K vs. $250K) and not renewed. Because of the reduced funding, we concentrated on the initial development of the expertise system. We developed an interactive semiconductor calculation tool needed for reliability analyses. We also were able to generate a basic functional system using Microsoft Siteserver Commerce Edition and Microsoft Sequel Server. The base system has the capability to store Office documents from multiple authors, and has the ability to track and charge for usage. The full outline of the knowledge model has been incorporated as well as examples of various types of content.
>This report describes the innovative modeling approach developed as a result of a 3-year Laboratory Directed Research and Development project. The overall goal of this project was to provide an effective suite of solvers for advanced production planning at facilities in the nuclear weapons complex (NWC). We focused our development activities on problems related to operations at the DOE's Pantex Plant. These types of scheduling problems appear in many contexts other than Pantex--both within the NWC (e.g., Neutron Generators) and in other commercial manufacturing settings. We successfully developed an innovative and effective solution strategy for these types of problems. We have tested this approach on actual data from Pantex, and from Org. 14000 (Neutron Generator production). This report focuses on the mathematical representation of the modeling approach and presents three representative studies using Pantex data. Results associated with the Neutron Generator facility will be published in a subsequent SAND report. The approach to task-based scheduling described here represents a significant addition to the literature for large-scale, realistic scheduling problems in a variety of production settings.
This report describes the design of PICO, a C++ framework for implementing general parallel branch-and-bound algorithms. The PICO framework provides a mechanism for the efficient implementation of a wide range of branch-and-bound methods on an equally wide range of parallel computing platforms. We first discuss the basic architecture of PICO, including the application class hierarchy and the package's serial and parallel layers. We next describe the design of the serial layer, and its central notion of manipulating subproblem states. Then, we discuss the design of the parallel layer, which includes flexible processor clustering and communication rates, various load balancing mechanisms, and a non-preemptive task scheduler running on each processor. We describe the application of the package to a branch-and-bound method for mixed integer programming, along with computational results on the ASCI Red massively parallel computer. Finally we describe the application of the branch-and-bound mixed-integer programming code to a resource constrained project scheduling problem for Pantex.
The Air Force's Electronic Systems Center has funded Sandia National Laboratories to develop an Automatic Target Recognition (ATR) System for the Air Force's Joint STARS platform using Mercury Computer systems hardware. This report provides general theory on the internal operations of the Real-Time ATR system and provides some basic techniques that can be used to reconfigure the system and monitor its runtime operation. In addition, general information on how to interface an image formation processor and a human machine interface to the ATR is provided. This report is not meant to be a tutorial on the ATR algorithms.
A pointing control system is developed and tested for a flying gimbaled telescope. The two-axis pointing system is capable of sub-microradian pointing stability and high accuracy in the presence of large host vehicle jitter. The telescope also has high agility--it is capable of a 50-degree retarget (in both axes simultaneously) in less than 2 seconds. To achieve the design specifications, high-accuracy, high-resolution, two-speed resolvers were used, resulting in gimbal-angle measurements stable to 1.5 microradians. In addition, on-axis inertial angle displacement sensors were mounted on the telescope to provide host-vehicle jitter cancellation. The inertial angle sensors are accurate to about 100 nanoradians, but do not measure low frequency displacements below 2 Hz. The gimbal command signal includes host-vehicle attitude information, which is band-limited. This provides jitter data below 20 Hz, but includes a variable latency between 15 and 25 milliseconds. One of the most challenging aspects of this design was to combine the inertial-angle-sensor data with the less perfect information in the command signal to achieve maximum jitter reduction. The optimum blending of these two signals, along with the feedback compensation were designed using Quantitative Feedback Theory.
The purpose of this research was to develop a science-based understanding of the early-time behavior of electric surface arcing in air at atmospheric pressure. As a first step towards accomplishing this, we used a kinetic approach to model an electron swarm as it evolved in a neutral gas under the influence of an applied electric field. A computer code was written in which pseudo-particles, each representing some number of electrons, were accelerated by an electric field. The electric field due to the charged particles was calculated efficiently using a tree algorithm. Collision of the electrons with the background gas led to the creation of new particles through the processes of ionization and photoionization. These processes were accounted for using measured cross-section data and Monte Carlo methods. A dielectric half-space was modeled by imaging the charges in its surface. Secondary electron emission from the surface, resulting in surface charging, was also calculated. Simulation results show the characteristics of a streamer in three dimensions. A numerical instability was encountered before the streamer matured to form branching.
Prediction of the evolution of microstructures in weapons systems is critical to meeting the objectives of stockpile stewardship in accordance with the Nuclear Weapons Test Ban Treaty. For example, accurate simulation of microstructural evolution in solder joints, cermets, PZT power generators, etc. is necessary for predicting the performance, aging, and reliability both of individual components and of entire weapons systems. A recently developed but promising approach called the ''Phase-Field Model'' (PFM) has the potential of allowing the accurate quantitative prediction of microstructural evolution, with all the spatial and thermodynamic complexity of a real microstructure. Simulating with the PFM requires solving a set of coupled nonlinear differential equations, one for each material variable (e.g., grain orientation, phase, composition, stresses, anisotropy, etc.). While the PFM is versatile and is able to incorporate the necessary complexity for modeling real material systems, it is very computationally intensive, and it has been a difficult and major challenge to formulate an efficient algorithmic implementation of the approach. We found that second order in space algorithm is more stable and leads to more accurate results. However, the computational requirements still remain high, so we have developed a single field algorithm to reduce the computations by 2 orders of magnitude. We have created a 3-D parallel version of the basic phase-field (PF model) and benchmarked it performance. Preliminary results indicate that we will be able to run very large problems effectively with the new parallel code. Microstructural evolution in a diffusion couple was simulated using PFM to simultaneously simulate grain growth, diffusion and phase transformation. Solute drag in a variable composition material, a process no other model can simulate, was successfully simulated using the phase-field model. The phase field model was used to study the evolution of fractal high curvature structures to show that these structures have very different morphological and kinetic behaviors than those of equi-axed structures.
Roll-isolated inertial measurement units are developed at Sandia for use in the instrumentation, guidance, and control of rapidly spinning vehicles. Roll-isolation is accomplished by supporting the inertial instrument cluster (gyros and accelerometers) on a single gimbal, the axis of which is parallel to the vehicle's spin axis. A rotary motor on the gimbal is driven by a servo loop to null the roll gyro output, thus inertially stabilizing the gimbal and instrument cluster while the vehicle spins around it. Roll-isolation prevents saturation of the roll gyro by the high vehicle spin rate, and vastly reduces measurement errors arising from gyro scale factor and alignment uncertainties. Nine versions of Sandia-developed roll-isolated inertial measurement units have been flown on a total of 27 flight tests since 1972.
In order to provide real-time data for validation of three dimensional numerical simulations of heterogeneous materials subjected to impact loading, an optically recording velocity interferometer system (ORVIS) has been adapted to a line-imaging instrument capable of generating precise mesoscopic scale measurements of spatially resolved velocity variations during dynamic deformation. Combining independently variable target magnification and interferometer fringe spacing, this instrument can probe a velocity field along line segments up to 15 mm in length. In high magnification operation, spatial resolution better than 10 {micro}m can be achieved. For events appropriate to short recording times, streak camera recording can provide temporal resolution better than 0.2 ns. A robust method for extracting spatially resolved velocity-time profiles from streak camera image data has been developed and incorporated into a computer program that utilizes a standard VISAR analysis platform. The use of line-imaging ORVIS to obtain measurements of the mesoscopic scale dynamic response of shocked samples has been demonstrated on several different classes of heterogeneous materials. Studies have focused on pressed, granular sugar as a simulant material for the widely used explosive HMX. For low-density (65% theoretical maximum density) pressings of sugar, material response has been investigated as a function of both impact velocity and changes in particle size distribution. The experimental results provide a consistent picture of the dispersive nature of the wave transmitted through these samples and reveal both transverse and longitudinal wave structures on mesoscopic scales. This observed behavior is consistent with the highly structured mesoscopic response predicted by 3-D simulations. Preliminary line-imaging ORVIS measurements on HMX as well as other heterogeneous materials such as foam and glass-reinforced polyester are also discussed.
VxInsight provides a visual mechanism for browsing, exploring and retrieving information from a database. The graphical display conveys information about the relationship between objects in several ways and on multiple scales. In this way, individual objects are always observed within a larger context. For example, consider a database consisting of a set of scientific papers. Imagine that the papers have been organized in a two dimensional geometry so that related papers are located close to each other. Now construct a landscape where the altitude reflects the local density of papers. Papers on physics will form a mountain range, and a different range will stand over the biological papers. In between will be research reports from biophysics and other bridging disciplines. Now, imagine exploring these mountains. If we zoom in closer, the physics mountains will resolve into a set of sub-disciplines. Eventually, by zooming in far enough, the individual papers become visible. By pointing and clicking you can learn more about papers of interest or retrieve their full text. Although physical proximity conveys a great deal of information about the relationship between documents, you can also see which papers reference which others, by drawing lines between the citing and cited papers. For even more information, you can choose to highlight papers by a particular researcher or a particular institution, or show the accumulation of papers through time, watching some disciplines explode and other stagnate. VxInsight is a general purpose tool, which enables this kind of interaction with wide variety of relational data: documents, patents, web pages, and financial transactions are just a few examples. The tool allows users to interactively browse, explore and retrieve information from the database in an intuitive way.
Proceedings of the 9th United Engineering Foundation Conference on Risk-Based Decisionmaking in Water Resources - Risk-Based Decisionmaking in Water Resources IX
Data mining involves the discovery and fusion of features from large databases to establish minimal probability of error (MPE) decision and estimation models. Our approach combines a weighted nearest neighbor (WNN) decision model for classification and estimation with genetic algorithms (GA) for feature discovery and model optimization. The WNN model is used to provide a mathematical framework for adaptively discovering and fusing features into near-MPE decision algorithms. The GA is used to discover weighted features and select decision points for the WNN decision model to achieve near-MPE decisions. The performance of the WNN fusion model is demonstrated on the first of two very different problems to demonstrate its robust and practical application to a wide variety of data-mining problems. The first problem involves the isolation of factors that cause hepatitis C virus (HCV) and requires the evaluation of large databases to establish the critical features that can detect with minimal error whether a person is at risk of having HCV. This requires discovering and extracting relevant information (features) from a questionnaire database and combining (fusing) them to achieve a minimal error decision rule. The primary objective of the research is to develop a practical basis for fusing information from questionnaires administered at hospitals to identify and verify features important to isolate risk factors for HCV. The basic problem involves creating a feature database from the questionnaire information, discovering features that provide sufficient information to reliably identify when a person is at risk under conditions with uncertainties caused by recording errors and evasive tactics of people answering the questionnaire. The results of this study demonstrate the WNN fusion algorithm ability to perform in supervised learning environments. The second phase of the research project is directed at the unsupervised learning environment. In this environment the feature data is presented without any classification. Clustering algorithms are developed to partition the feature data into clusters based upon similarity measure models. After the feature data is clustered and classified the supervised WNN fusion algorithms are used to classify the data based upon the minimal probability of error decision rule.
This report provides a summary of the LDRD project titled: An Electromagnetic Imaging System for Environmental Site Reconnaissance. The major initial challenge of this LDRD was to develop a ground penetrating radar (GPR) whose peak and average radiated power surpassed that of any other in existence. Goals were set to use such a system to detect the following: (1) disrupted soil layers where there is potential for buried waste, (2) buried objects such as 55-gallon drums at depths up to 3 m, and (3) detecting contaminated soil. Initial modeling of the problem suggested that for soil conditions similar to Puerto Rican clay loam, moisture content 10 percent (conductivity = 0.01 mhos at 350 MHz), a buried 55-gallon drum could be detected in a straightforward manner by an UWB GPR system at a depth of 3 meters. From the simulations, the highest attenuation ({minus}50 dB) was the result of scattering from a 3-m deep vertically orientated drum. A system loss of {minus}100 dB is a typical limit for all kinds of radar systems (either direct time-domain or swept frequency). The modeling work also determined that the waveshape of the pulse scattered off the buried drum would be relatively insensitive to drum orientation, and thus easier to detect with the GPR system.
Over the past several years, extensive databases have been developed for the S-N behavior of various mate&& used in wind turbine blades, primarily fiberglass composites. These data are typically presented both in their "raw" form and curve fit to define their average properties. For design, confidence limits must be placed on these descriptions. In particular, most designs call for the "95195" design values; namely, with a 95 percent level of confidence, thedesiguerisassuredthat95percentofthematerial will 'meet or exceed the design value. For such material properties as the ultimate streng& the procedures for estimating its value at a particular confidence level is wellffiedifthemeasured values follow a normal or a log-normal distribution. Namely, based upon the number of sample points and their standard deviation, a commonly-found table may be used to determine the survival percentagea t a particular confidencel evel with respect to its mean value. The same is true for fatigue data at a constaut stress level (the number of cycles to failure N at stress level SI). However, when the stress level is allowed to vary, as with a typical S-N fatigue curve, the proceduresf or determmingc onfidencel imits are not as well delked. This paper outlines techn.iques for determimng confklence limits of fatigue data Different approachesto estimating the 95195l evel are compared. Data from the MSUIDOE and the FACT fatigue databam are used to illustrate typical results.
A MEMS test structure capable of measuring friction between polysilicon surfaces under a variety of test conditions has been refined from previous designs. The device is applied here to measuring friction coefficients of polysilicon surfaces under different environmental, loading, and surface conditions. Two methods for qualitatively comparing friction coefficients (μ) using the device are presented. Samples that have been coated with a self-assembled monolayer of the lubricating film perfluorinated-decyltrichlorosilane (PFTS) have a coefficient of friction that is approximately one-half that of samples dried using super-critical CO2 (SCCO2) drying. Qualitative results indicate that μ is independent of normal pressure. Wear is shown to increase μ for both supercritically dried samples and PFTS coated samples, though the mechanisms appear to be different. Super critically dried surfaces appear to degrade continuously with increased wear cycles, while PFTS coated samples reach a steady state friction value after about 105 cycles.
Hydrogen atom-hydrogen atom scattering is a prototype for many of the fundamental principles of atomic collisions. In this work we present the formalism and the predictions of a time-dependent self-consistent-field description of the H + H system for scattering in the intermediate energy regime of 1-100 keV. Because of the unrestricted nature of the numerical orbital description, this method includes the effects of an unlimited basis set within each orbital. Electron exchange and a limited amount of electron correlation are included as well. We solve numerically coupled three-dimensional Schrodinger equations for the two-electron orbitals in singlet and triplet symmetries. Excitation and ionization cross sections are computed and compared with other theory and experiment. The results capture many features of the problem but illustrate a need for more quantitative experimental information concerning the H + H system in this energy range.
The cookoff of energetic materials involves the combined effects of several physical and chemical processes. These processes include heat transfer, chemical decomposition, and mechanical response. The interaction and coupling between these processes influence both the time-to-event and the violence of reaction. The prediction of the behavior of explosives during cookoff, particularly with respect to reaction violence, is a challenging task. To this end, a joint DoD/DOE program has been initiated to develop models for cookoff, and to perform experiments to validate those models. In this paper, a series of cookoff analyses are presented and compared with data from a number of experiments for the aluminized, RDX-based, Navy explosive PBXN-109. The traditional thermal-chemical analysis is used to calculate time-to-event and characterize the heat transfer and boundary conditions. A reaction mechanism based on Tarver and McGuire's work on RDX{sup 2} was adjusted to match the spherical one-dimensional time-to-explosion data. The predicted time-to-event using this reaction mechanism compares favorably with the validation tests. Coupled thermal-chemical-mechanical analysis is used to calculate the mechanical response of the confinement and the energetic material state prior to ignition. The predicted state of the material includes the temperature, stress-field, porosity, and extent of reaction. There is little experimental data for comparison to these calculations. The hoop strain in the confining steel tube gives an estimation of the radial stress in the explosive. The inferred pressure from the measured hoop strain and calculated radial stress agree qualitatively. However, validation of the mechanical response model and the chemical reaction mechanism requires more data. A post-ignition burn dynamics model was applied to calculate the confinement dynamics. The burn dynamics calculations suffer from a lack of characterization of the confinement for the flaw-dominated failure mode experienced in the tests. High-pressure burning rates are needed for more detailed post-ignition studies. Sub-models for chemistry, mechanical response and burn dynamics need to be validated against data from less complex experiments. The sub-models can then be used in integrated analysis for comparison with experimental data taken during integrated tests.
Since the discovery of surfactant-templated silica by Mobil scientists in 1992, mesostructured silica has been synthesized in various forms including thin films, powders, particles, and fibers. In general, mesostructured silica has potential applications, such as in separation, catalysis, sensors, and fluidic microsystems. In respect to these potential applications, mesostructured silica in the form of thin films is perhaps one of the most promising candidates. The preparation of mesostructured silica films through preferential solvent evaporation-induced self-assembly (EISA) has recently received much attention in the laboratories. However, no amphiphile/silica films with reverse mesophases have ever been made through this EISA procedure. Furthermore, templates employed to date have been either surfactants or poly(ethylene oxide)-b-poly(propylene oxide)-b-poly(ethylene oxide) triblock copolymers, such as pluronic P-123, both of which are water-soluble and alcohol-soluble. Due to their relatively low molecular weight, the templated silica films with mesoscopic order have been limited to relatively small characteristic length scales. In the present communication, the authors report a novel synthetic method to prepare mesostructured amphiphilic/silica films with regular and reverse mesophases of large characteristic length scales. This method involves evaporation-induced self-assembly (EISA) of amphiphilic polystyrene-block-poly(ethylene oxide) (PS-b-PEO) diblock copolymers. In the present study, the PS-b-PEO diblocks are denoted as, for example, PS(215)-b-PEO(100), showing that this particular sample contains 215 S repeat units and 100 EO repeat units. This PS(215)-b-PEO(100) diblock possesses high molecular weight and does not directly mix with water or alcohol. To the authors knowledge, no studies have reported the use of water-insoluble and alcohol-insoluble amphiphilic diblocks as structure-directing agents in the synthesis of mesostructured silica films through EISA. It is believed that the present system is the first to yield amphiphile/silica films with regular and reverse mesophases, as well as curved multi-bilayer mesostructures, through EISA. The ready formation of the diblock/silica films with multi-bilayer vesicular mesostructures is discussed.
Heat pipes are often proposed as cooling system components for small fission reactors. SAFE-300 and STAR-C are two reactor concepts that use heat pipes as an integral part of the cooling system. Heat pipes have been used in reactors to cool components within radiation tests (Deverall, 1973); however, no reactor has been built or tested that uses heat pipes solely as the primary cooling system. Heat pipe cooled reactors will likely require the development of a test reactor to determine the main differences in operational behavior from forced cooled reactors. The purpose of this paper is to describe the results of a systems code capable of modeling the coupling between the reactor kinetics and heat pipe controlled heat transport. Heat transport in heat pipe reactors is complex and highly system dependent. Nevertheless, in general terms it relies on heat flowing from the fuel pins through the heat pipe, to the heat exchanger, and then ultimately into the power conversion system and heat sink. A system model is described that is capable of modeling coupled reactor kinetics phenomena, heat transfer dynamics within the fuel pins, and the transient behavior of heat pipes (including the melting of the working fluid). The paper focuses primarily on the coupling effects caused by reactor feedback and compares the observations with forced cooled reactors. A number of reactor startup transients have been modeled, and issues such as power peaking, and power-to-flow mismatches, and loading transients were examined, including the possibility of heat flow from the heat exchanger back into the reactor. This system model is envisioned as a tool to be used for screening various heat pipe cooled reactor concepts, for designing and developing test facility requirements, for use in safety evaluations, and for developing test criteria for in-pile and out-of-pile test facilities.
An integrated approach, combining the continuum theory of sintering and Potts model based mesostructure evolution analysis, is used to solve the problem of bi-layered structure sintering. Two types of bi-layered structures are considered: layers of the same material with different initial porosity, and layers of two different materials. The effective sintering stress for the bi-layer powder sintering is derived, both at the meso- and the macroscopic levels. Macroscopic shape distortions and spatial distributions of porosity are determined as functions of the dimensionless specific time of sintering. The effect of the thickness of the layers on shrinkage, warpage, and pore-grain structure is studied. Ceramic ZnO powders are employed as a model experimental system to assess the model predictions.
This paper presents an analysis technique used to generate the structural response at locations not measured during the ejection of a captive-carried store. The ejection shock event is complicated by the fact that forces may be imparted to the store at eight distinct locations. The technique derives forcing functions by combining the initial field test data for a limited number of measurement locations with Frequency Response Functions (FRFs) measured using a traditional modal-type impact (tap) test at the same locations. The derived forcing functions were then used with tap test FRFs measured at additional locations of interest to produce the desired response data.
Operating in a potentially hostile city is every soldier's nightmare. The staggering complexity of the urban environment means that deadly threats--or non-combatants-may lurk behind every corner, doorway, or window. Urban operations present an almost unparalleled challenge to the modern professional military. The complexity of urban operations is further amplified by the diversity of missions that the military will be called upon to conduct in urban terrain. Peace-making and peace-keeping missions, urban raids to seize airports or WMD sites or to rescue hostages, and extended urban combat operations all present different sorts of challenges for planners and troops on the ground. Technology almost never serves as a magic bullet, and past predictions of technological miracles pile high on the ash heap of history. At the same time, it is a vital element of planning in the modern age to consider and, if possible, take advantage of emerging technologies. We believe that technologies can assist military operations in urbanized terrain (MOUT) in three primary areas, which are discussed.
This paper summarizes an evaluation of a treatment loop designed to upgrade the quality of spent rinse waters discharged from 10 wet benches located in the fab at Sandia's Microelectronics Development Laboratory (MDL). The goal of the treatment loop is to make these waters, presently being discharged to the fab's acid waste neutralization (AWN) station, suitable for recycling as feed water back into the fab's ultrapure water (UPW) plant. The MDL typically operates 2 shifts per day, 5 days per week. Without any treatment, the properties of the spent rinse waters now being collected have been shown to be compatible with recycling about 30% (50/168) of the time (weekends primarily, when the fab is idling) which corresponds to about 12% of the present water discharged from the fab to the AWN. The primary goal of adding a treatment loop is to increase the percentage of recyclable water from these 10 wet benches to near 100%, increasing the percentage of total recyclable water to near 40% of the total present fab discharge to the AWN. A second goal is to demonstrate compatibility with recycling this treated spent rinse water to the present R/O product water tank, reducing both the present volume of R/O reject water and the present load on the R/O. The approach taken to demonstrate achieving these goals is to compare all the common metrics of water quality for the treated spent rinse waters with those of the present R/O product water. Showing that the treated rinse water is equal or superior in quality to the water presently stored in the R/O tank by every metric all the time is assumed to be sufficient argument for proceeding with plans to incorporate recycling of these spent rinse waters back into MDL's R/O tank.
The use of plasma spray to deposit thin metal-sulfide cathode films is described in this paper. Conventional electroactive stack components in thermal batteries are constructed from pressed-powder parts that are difficult to fabricate in large diameters in thicknesses <0.010. Plasma-sprayed electrodes do not steer from this difficulty, allowing greater energy densities and specific energies to be realized. Various co-spraying agents have been found suitable for improving the mechanical as well as electrochemical properties of plasma-sprayed cathodes for thermal batteries. These electrodes generally show equal or improved performance over conventional pressed-powder electrodes. A number of areas for future growth and development of plasma-spray technology is discussed.
When embarking on an experimental program for purposes of discovery and understanding, it is only prudent to use appropriate analysis tools to aid in the discovery process. Due to the limited scope of experimental measurement analytical results can significantly complement the data after a reasonable validation process has occurred. In this manner the analytical results can help to explain certain measurements, suggest other measurements to take and point to possible modifications to the experimental apparatus. For these reasons it was decided to create a detailed nonlinear finite element model of the Sandia Microslip Experiment. This experiment was designed to investigate energy dissipation due to microslip in bolted joints and to identify the critical parameters involved. In an attempt to limit the microslip to a single interface a complicated system of rollers and cables was devised to clamp the two slipping members together with a prescribed normal load without using a bolt. An oscillatory tangential load is supplied via a shaker. The finite element model includes the clamping device in addition to the sequence of steps taken in setting up the experiment. The interface is modeled using Coulomb friction requiring a modest validation procedure for estimating the coefficient of friction. Analysis results have indicated misalignment problems in the experimental procedure, identified transducer locations for more accurate measurements, predicted complex interface motions including the potential for galling, identified regions where microslip occurs and during which parts of the loading cycle it occurs, all this in addition to the energy dissipated per cycle. A number of these predictions have been experimentally corroborated in varying degrees and are presented in the paper along with the details of the finite element model.
In this paper the authors investigate the relationship between glassy and ferromagnetic phases in disordered Ising ferromagnets in the presence of transverse magnetic fields, {Lambda}. Iterative mean field simulations probe the free energy landscape and suggest the existence of a glass transition line in the {Lambda}, temperature T plane well within the ferromagnetic phase. New experimental field-cooled and zero-field-cooled data on LiHo{sub x} Y{sub 1{minus}x}F{sub 4} provide support for our theoretical picture.
As a joint is loaded, the tangent stiffness of the joint reduces due to slip at interfaces. This stiffness reduction continues until the direction of the applied load is reversed or the total interface slips. Total interface slippage in joints is called macro-slip. For joints not undergoing macro-slip, when load reversal occurs the tangent stiffness immediately rebounds to its maximum value. This occurs due to stiction effects at the interface. Thus, for periodic loads, a softening and rebound hardening cycle is produced which defines a hysteretic, energy absorbing trajectory. For many jointed sub-structures, this hysteretic trajectory can be approximated using simple polynomial representations. This allows for complex joint substructures to be represented using simple non-linear models. In this paper a simple one dimensional model is discussed.
Unattended monitoring systems are being studied as a means of reducing both the cost and intrusiveness of present nuclear safeguards approaches. Such systems present the classic information overload problem to anyone trying to interpret the resulting data not only because of the sheer quantity of data but also because of the problems inherent in trying to correlate information from more than one source. As a consequence, analysis efforts to date have mostly concentrated on checking thresholds or diagnosing failures. Clearly more sophisticated analysis techniques are required to enable automated verification of expected activities level concepts in order to make automated judgments about safety, sensor system integrity, sensor data quality, diversion, and accountancy.
Multiple techniques have been developed to model the temporal evolution of infectious diseases. Some of these techniques have also been adapted to model the spatial evolution of the disease. This report examines the application of one such technique, the SEIR model, to the spatial and temporal evolution of disease. Applications of the SEIR model are reviewed briefly and an adaptation to the traditional SEIR model is presented. This adaptation allows for modeling the spatial evolution of the disease stages at the individual level. The transmission of the disease between individuals is modeled explicitly through the use of exposure likelihood functions rather than the global transmission rate applied to populations in the traditional implementation of the SEIR model. These adaptations allow for the consideration of spatially variable (heterogeneous) susceptibility and immunity within the population. The adaptations also allow for modeling both contagious and non-contagious diseases. The results of a number of numerical experiments to explore the effect of model parameters on the spread of an example disease are presented.
This Safety Analysis Report (SAR) is prepared in compliance with the requirements of DOE Order 5480.23, Nuclear Safety Analysis Reports, and has been written to the format and content guide of DOE-STD-3009-94 Preparation Guide for U. S. Department of Energy Nonreactor Nuclear Safety Analysis Reports. The Hot Cell Facility is a Hazard Category 2 nonreactor nuclear facility, and is operated by Sandia National Laboratories for the Department of Energy. This SAR provides a description of the HCF and its operations, an assessment of the hazards and potential accidents which may occur in the facility. The potential consequences and likelihood of these accidents are analyzed and described. Using the process and criteria described in DOE-STD-3009-94, safety-related structures, systems and components are identified, and the important safety functions of each SSC are described. Additionally, information which describes the safety management programs at SNL are described in ancillary chapters of the SAR.
The magnetic implosion of a high-Z quasi-spherical shell filled with DT fuel by the 20-MA Z accelerator can heat the fuel to near-ignition temperature. The attainable implosion velocity on Z, 13-cm/{micro}s, is fast enough that thermal losses from the fuel to the shell are small. The high-Z shell traps radiation losses from the fuel, and the fuel reaches a high enough density to reabsorb the trapped radiation. The implosion is then nearly adiabatic. In this case the temperature of the fuel increases as the square of the convergence. The initial temperature of the fuel is set by the heating of an ion acoustic wave to be about 200-eV after a convergence of 4. To reach the ignition temperature of 5-keV an additional convergence of 5 is required. The implosion dynamics of the quasi-spherical implosion is modeled with the 2-D radiation hydrodynamic code LASNEX. LASNEX shows an 8-mm diameter quasi-spherical tungsten shell on Z driving 6-atmospheres of DT fuel nearly to ignition at 3.5-keV with a convergence of 20. The convergence is limited by mass flow along the surface of the quasi-spherical shell. With a convergence of 20 the final spot size is 400-{micro}m in diameter.
A theory is developed for the response of moderately porous solids (no more than {approximately}20% void space) to high-strain-rate deformations. The model is consistent because each feature is incorporated in a manner that is mathematically compatible with the other features. Unlike simple p-{alpha} models, the onset of pore collapse depends on the amount of shear present. The user-specifiable yield function depends on pressure, effective shear stress, and porosity. The elastic part of the strain rate is linearly related to the stress rate, with nonlinear corrections from changes in the elastic moduli due to pore collapse. Plastically incompressible flow of the matrix material allows pore collapse and an associated macroscopic plastic volume change. The plastic strain rate due to pore collapse/growth is taken normal to the yield surface. If phase transformation and/or pore nucleation are simultaneously occurring, the inelastic strain rate will be non-normal to the yield surface. To permit hardening, the yield stress of matrix material is treated as an internal state variable. Changes in porosity and matrix yield stress naturally cause the yield surface to evolve. The stress, porosity, and all other state variables vary in a consistent manner so that the stress remains on the yield surface throughout any quasistatic interval of plastic deformation. Dynamic loading allows the stress to exceed the yield surface via an overstress ordinary differential equation that is solved in closed form for better numerical accuracy. The part of the stress rate that causes no plastic work (i.e-, the part that has a zero inner product with the stress deviator and the identity tensor) is given by the projection of the elastic stressrate orthogonal to the span of the stress deviator and the identity tensor.The model, which has been numerically implemented in MIG format, has been exercised under a wide array of extremal loading and unloading paths. As will be discussed in a companion sequel report, the CKP model is capable of closely matching plate impact measurements for porous materials.